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Summary 

We  have  investigated  many  topics  in  the  theory  of  semiconductors  and 
alloys  including  (1)  Generalized  Brooks'  formula  and  the  electron  mobility  in 
SiGe  alloys,  (2)  bond-length  distribution,  lattice  relaxation,  bond  energies 
and  mixing  enthalpies  in  alloys,  (3)  hardness  and  dislocation  energy,  (4)  SiGe 
alloys  -  band  structures  and  core-exciton  binding  energy  and  linewidth,  (5) 
alloy  statistics  and  raicroclustering  and  (6)  sensitivity  of  deep  levels  to 
band  structure  and  potential.  The  important  results  are  outlined  below.  The 
details  are  discussed  in  eight  publications  enclosed  with  this  report. 

Generalized  Brooks'  Formula  and  the 
Electron  Mobility  in  SiGe  Alloy 

Although  Brooks'  formula  has  been  widely  used  for  calculating  the  alloy¬ 
scattering  limited  electron  mobility,  we  showed  that  this  formula  is  only 
valid  for  a  direct-gap  semiconductor.  Besides,  there  are  questions  about  the 
scattering  parameter  and  the  effective  mass.  We  generalized  the  formula  for 
indirect-gap  alloys  with  multiple  bands  and  applied  it  to  SiGe  alloy.  Our 
results,  correlated  well  with  experiments,  showed  that  the  electron  mobility 
drops  fast  with  alloying.  The  mobility  has  a  dip  at  15%  Si  concentration, 
corresponding  to  a  transition  from  the  X(A )  band  edge  to  the  L  edge. 

Bond  Lengths,  Lattice  Relaxation,  and  Mixing  Enthalpies 
in  Semiconductor  Alloys 

We  treated  the  problem  with  a  model  which  combines  Harrison's  bonding 
theory  with  the  valence  force  field  model  and  an  elastic  continuum.  While  the 
local  strain  is  the  main  driving  force  for  the  bimodal  bond-length  distribu¬ 
tion  in  pseudo-binary  alloys  found  in  the  EXAFS  experiment,  we  found  that  the 
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chemical  shifts  arising  from  different  bond  lengths  and  polarities  of  the 
constituent  bonds  can  have  a  significant  contribution  to  the  mixing  enthalpies, 
In  fact,  the  sizable  negative  values  of  the  chemical  shifts  in  the  cation- 
substitutional  alloys,  e.g.  Ga^Inj-xAs,  may  be  important  for  stabilizing  the 
mixture.  We  also  deduced  a  simple  criterion  for  separating  miscible  from 


immiscible  alloys. 


Hardness  and  Dislocation  Energy 


We  showed  that,  if  the  dislocation  energy  in  semiconductor  is  dominated 
by  interactions  among  dislocations,  then  the  hardness  is  an  intrinsic  material 
property,  independent  of  the  size  of  the  indentation  and  the  force  applied. 

The  fact  that  dislocation  energy  is  proportional  to  hardness  suggests  that 
alloy  hardening  is  a  mechanism  for  reducing  dislocation  densities. 


SiGe  Alloys  -  Band  Structure  and  Core-Excitons 


SiGe  binary  alloy  has  regenerated  research  interests  because  of  its 
potential  application  for  high-mobility  devices  in  the  strained  superlattice 
configuration.  We  have  applied  our  techniques  to  obtain  high-quality  band 
structures  for  Si  and  Ge  and  have  performed  alloy  calculations  for  Si^Ge]_j{. 
The  results  have  been  checked  against  optical  spectra  and  used  for  mobility 
study  mentioned  earlier.  Another  interesting  result  is  that  the  alloy  band 
parameters  allowed  us  to  correlate  the  Si  2p  core-exciton  binding  energy  with 
the  linewidth.  The  observed  minimum  in  the  linewidth  near  X  =  0.15  can  be 
explained  as  the  result  of  a  competetion  between  intrinsic  broadening  due  to 
screening  and  extrinsic  alloy  broadening.  The  most  reasonable  binding  energy 
in  pure  Si  was  found  to  be  0.15  eV.  This  work  thus  helps  resolve  the  contro¬ 
versy  about  the  unusually  large  binding  energy  for  the  core  exciton  in  Si. 


2 


Statistics  and  Micro-clustering  in  Semiconductor  Alloys 


One  of  our  major  efforts  was  to  develop  a  statistical  theory  for  semi¬ 
conductor  alloys  which  covers  three  aspects  -energetics,  statistics  and  phase 
diagram.  Because  of  the  smallness  of  the  mixing  enthalpies,  typically  several 
k  Cal/mol,  the  present  first-principle  theory  is  not  accurate  enough  for  this 
purpose.  We  found  that  a  combination  of  Harrison's  model  and  the  valence 
force  field  model  provides  a  simple  and  adequate  method  for  calculating  the 
mixing  energies.  We  then  generalized  Guggenheim's  quasi-chemical  aoproxima- 
tion  to  treat  the  tetrahedral  clusters.  Our  results  showed  that  the  arrange¬ 
ment  of  atoms  can  have  a  appreciable  deviation  from  the  random  distribution. 
The  distribution  is  governed  by  the  mixing  energy  which  is  the  sum  of  the 
strain  energy  and  chemical  shifts  mentioned  earlier.  The  non-random  distribu¬ 
tion  will  have  profound  effects  on  the  band  structures,  mobilities,  mechanical 
properties  and  growth  mechanisms,  and  will  reflect  in  many  measurable 
properties  such  as  phonon  spectra,  EXAFS,  deep  level  spectra  etc.  The  study 
of  these  effects  is  one  of  our  current  activities. 
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A  formula  for  alloy-scattering-limited  electron  mobility  in  semiconductors  is  obtained  for 
indirect  gap  systems  with  multiple  band  minima.  All  the  input  parameters  needed  are  defined 
explicitly.  The  drift  mobility  of  Si,  Ge ,  _  ,  which  has  a  dip  at  x  —  0. 1 3  and  a  broader  minimum  at 
XT  ~  0.5  is  calculated  by  adding  alloy  scattering  toother  scattering  mechanisms  and  correlates  well 
with  the  measured  Hail  mobility. 


The  electron  and  hole  mobilities  in  semiconductors  are 
determined  by  the  band  structure  and  various  scattering 
mechanisms,  predominately  impurity  and  phonon  scatter¬ 
ing.  For  alloys,  the  mobility  is  also  affected  by  disorder  aris¬ 
ing  from  aperiodic  atomic  potentials  and  atomic  positions. 
Many  years  ago.  Nordheim'  and  Brooks’  obtained  an 
expression  for  ailoy-scattering-limited  electron  mobilities  in 
metals  and  semiconductors,  respectively.  Brooks*  well- 
known  formula  reads 

1 

P*  ~~  ;  ;  ijjT  - > 

3jc(1  {AE]\KT 

where  .V,,  is  the  number  of  atoms  per  unit  volume,  m*  is  a 
band-edge  effective  mass,  x  is  the  fractional  concentration  of 
one  of  the  species,  and  Jf*  is  an  energy  parameter  character¬ 
izing  the  alloy  potential  fluctuations.  Although  this  formula 
has  been  widely  and.  to  some  extent,  successfully  used  for 
direct  gap  matenals.  '"’  the  identification  of  the  alloy  disor¬ 
der  parameter  AE  remains  uncertain.  Various  suggestions 
have  previously  been  made  for  J£,  e.g.,  and  band-edge  dis¬ 
continuity"  or  band-gap  differences.’  Any  of  these  simple 
choices  is  bound  to  fail  when  one  applies  Eq.  (11  to  more 
complicated  indirect  gap  systems  such  as  Si,  Ge,  _  ,  alloys, 
w  here  one  encounters  conduction-band  minima  transferring 
between  the.V  and  L  points  of  the  Brillouin  zone.  For  exam¬ 
ple.  if  Jf  is  taken  to  be  the  difference  in  corresponding  band 
edges,  then  one  finds  that  AE  ~0. 1  eV  for  the  X  lA  )  valley 
and  —  1.2  eV  for  the  L  valley.  The  values  that  fit  the  experi¬ 
ment  are  about  half  this  value  for  £  and  —0.5  eV  for  X"  The 
purpose  of  this  letter  is  to  resolve  the  identity  of  for  indi¬ 
rect  gap  materials. 

Moreover,  there  is  a  problem  with  the  m*  that  enters 
Eq.  ( 1 1.  For  direct  gap  alloys,  the  band-edge  effective  mass  at 
r  naturally  enters  Eq.  |1|.  For  the  indirect  gap  alloys,  the 
effective  mass  is  anisotropic  and  hence  an  appropriate  mass 
must  be  chosen.  Previous  authors"  ’  have  chosen  m*  to  be 
the  effective  conductivity  mass  m*.  We  shall  show  that  dif¬ 
ferent  masses  enter  for  different  cases. 

The  first  unambiguous  assignment  for  J£in  a  direct  gap 
alloy  was  given  by  Hass  etal '  To  estimate  the  limiting  elec¬ 
tron  mobility  in  Hg,  Cd,  ,  Te  based  on  a  tight-binding  (TD) 
band  description,  they  defined  AE  to  be  f.AE,,  where  f  is 
the  s  fraction  in  the  density  of  states  and  J£  is  the  difference 
between  the  t  atomic  term  values  of  the  Hg  and  Cd  atoms.  By- 
extending  this  approach  to  alloys  with  indirect  gaps  and 


multiple  bands,  we  show  that  all  the  uncertainties  identified 
above  are  resolved.  Our  generalized  Brooks’  formula  will 
then  be  applied  to  Si,  Ge,  systems  to  explain  their  ob¬ 
served  mobility.’* 

Because  Brooks’  formula  has  never  been  derived  expli¬ 
citly  in  the  literature,  we  rederive  it  first  and  then  generalize 
it.  Consider  the  case  of  a  single  band  with  an  isotropic  effec¬ 
tive  mass.  The  dc  electronic  conductivity  based  on  the  linear 
response  theory’  is  given  by 

where  the  energy-dependent  (7\6)  in  the  weak  alloy  scattering 
limit  is 

off)  =  (e’/3)ii’(f  IZ)  (f  )r1e).  (3) 

D  (ej  is  the  density  of  states  (DOS)  per  unit  volume  for  both 
spins,  so  D  If)  =  IN^pie),  with  -V^  =  .Vo/2  being  the  number 
of  unit  cells  per  volume  Ifor  the  diamond  structure,  half  the 
number  of  atoms  .V„  per  unit  volume;  andp(f )  being  the  DOS 
per  unit  cell  per  spin.  The  mean  square  velocity  u’(e)  for 
carriers  with  energy  e  is  given  by 

.;=(e)  =  S[e-eik]] 

k  pif] 

The  scattering  lifetime  for  carriers  with  energy  e,  tie),  is  re¬ 
lated  to  the  alloy  broadening  J  (f  |  by  rtel  =  H/IA  (e),  where 
the  energy  J  (el  is  the  imaginary  part  of  the  self-energy  in  the 
averaged  alloy  Green's  function.  For  weak  scattering  J  (el  is 

J  (f)  =  rr.xi  1  —  x)  (J£  )‘/?(el,  (5) 

where  in  a  tight  binding  (TB)  description  J  £  is  the  difference 
in  the  term  values  of  the  constituents.  Then  the  mobility  is 
p ,  =  cr/ne  with  the  electron  density  given  by 


For  a  nondegenerate  semiconductor.y'le)  is  the  Boltzmann 
distribution  and/|fia  e  '*  Furthermore,  for  a  para¬ 

bolic  bande'kl  =  frk  -/2m*.p{€\  =  (2m*)'''f then 
all  the  above  equations  can  be  combined  to  arrive  at  Eq.  1 1 1. 

For  a  real  semiconductor  alloy  in  a  TB  description,  the 
alloy  scattering  can  be  characterized  by  two  parameters  J 
and  AEj..  the  differences  in  s  and p  term  values  between  two 
substitutional  atoms.  Then  an  effective  broadening  is  given 
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TAULE  I.  CJculalioii  paramcicrs 


Ale)  =  {A,p,  +^^p^]/p,  (7) 

where^,  ,/>,  are  partial  density  of  states  IPDOS)  and  J,  and 
4,  are  similar  to  Eq.  (5),  with p  replaced  by  p^  and  p^,  re¬ 
spectively.  For  Hg,  _,Cd,Te,  the  s  disorder  is  predomin- 
eni'  ‘“  and  one  can  neglect  Defining  p,  =/p  (and 
Pp  ~fpP)>  one  amves  at 

Thus  as  was  pointed  out  by  Hasset  al.,^  f,AE^  plays  the  role 
of  in  this  special  case  where  AE^,  can  be  neglected. 

For  an  alloy  with  a  single  indirect  gap  minimum,  one 
has  to  consider  both  s  and p  contributions  to  the  alloy  broad¬ 
ening  and  the  masses  that  enter  p  and  o'.  Again.  Eqs.  (2H<>I 
can  be  combined  to  yield 


_ _ |e1g^V,.v2lr) _ 

[  3:r(  1  -  ) ' ' -( A-  7- ) " X  ( /r J £ :  +  J £ ; )  ]  ■ 

18) 

where  m*  and  m*  are  respectively  the  longitudinal  and  the 
transverse  mass  at  the  band  edge,  and  is  the  number  of 
equivalent  minima,  e.g.,  6  for  Si.  The  conductivity  mass  m* 
comes  from  averaging  y*  in  Eq.  (2)  and  is  given  by  3(2/mf 
-b  l/mf)"'.  Equation  (8)  clearly  identifies  the  masses  and 
the  energy  parameter  that  enter  Brooks’  formula. 

Ne.xt  we  consider  a  still  more  complicated  case  where 
the  contribution  to  the  mobility  comes  from  more  than  one 
band.  For  e.sample.  in  Si,  Ge,  _ ,  the  ylT  and  L  minima  cross 
near  .r  =  0. 1 5. "  There  are  now  two  contributions  to  the  net 
conductivity,  so  <7  =  la, ,  where  i  is  X  or  L.  The  quantities 
D,  ifi.  and  .V,  jf )  now  take  different  values  for  different 
bands.  The  structure  of  r,  (el  requires  more  careful  consider¬ 
ation.  The  complication  comes  from  the  fact  that  the  effec¬ 
tive  broadening  A  is  still  given  by  Eq.  |7),  but  p,,Pp,  and p 
contain  contributions  from  both  the  bands.  The  proper  ex¬ 
pressions  are p  —  1,p,  S  and p„  =  2,/,,  N  [.p,  where  i  —  X 
or  L,  a  =  s  or p,  and  .V^  =  6,  =  4.  The  equation  for  J  is 

A  {€)  =  irx\l  -  x]'^(^Y/„,X',p,{e]AE„ 

The  mobility  associated  with  the  /th  band  is  defined  as 
p,  =  a,/{n,e),  then 

,  =  - ! -  /  ^  ,,0) 

3x(l-.x)  [m>,*(2mr)''-], 


/,  =  r  - 


de.  (11) 


Thus,  the  generalized  formula  no  longer  has  the  explicit 
rand  F  dependences  of  the  original  Brooks' form.  However, 
all  the  quantities  needed — the  masses,  the  scattering  param¬ 
eters  J£,  .  the  band  gaps,  and  the  fractions/,, — can  be  eval¬ 
uated  theoretically  without  resorting  to  experimentally  fit¬ 
ted  parameters.  To  demonstrate,  we  shall  apply  Eq.  (lOI  to 
Si,  Ge,  , .  The  band  quantities  are  obtained  from  our  CPA 
calculation, ' '  We  found  that  the  etfective  mas.ses  vary  weak- 
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ly  with  the  concentration,  so  m*  and  m*  are  assumed  to  be 
constant  and  assigned  the  values  0.97  and  0.19  for  the  X 
minima  and  1.64  and  0.082  for  the  L  minima,  respectively. 
The  calculated  energy  gaps  for  the  X  (J  )  follows  the  func¬ 
tional  form  E'*  '  =  a  -i-  bx  -t-  ex'  and  for  L  is  given  by  E'^^-' 
=  A  -i-  Bx  +  Cr'.  .411  the  parameters  ofour  calculations  are 
listed  in  Table  I. 

To  correlate  the  calculation  with  the  measured  mobili¬ 
ties.  we  need  to  have  an  estimate  of  scattering  rates  l/r„  due 
to  impurities  and  phonons.  A  crude  approximation  is  to  as¬ 
sume  l/r„  for  a  given  valley  to  be  the  same  as  the  appropriate 
constituent’s  values  and  add  to  it  the  alloy  scattering  rate  1/ 
.  Then  the  average  mobility  and  the  mobility  from  the  ith 
minimum  in  the  alloy  are 

lu,)"'  =  (|t4?)-'  (12) 

p  *  is  given  by  Eq.  (10)  and  p°  are  the  measured  drift  mobili¬ 
ties  for  Si  or  Ge,'*  The  drift  mobility,  calculated  from  Eq. 
(12),  is  plotted  as  a  function  of  alloy  concentrations:  in  Fig.  1. 

For  .r<0.05  and  s:>0.20,  the  energy  difference  between 
the  X  and  L  edges  is  large  enough  so  there  is  a  negligible 
contribution  to  the  mobility  from  the  higher  minima.  In  the 
Si,  Gei  _  ,  system,  the  s  scattering  is  predominent.  Because 


FIG  I  Cjicul.ited  sohd  lmi'‘  electron  dnfl  mobility  and  the  ecpenmental 
Hall  mobility  idached  iiiiei  tVom  Rets  6  and  8  are  plotted  as  a  function  of 
alloy  conccnfralmn 
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the  L  edges  have  a  larger  i  eontent,  alloy  scattenng  reduces 
the  average  mobility  substantially  for  small  x.  Even  though 
the  j  content  is  almost  the  same  for  all  20  at  the-1'  edge, 
the  mobility  still  decreases  to  ;ic  — 0.5  as  shown  in  Fig.  1  be¬ 
cause  of  the  a:(  1  —  x)  term  in  Eq.  1 10). 

An  interesting  feature  is  obtained  for  the  compositions 
0.13sx  iO.  18.  The  average  mobility  attains  a  local  mini¬ 
mum  near  x  =  0.14  and  a  smaller  maximum  at  t  =  0.17. 
This  feature  occurs  because  of  the  X  \o  L  crossover.”  For 
x<0. 14,  the  major  contribution  xofi  comes  from  L  minima. 
Near  x  =  0. 14,  the  density  of  states  increases  because  the  X 
and  L  minima  merge.  So  the  alloy  scattering  increases  there 
and  the  average  mobility  decreases.  For.r>0.14,  the.Ybands 
have  the  lower  minima.  .4s  the  s  content  is  small  at  the  X 
minima,  the  reduced  alloy  scattering  increases  the  average 
mobility.  For  larger  values  of .r,  the  x(  1  —  x)  term  takes  over 
and  the  mobility  varies  as  shown.  The  values  of  measured 
Hall  mobility  in  Si,  Ge,  _  ,  systems  are  also  plotted  in  Fig.  1. 
The  interesting  feature  near  x  =  0. 14  is  clearly  seen.  Since 
the  experimental  drift  mobility  fj.^  for  Si.Ge,  _  ,  is  not 
available  and  the  generalization  of  Eq.  1 10)  to  Hall  mobility 
fj.„  IS  less  clear,  we  present  thecalculated^Q  and  experimen¬ 
tal  iRef.  7,3)  here.  While  we  do  not  expect  quantitative 
agreement,  because  can  range  from  1  to  2, '  we  do 

expect  them  to  display  the  same  qualitative  .t  dependence.  It 
is  rewarding  to  note  the  similarity  in  the  trend  in  Fig.  1. 
Previous  authors  explained  the  dip  in  the  mobility  curve  by 
including  intervalley  scattering  with  an  arbitrary  adjustable 
coupling  constant.'’  Our  calculations  automatically  include 
that  portion  of  intervalley  scattering  that  results  from  alloy 
disorder  with  a  coupling  constant  set  by  the  atomic  proper¬ 


ties  of  the  constituents.  However,  the  additon  of  intervafiey 
scattenng  mediated  by  phonons  and  impurities  is  expected 
to  increase  the  dip  near  crossover. 

In  summary,  an  expression  for  alloy-scattering-limited 
charge  carrier  mobilities  is  derived  for  indirect  gap  alloys 
with  multiple  bands.  This  expression  reduces  to  Brooks’  for¬ 
mula  for  direct  gap  alloys.  The  quantities  m*  and  can  be 
calculated  exactly.  Alloy  scattering  accounts  for  the  ob¬ 
served  mobility  features  in  the  Si,  Ge,  _ ,  alloy,  including 
the  anomaly  near  the  L  to  X  [A]  crossover. 
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FIG.  1.  Variation  of  the  band  gap  (solid  line)  and  the  Si 
2/7  core-exeiton  level  with  .r  in  Si,Gei-,  alloys.  The  energy 
is  measured  from  the  top  of  the  valence  band.  The  dashed 
curves  a.  h.  and  c  represent  exciton  levels  calculated  with 
£'j'*  =  0.1, 0.15.  and  0.3  eV.  respectively. 

by  replacing  the  short-range  Coulomb  potential  with  a 
spherical  square  well  of  variable  depth  and  a  screened 
Coulomb  tail.  Strinaii's  results  can  be  used  to  e.stimaie 
A,  corresponding  to  the  calculated  f*.  A/  decrea.ses 
rapidly  with  £*,  then  saturates  for  larger  £/,. 

The  contribution  to  the  natural  linewidth  from  the 
alloy  broadening  is  calculated  by  a  consideration  of  the 
electron  part  of  the  exciton  wave  function.  The  '/», 
is  expanded  in  a  linear  combination  of  the  spart  of  the 
conduction-band  wave  functions 

=  (7) 

n.s 

We  found  that  alloy  scattering  is  only  moderate  and  s 
scattering  is  dominant;  thus,  the  alloy  broadening 
A  ^  (  £)  is  well  approximated  by 

A4(£)  =  .v(1 -,v)S;lm£,(£).  (8) 

where  S,  is  the  difference  between  E^‘  and 
Hence,  the  alloy-broadening  contribution  to  A  is  relat¬ 
ed  to  the  alloy  broadening  of  the  band  states. 
A,(k,E) 

A,  =  4r2)(';',(k)lA^(£)|.l;,(k)) 

k 

■  k  n 

=  J^>,{E)^^{E)ilE 

=  V  ( 1  -  .X  Ifi'TT  J'/7;(f)i/£,  (9) 

The  integral  in  liq.  (9)  is  evaluated  mmierically 
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FIG  2.  Variation  of  A  (solid  lines)  and  A/  (dashed  lines) 
with  .vfor  three  £»  values. 

The  calculated  A.  which  is  the  sum  of  A/  and  A  <.  is 
plotted  against  .v  in  Fig.  2  for  three  values  of  £j'.  In  all 
three  panels,  the  dashed  curve  represents  A/  and  the 
solid  line  represents  A.  It  is  seen  from  Fig.  1  that  the 
exciton  level  follows  the  .V  edge  of  the  conduction 
band.  Hence  the  binding  energy  £;,.  relative  to  the 
conduction  band  edge,  remains  almost  constant  (for  a 
given  E^)  until  the  minimum  switches  from  the  V 
edge  to  the  L  edge.  Because  of  the  change  in  the  slope 
of  £,.  Ef,  decreases  rapidly  when  L  becomes  the 
minimum.  Correspondingly.  A/  varies  slowly  until  the 
.V  to  L  crossover  and  then  increases  rapidly.  This 
feature  is  clearly  seen  in  Fig.  2. 

For  £(l’  =  0  15,  the  A/  and  A^  are  comparable  near 
V  =  0  50.  and  A,  dominates  for  all  small  x  and  large  x. 
These  two  competing  mechanisms  give  a  relative 
minimum  ne.ir  v=:0.15.  a  broader  maximum  near 
.X  =0.50.  and  a  smaller  minimum  for  pure  silicon.  .As 
E^  is  decreased,  the  relative  minimum  is  shifted  to 
larger  x.  e  g  .  the  minimum  shifts  to  .x  =0.20  for 
£j’^()  10  cV.  For  7.(1'  =  0  15  eV.  the  pi'siiion  of  the 
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The  dislocation  energies  and  hardness  of  semiconductors  are  calculated  by  an  extension  of 
Harrison's  method.  It  is  demonstrated  in  agreement  with  experiment  that  dislocation  energies  per 
unit  length  are  proportional  to  d  “  d  ~  where  d  is  the  bond  length  and  hardness  is  proportional 
tod  “’-d  The  hardness  is  related  to  the  interaction  energies  among  dislocations.  It  is  argued 
that  dislocation  densities  of  semiconductors  will  be  reduced  if  they  are  alloyed  with  a  second 
constituent  that  has  a  shorter  bond  length.  Experimental  evidence  supporting  this  strategy  is 
noted. 
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Dislocations  in  semiconductors  are  detrimental  to  de¬ 
vice  function;  they  serve  as  channels  for  impurity  migration 
and  trapping,  which  cause  nonuniform  doping  and  degrades 
p-n  junctions.’  They  also  decrease  the  material’s  resistance 
to  plastic  deformation.  The  aim  of  this  letter  is  to  provide 
insights  into  the  underlying  physical  mechanisms  control¬ 
ling  dislocations  and  semiconductor  hardness,  and  then  to 
suggest  strategies  for  decreasing  dislocation  densities.  It  is 
well  established  that  the  hardness  of  tetrahedrally  coordi¬ 
nated  semiconductor  matenals — groups  IV,  III-V,  and  II- 
VI  compounds — exhibits  a  sharp  variation  with  their  near¬ 
neighbor  distance  d.  approximately  proportional  to  d  for 
one  group  of  seven  compounds.*  Thus,  semiconducton  with 
small  lattice  constants  tend  to  be  harder  materials.  These 
same  materials  have  larger  stiffness  coefficients'*  and  have 
fewer  dislocations  in  as-grown  crystals.^ 

The  shear  coefficients  (combinations  of  C44  and 
C||  —  C,;  in  the  Schoenflies  notation!’  depend  on  crystal 
orientation  and  I  in  Harrison's  notation”)  are  proportional  to 
V\/d''\V\  1^3)^’’,  where  1^2®^  is  the  covalent  and  Kj 

is  the  ionic  energy.  The  metallic  interaction  modifies  the 
functional  dependence  of  the  shear  coefficient  on  Kj  and  Kj, 
but  introduces  no  explicit  dependence  on  the  hopping  inte¬ 
grals.  denoted  K,  by  Harrison.*'"  In  a  pure  covalent  material, 
the  bond  energy  is  proportional  to  Vj  (or  d  "  ’),  and  the  bond 
volume  is  oc  d  *;  hence,  in  this  case,  the  shear  coefficient  var¬ 
ies  asd  In  the  limit,  Fj>  Kj,  C, ,  —  C,2<r  d  ~  For  most 
polar  semiconductors,  d  is  a  good  approximation. 

Hardness  is  determined  by  applying  a  known  force  F  to 
a  probe  of  a  prescribed  shape  driving  it  into  the  surface  of  the 
sample.’  The  area  A  of  the  resulting  indentation  is  measured, 
and  the  hardness  is  the  force  per  unit  indented  area.  Many 
dislocations  must  be  formed  to  allow  the  probe  to  indent  the 
semiconductor.  If  the  indenter  is  a  rectangular  pyramid, 
then  the  he  dnessis/f  =  F/A  =  Fh  /Ah  =  '*'here<'r 

is  the  work  required  to  cause  the  indenter  to  penetrate  to  a 
depth  h.  A  side  view  of  the  indentation  in  a  cut  through  its 
center  is  illustrated  schematically  in  Fig.  1.  The  top  of  the 
indentation  has  side  length  fV;  thus,  A  =  The  Burger’s 
vector  has  magnitude  b,  proportional  to  the  bond  length  d. 
The  number  of  dislocations  M/,  required  to  accommodate  an 


indentation  to  depth  A  is  iV*  =  A  /A  =  J  If'  cos  i?  /b,  where  i? 
is  the  angle  between  the  normal  to  the  tip  of  the  indenter  and 
a  side.  Figure  I  also  shows  a  model  of  one  possible  configura¬ 
tion  of  the  dislocations.  The  edges  of  the  extra  atom  planes 
that  are  driven  from  the  indented  volume  into  the  bulk  of  the 
semiconductor  are  shown  as  lines  terminated  by  dots.  The 
dotted  ends  of  these  lines  are  the  positions  of  the  disloca¬ 
tions,  which  are  perpendicular  to  the  plane  of  the  figure.  The 
planes  driven  to  the  sides  each  have  a  finite  extent  and  a 
trapezoidal  shape.  The  planes  driven  down  under  the  indent¬ 
er  have  a  square  shape. 

Much  of  the  work  done  on  the  indenter  goes  into  the 
energy  to  form  the  indicated  dislocation  configuration,  al¬ 
though  some  certainly  goes  into  heat.  There  are  two  major 
contributions  to  this  formation  energy.  The  first  is  the  enegy 
needed  to  generate  each  dislocation  as  an  isolated  entity,  and 
the  second  is  the  interaction  energy  among  these  disloca¬ 
tions.  Because  the  interaction  term  dominates  If,  approxi¬ 
mations  made  to  simplify  the  first  term  are  relatively  unim¬ 
portant.  The  extra  planes  driven  to  the  sides  of  the 
indentation  have  a  finite  extent;  accordingly  (in  this  idealized 
picture),  there  are  both  edge  dislocations  at  their  base  and 
screw  dislocations  associated  with  their  termination.  The 
square  planes  driven  below  the  indentation  have  edge  dislo¬ 
cations  around  the  sides  and  screw  dislocations  at  the 
comers  to  make  the  turns.  Moreover,  there  are  interactions 


FIG.  1.  Schematic  representation  of  an  idealized  minimum-energy  disloca¬ 
tion  configuration  produced  by  a  square  cross-section  indenter  in  a  hard¬ 
ness  measurement.  If'is  the  side  length  and  A  is  the  depth  of  the  indentation 
In  this  ideal  ca.se.  half  the  matenal  from  the  indented  region  is  displaced 
along  the  glide  planes  (indicated  by  the  dashed  arrotvsl  to  the  sides  and  half 
is  displaced  below  the  indenter. 
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among  the  dislocations,  which  can  produce  a  minimum-en¬ 
ergy  configuration.  For  the  arrangement  depicted  in  Fig.  1, 
it  always  costs  energy  to  position  a  second  dislocation  on  a 
parallel  glide  plane  to  one  already  present.  However,  the 
magnitude  of  this  extra  energy  can  be  minimized  and,  for 
proper  configurations  of  the  dislocations,  there  are  attractive 
forces  along  the  glide  planes  that  will  tend  to  position  the 
dislocations  into  the  minimum-energy  configuration.  The 
minimum  configuration  arises  when  half  of  the  atoms  from 
the  indented  volume  go  respectively  to  the  side  and  below  the 
indenter.  Then,  in  both  regions,  the  maximum  angle  made 
between  successive  close-spaced  dislocation  lines  and  their 
glide  planes  is  i?  =  w/2  —  t?  /2,  as  shown  in  the  figure.  This  is 
the  minimum  realistic  energy  configuration.  If  the  disloca¬ 
tions  are  separated  more  than  shown  in  Fig.  1,  then  there  is 
more  volume  of  strained  material  and  the  interaction  energy 
would  be  larger  still. 

An  approximate  expression  for  the  energy  required  to 
indent  the  material  is’ 


er  =  24f  f  E,L,  -h  1  f  •£,  [min(L,^)l  (1) 

Li  „  I  2  „  I 

In  the  first  term,  using  an  isotropic  medium  approximation 
and  neglecting  core  terms,  the  energy  per  unit  length  to  form 
an  edge  dislocation  is’ 


Gb-  ,  R 
- In — , 

4)ri  1  -  V)  to 


(2) 


the  shear  coefficient  is  G,  the  Burger's  vector  b  =  d  /3  for  an 
indentation  along  a  ( 1 OOV  axis,  the  range  of  the  elastic  defor¬ 
mation  of  a  dislocation  R  is  taken  equal  to  tV  (for  want  of  a 
better  approximation!,  r^  —  d  is  the  dislocation  core  radius, 
the  Poisson  ratio  is  v~0.2  for  most  semiconductors,  and  L, 
is  the  length  of  the  I'th  dislocation.  In  the  second  term,  is 
the  interaction  energy  per  unit  length  between  dislocations  i 
and  j.  Assuming  they  have  parallel  glide  planes  and  their 
Burger’s  vectors  have  the  same  sign,  £,,  is  given  by’ 


where  is  the  separation,  and  <p,^  is  the  angle  that  a  line 
perpendicular  to  and  joining  the  dislocations  makes  with  the 
glide  plane,  as  shown  in  Fig.  1.  Because  the  various  disloca¬ 
tions  in  a  region  have  different  lengths,  the  net  interaction 
energy  is  approximated  by  multiplying  the  energy  per  unit 
length  by  the  length  of  the  shorter  one.  The  upper  limits  on 
the  sums  .Vare  the  number  of  dislocations  in  one  region  (side 
or  bottom  I  associated  with  one  edge.  For  the  minimum-ener¬ 
gy  configuration,  .V  =  .V,/2.  The  four  that  multiplies  the 
bracket  accounts  for  the  four  sides,  and  the  two  for  the  two 
regions  for  each  side. 

We  now  encounter  our  first  surprise.  As  we  can  see 
from  Eq.  I2land  (3),  and  £  have  comparable  magnitudes. 
Because  there  are  approximately  .V  ^  terms  in  the  interaction 
energy  sum,  only  .V  terms  in  the  formation  energy  sum,  and 
iVis  a  large  number  li\>  1),  the  interaction  energy  completely 
dominates  the  hardness.  In  fact.  A’  is  typically  of  the  order 
10’.  Hence,  terms  owing  to  screw  dislocations,  core  energies, 
heat  dissipation  as  the  dislocations  propagate  to  their  places, 
and  other  effects  associated  with  the  first  term  are  unimpor¬ 


tant.  However,  care  must  be  taken  with  the  interaction 
terms;  Eq.  11)  neglects  a  number  of  secondary  interactions, 
some  positive  and  others  negative;  these  will  be  added  later. 
The  principal  neglected  terms  are  the  interactions  between 
the  dislocations  in  the  different  regions  on  each  side  (posi¬ 
tive!,  and  the  interactions  between  adjacent  side  and  bottom 
regions  (positive),  and  the  interaction  between  the  opposite 
sides  in  the  bottom  regions  (negative).  Comparison  of  the 
results  with  expenment  will  indicate  how  impiortant  these 
neglected  terms  are  likely  to  be.  TTie  length  £,  is  IF  (A  —  ij/A' 
for  /  from  1  to  .V’.  This  is  the  largest  length  of  the  side-inserted 
planes  and  its  choice  partially  accounts  for  interacitons 
between  the  otherwise  neglected  screw  dislocations.  The  dis¬ 
tance  =  r,  —  is  given  by  l^llb  i  —  j\  for  /  and  j  ranging 
from  1  to  S  for  the  minimum-energy  configuration  and  a 
tetrahedrally  bonded  semiconductor.  Finally,  in  the  indicat¬ 
ed  configuration,  =if  =  t/2  —  t?  /2.  Inserting  these  ex¬ 
pressions  into  Eq.  ( 1 )  and  retaining  only  terms  of  order  ,V  ’ 
yields 


H 


min 


C  cot  t?  cot  t?  \ 

6m  1  -  V)  V  2 


4 

-h  -I-  sin 


2 


(4) 


One  can  also  get  a  number  for  the  hardness  of  a  disloca¬ 
tion  in  which  all  the  material  is  pushed  along  the  same  glide 
plane,  e.g..  to  the  side,  to  the  bottom,  or  normal  to  the  face  of 
the  indenter  la  possibility  not  depicted  in  Fig.  1 ).  In  this  case, 
the  factor  of  2  in  front  of  Eq.  Ill  is  removed.  .V  =  .\\  and 
d)  =  ir/2  —  t?.  Then  a  higher  nonequilibrium  hardness  in  the 
context  of  this  model  (denoted  H,),  is  obtained. 


G  cot  I? 


The  proper  answer  for  most  materials,  and  depending  on 
crystal  onentation.  probably  lies  somewhere  between 
and  H,.  For  an  indenter  with  i?  =  n-/4,  we  have 
=  0.0969(7/(1  —V)  and  =2.39.  Harrison*  has 

shown  that  one  coninbution  to  the  shear  coefficient  (actually 
C, ,  —  Cjj)  is  G  =  2.38  frmal/md^,  where  m  is  the  free- 
electron  mass,  is  the  covalence,  a^.  =  FVIFj  +  Fj  )'■'-, 
and  d  is  the  bond  length.  We  will  approximate  G  by  this 
expression.  Using  this  G  and  v  =  0.2,  and  changing  the  di¬ 
mensions  to  those  in  terms  of  which  experimental  hardness 
numbers  are  customarily  quoted  gives  H„,„  =  2.38  X  lO* 
\al/d ')  kgm/mm^  where  if  is  in  angstroms.  Calculated  val¬ 
ues  of  and  are  plotted  against  experimental  results 
in  Fig.  2  for  a  number  of  semiconductors. 

Figure  2  has  the  theoretical  and  //,  values  con¬ 
nected  by  arrows  from  to  //,  for  each  compound,  plot¬ 
ted  as  a  function  of  the  corresponding  experimental  val¬ 
ues.''*  If  the  theory  were  perfect  and  the  exp  .nmental  values 
w'ere  accurate,  the  points  would  fall  on  the  indicated  unity 
slope  line.  Several  conclusions  can  be  drawn.  Firstly,  the 
order  of  magnitude  of  the  predicted  and  measured  values  are 
the  same,  a  result  obtained  with  no  adjustable  parameters  in 
the  theory.  Secondly,  the  trends  from  one  compound  to  an¬ 
other  are  properly  given  by  the  theory.  Although  the 
values  are  generally  too  small,  they  fit  the  soft  materials  bet¬ 
ter,  and  the  //,  values  fit  the  harder  materials  better.  Thirdly, 
from  Eqs.  (1)  and  (3),  H  is  given  in  a  rough  but  revealing 
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FIG.  2.  Theoreticai  vs  oxpMnmentai  hardness  of  several  semiconductors  in 
Refs.  2  and  9.  The  two  theoretical  values  for  each  semiconductor  are  con* 
necied  by  an  arrow  from  //„,„  to  //,.  Perfect  agreement  would  correspond 
tc  points  being  on  the  unity-slope  solid  line. 

approximation  by  4.V’£^Z  /Ah,  where  the  averge  disloca¬ 
tion  length  Z  =  W'/2  and  £,,  is  the  average  dislocation  pair 
interaction  energy  per  unit  length.  Notice  that  .V  (or  IV]  can¬ 
cels  from  this  expression;  thus,  //is  independent  of  W  (orF) 
and  therefore.  H  is  truly  a  measure  of  the  properties  of  the 
material.  This  result  would  not  be  found  if  the  dislocation 
energies  ( « .V )  were  to  dominate  H  rather  than  the  pair  inter¬ 
action  energy  (<i.V‘).  Finally,  the  Berger’s  vector  cancels 
from  the  leading  term  and  appears  only  in  the  argument  of 
the  logarithm  in  Eq.  |3|.  Thus,  the  answers  are  also  insensi¬ 
tive  to  Its  choice. 

Dislocations  are  often  found  in  materials  as  they  are 
grown.  Their  density  is  determined  by  the  thermal  and  me¬ 
chanical  stresses  to  which  they  are  subjected  in  the  growth 
process.  A  dislocation  constitutes  a  metastable  excitation 
relative  to  the  perfect  crystal  ground  state.  At  the  elevated 
growth  temperatures  and  temperature  gradient  behind  the 
growth  front,  the  number  of  dislocations  present  is  con¬ 
trolled  by  the  relative  rate  at  which  vacancies  anneal  or  con¬ 
dense  into  dislocations. The  dislocation  formation  rate  will 
be  slower  in  a  material  grown  at  the  same  temperature  if  £,  is 
higher.  If  an  alloy  is  formed  from  a  material  of  interest  and  a 
second  constituent  with  a  shorter  bond  length,  one  expects 
the  average  bond  energy  (and  thus  both  the  melt  temperature 
and  vacancy  formation  energy!  to  increase  proportional  to  a 
low  inverse  power  of  the  average  bond  length.*  Hence,  the 
equilibrium  vacancy  density  just  below  the  melting  point 
tends  to  be  the  same  in  lowest  order  for  all  materials,  inde¬ 
pendent  of  the  bond  lengths  of  the  constituents.  However, 
the  shear  coefficient  and  dislocation  energy  per  unit  length 
will  increase  with  much  higher  inverse  powers  of  the  bond 
length.  Consequently,  dislocation  densities  should  be  re¬ 
duced  in  such  alloys  relative  to  those  found  in  the  longer 
bond  length  pure  constituent.  This  expectation  is  conhrmed 


in  the  recent  work  reported  on  Zn ,  _  ^  Cd,  Te  bulk  material.* 
The  best  CdTe  that  has  ever  been  grown  has  dislocation  den¬ 
sities  in  excess  of  5  X  1(F  cm  ”  ^  The  addition  of  only  4%  Zn 
reduced  the  dislocation  count  to  less  than  5  X  10*  cm  ”  The 
ZnTe  bond  length  is  2.643  A,  while  the  CdTe  bond  length  is 
2.805  A,  a  6%  difference.  This  6%  difference  in  bond  length 
translates  into  a  2%  difference  in  the  dislocation  energy  per 
unit  length  for  1  —  jc  =  0.04.  Dislocation  energies  per  unit 
length  are  typically  10  eV  per  lattice  spacing;  accordingly,  a 
2%  increase  can  be  expected  to  slow  their  formation  rate 
considerably. 

The  argument  just  presented  naturally  leads  to  a  strate¬ 
gy  for  decreasing  dislocations  in  other  semiconductors.  If  an 
alloy  is  made  of  the  material  of  interest  with  another  com¬ 
pound  with  a  shorter  bond  length,  then  the  dislocation  den¬ 
sity  should  be  reduced.  For  example,  this  suggests  that  the 
addition  of  a  small  amount  of  GaP  [d  =  2.359  A)  may  signifi¬ 
cantly  reduce  the  dislocation  density  of  bulk  grown  GaAs 
{d  =  2.448  A).  It  has  been  demonstrated  that  the  addition  of 
approximately  1%  GaN  id  =  1.946  A)*  *  or  of  a  10‘*  cm~’ 
BAs  {d  =  2.069  A)  concentration ' '  to  GaAs  can  yield  a  large 
volume  of  dislocation-free  material.  An  InAs  additive  with 
its  longer  bond  length  [d  =  2.623  A)  serves  the  same  function 
indirectly,  by  causing  GaAs  bonds  in  its  neighborhood  to  be 
compressed.  This  indirect  mechanism  should  be  less  effec¬ 
tive  than  substituting  short  bond  length  additives. 

We  have  demonstrated  that  the  dislocation  energies 
and  hardness  of  tetrahedrally  bonded  semiconductors  are 
rapid  functions  of  the  reciprocal  of  the  bond  length.  This 
rapid  d  dependence  of  dislocation  energies  provides  a  ratio¬ 
nale  for  the  dramatic  decrease  of  the  dislocation  density  in 
bulk  grown  Zn^  q,  Cd,,  Te  material  relative  to  that  found  in 
CdTe,  and  suggest  means  for  accomplishing  the  same  ends  in 
other  materials. 
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The  sensitivity  of  defect  energy  levels  in  semiconductors  to  the  host  band  structures  and  impurity 
potentials  has  been  studied  for  approximately  30  impunties  in  CdTe  using  four  different  band- 
structure  models.  The  discrepancies  in  the  defect  levels  between  two  different  sets  of  band  struc¬ 
tures  and  impurity  potentials  are  found  to  range  from  less  than  O.I  eV  to  the  whole  band  gap  (1.6 
eVi.  The  band-structure  effects  are  analyzed  here  in  terms  of  detailed  partial  densities  of  states. 

Examples  of  contradictory  predictions  from  different  band  structures  are  illustrated,  and  ways  to 
improve  the  theory  are  suggested. 


I.  INTRODUCTION 

In  several  of  our  recent  papers,'”^  we  have  applied  a 
method  to  calculate  the  band  structure  of  semiconductors 
that  is  both  efficient  and  accurate.  Because  the  procedure 
involves  casting  the  basis  functions  into  orthonormal  local 
orbitals"  (OLO),  our  method  has  the  advantages  common 
to  empirical  tight-binding  (ETB)  calculations,*" except 
that  the  Hamiltonian  matrix  elements  to  all  ranges  are  re¬ 
tained.  The  inclusion  of  these  higher  coefficients  makes  it 
possible  to  produce  excellent  band  structures  including 
conduction  bands  and  effective  masses.  The  method  also 
yields  wave  functions  for  optical  property  calculations.’ 
Moreover,  its  OLO  description  also  permits  its  extension, 
through  the  coherent-potential  approximation,  to  al¬ 
loys.'”* 

The  recent  attention  focused  on  defects  in  semiconduc¬ 
tors  motivated  us  to  apply  our  method  to  this  problem. 
The  theories  of  defects  have  ranged  from  very  sophisticat¬ 
ed  self-consistent  density-functional  theory*'”'*  (SCDF) 
to  simple  ETB  calculations.  It  is  generally  recognized 
that  SCDF  is  as  accurate  in  defects  for  the  ground-state 
properties  as  it  is  for  pure  semiconductors,  but  less  certain 
in  assigning  excited  energy  levels.  ETB,  because  it  can 
produce  results  for  many  systems  in  one  study,  claims  to 
predict  the  trends  of  deep  levels'"  even  if  the  accuracy  for 
a  given  impurity  may  be  poor.  However,  this  contention 
remains  to  be  verified. 

To  assess  this  concern,  we  ask  the  following  question: 
“How  sensitive  are  defect  levels  to  host  band  structures 
and  impurity  potentials?"  To  this  end.  we  have  adopted 
the  simple  yet  nontrivial  defect  model,  t’r.it  of  site- 
diagonal  substitutional  defects  often  used  in  El  B  studies. 
CdTe  was  selected  in  this  study  because  its  band  structure 
has  been  examined  in  great  detail  by  us.  and  there  are 
three  published  band-structure  models*”'"  '..  at  we  could 
easily  generate  for  comparison.  There  is  also  a  consider¬ 
able  body  of  experimental  data  on  deep  states  in  this  sys- 


II.  CALCULATIONAL  PROCEDURE 

In  the  simple  site-diagonal  substitutional  defect  model, 
the  impurity  energy  levels  E  are  determined  by  the  equa¬ 
tion 

l-u„g„(£)=0,  (1) 

where  a  designates  the  symmetry  of  a  local  state,  e.g.,  P 

r?,  and  Pg  on  an  atomic  site  in  the  zinc-blende  structure, 
and  ga  is  the  real  part  of  the  diagonal  matrix  element  of 
the  host-crystal  Green  function,  can  be  calculated 
from  the  partial  density  of  states  (PDOS)  by 

gJE)=Jpa(€)AE-€)de.  (2) 

The  PDOS  is  given  by 

,  (3) 

n.k 

where  e„(k)  are  band  energies  and  a“(k)  are  the  probabili¬ 
ty  amplitudes  of  the  band  state  in  the  Bloch  basis  con¬ 
struct^  from  the  OLO  labeled  by  a.  The  Brillouin-zone 
integration  in  Eq.  (3)  is  calculated  using  an  accurate  ray 
scheme.'* 

Because  a  principal  concern  of  this  paper  is  the  sensi¬ 
tivity  of  impurity  levels  to  the  host  band  structures,  we 
should  emphasize  the  difference  between  our  method  and 
ETB.  Our  method  consists  of  four  steps. 

(1)  We  start  with  four  Gaussian  orbitals  per  atom  and 
empirical  pseudopotentials,'’  and  compute  the  Hamiltoni¬ 
an  matrix  f/(k)  and  overlap  matrix  S(k)  as  was  done  by 
Kane*"  and  Cliadi.*' 

(2)  The  Gaussian  orbitals  are  transformed  into  OLO," 
so  //{k)  is  transformed  into  H^^ik)  and  S  into  the  identity 
matrix.  The  band  structures  calculated  from  /folk*  are 
accurate  to  5Cc  as  compared  to  more  sophisticated 
methods  using  the  same  potential.' 

(3)  A  spin-orbit  Hamiltonian  in  the  OLO  basis'*  is  in¬ 
corporated  to  deal  with  this  interaction. 

(4)  To  compensate  for  the  effects  of  truncated  basis  and 
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nonlocal  potentials,  a  penurbation  Hamiltonian  Hx  is 
added.  //|  has  the  same  form  as  a  truncated  ETB  Hamil¬ 
tonian.  The  parameters  in  Hx  are  adjusted  to  fine  tune 
the  important  band  energies  and  effective  masses. 

Although  both  ETB  and  our  methods  are  empirical, 
there  are  two  major  differences. 

(1)  While  most  ETB  retains  the  H  matrix  elements  only 
to  the  first-  or  second-neighbor  shell,  ours  extends  to  all 
ranges,  so  that  the  high  Fourier  components  needed  to 
produce  the  sharp  band  curvatures  are  properly  given. 

(2)  Our  method  can  directly  generate  wave  functions  for 
calculation  of  other  propenies. 

Thus,  while  our  method  yields  more  accurate  band  struc¬ 
tures,  it  retains  much  of  the  advantage  of  ETB.  namely 
the  computational  speed  and  a  simple  direct-space 
description  of  the  Hamiltonian. 


III.  BAND  STRUCTURES  AND  PARTIAL 
DENSITIES  OF  STATES 

Figure  1  depicts  the  four  band  structures  to  be  con¬ 
sidered  for  CdTe.  Our  result  is  in  panel  (a);  panels  (b) 
(Ref.  8)  and  (c)  (Ref.  9)  are  two  ETB  band  structures  with 
the  Hamiltonian  matrix  elements  truncated  at  second 
neighbors.  (Because  different  parameters  were  selected, 
these  two  band  structures  are  not  identical.)  Panel  (d) 
(Ref.  10)  results  from  the  use  of  five  basis  orbitals  per 
atom;  the  extra  one  is  an  excited  s  state.  All  these  band 
structures  are  adjusted  to  have  the  proper  fundamental 
band  gap  of  1.6  eV.  The  principal  differences  one  sees  on 
first  inspection  are  in  the  band  curvatures,  especially  the 
conduction  bands.  The  effective  mass  at  the  bottom  of 
the  conduction  band  in  panel  (a)  is  0.1  times  the  free- 
electron  mass,  in  agreement  with  experiment,'^  while  in 
other  panels  it  is  more  than  twice  as  large. 

Figure  2  shows  the  densities  of  states  (DOS)  for  each  of 
the  band  structures  in  Fig.  1.  WTiile  the  valence  bands  at 
least  exhibit  general  common  features,  the  conduction 
bands  are  almost  unrecognizable  as  representing  the  same 
compound.  In  panels  (c)  and  (d),  for  example,  there  is  a 
second  band  gap  above  the  fundamental  gap.  Also  note 
that  there  are  two  extra  narrow  peaks  associated  with  the 
two  extra  excited  s  orbitals  (one  for  Cd  and  the  other  for 
Te)  included  in  the  calculation. 

To  analyze  the  band  effects  on  defect  levels  [see  Eqs.  (1) 
and  (2  ],  the  DOS  is  further  decompo.sed  into  partial  den¬ 
sities  of  states  for  Fols),  F-lp'^-),  and  r,(p^^-)  states  on 
the  Cd  and  Te  sites,  as  shown  in  Figs.  3—6.  The  Tg 
PDOS  are  not  shown  because  they  are  nearly  the  same  as 
r7  with  only  a  slight  upward  energy  shift.  These  PDOS 
show  how  the  “atomic"  levels  evolve  into  band  states. 
These  curves  contain  useful  information  about  many 
properties,  e.g.,  the  relation  between  the  crystal  bonding 
and  atomic  energies,  and  how  potential  disorder  in  alloys 
affects  different  parts  of  the  bands, in  addition  to  de¬ 
fect  levels  studied  here. 

The  FetCd)  PDOS  shown  in  Fig.  3  split  between  the 
conduction  and  valence  bands.  It  is  generally  assumed 
that  the  cation  s  states  in  III-V  and  II-VI  compounds 
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FIG.  1.  Four  band  structures  of  CdTe  used  for  comparative 
studies:  (a)  present  work,  ibi  Ref.  8.  (c)  Ref.  and  d)  Ref.  10. 


FIG.  2.  Densities  of  states  calculated  fo'm  the  four  band 
structures  in  Fig.  1. 
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evolve  into  the  conduction  bands,  while  the  anion  p  states 
make  up  most  of  the  major  valence  bands  just  below  the 
gap.  Thus  it  is  perhaps  a  surprise  to  see  a  prominent  peak 
denved  from  the  cation  s  states  at  the  bottom  of  the  ma¬ 
jor  valence-band  structure.  However,  this  is  a  general 
feature  for  all  sp^-based  compound  semiconductors. 
These  are  the  states  responsible  for  the  first  observed 
breakdown  of  the  virtual-crystal  approximation  for  a 
semiconductor  alloy:  Hgi.j^Cd^Te  (which  is  caused  by 
the  large  s-energy  shift  between  the  Cd  and  Hg  sites).'*’^'^^ 

A  more  detailed  examination  draws  attention  to  some 
important  differences  among  the  four  panels  in  Fig.  3: 
the  valence-band  peak  in  panel  (c)  is  about  2  eV  higher 
than  the  rest,  and  it  is  also  high  compared  to  experi¬ 
ment."  Our  conduction-band  PDOS  in  panel  (a)  is 
broader  than  the  others.  The  ratio  of  the  integrated 
PDOS  in  the  conduction  bands  to  that  in  the  valence 
bands  in  our  model  is  larger  than  those  in  other  panels. 
Also  our  PDOS  just  below  the  valence-band  edge  is  obvi¬ 
ously  smaller  than  that  found  in  other  models. 

Figure  4  shows  that  the  Cd  p  states  are  concentrated  in 
the  conduction-band  states.  This  is  particularly  true  in 
panel  (a),  where  their  contribution  to  the  valence-band 
states  shrinks  almost  to  nothing.  In  other  panels,  there 
are  still  sizable  {~2Q%)  valence-band  states.  In  contrast. 


all  four  panels  in  Fig.  5  show  that  the  Te  s  states  are  con¬ 
fined  to  the  deep  valence-band  states,  as  generally  recog¬ 
nized.  Finally,  Fig.  6  shows  that  the  Te  p  states  dominate 
the  upper  valence-band  states.  Panel  (a)  has  much  less 
conduction-band  content  than  the  other  three  panels.  As 
we  will  see,  these  differences  can  result  in  quantitatively 
or  even  qualitatively  different  predictions  about  the  deep 
levels. 

IV.  IMPURITY-LEVEL  DETERMINATION 

A  convenient  way  to  study  the  impurity  energy  levels 
using  Eq.  (1)  is  to  rewrite  it  as  !;„  =  l/gaiE)  and  plot  E  as 
a  function  of  v.  Once  this  E-u  curve  is  deduced  for  each 
a,  the  deep  levels  for  a  given  impunty  can  be  read  off 
the  curve  by  drawing  a  vertical  line  at  the  appropriate 
value  of  Ua  for  the  impurity.  We  set  the  zero  of  energy  at 
the  top  of  the  valence  bands.  Because  the  gap  is  1.6  eV, 
we  will  focus  on  levels  in  the  energy  range  from  —0.5  to 
2.0  eV. 

Calculations  have  been  performed  for  all  neutral  impur¬ 
ities  listed  in  Table  I.  Because  we  do  not  believe  that 
there  exists  a  uniformly  accepted  table  for  v  we  have 
adopted  a  table  that  we  used  for  structural  studies.*^'*'* 
Table  I  lists  the  term  values,  which  we  obtained  from  to- 


TABLE  I.  r-  and  p-state  correlated  term  values  in  units  of  — eV.  The  top  entry  is  the  r-state,  the 
second  the  pi/2-state,  and  the  third  the  pi/j-state  energy.  (All  energies  are  negatice.) 
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FIG.  7.  The  E-v  curves  for  the  P 4  states  on  a  Cd  site. 

tal  energy  differences  between  atomic  configurations  cal¬ 
culated  using  the  norm-conserved  pseudopotentials’’  and 
self-consistent  charge-density-functional  theory,  with  the 
first  ionization  energies  adjusted  to  be  the  experimental 
values.’”  These  term  values  are  found  to  yield  consistent¬ 
ly  better  structural  properties^  in  Harrison’s  theory’’  ** 
than  those  based  on  Mann’s  values^’  adopted  by  Har- 
nson.’*  The  impurity-potential  parameters  will  then  be 
taken  as  the  difference  of  the  term  values  between  the  im- 
punty  atom  and  Cd  (or  Te).  To  study  the  sensitivity  of 
£3  to  L'a,  we  shift  i)<,  by  ±0.5  eV  and  compute  the  corre¬ 
sponding  changes  in  the  energy  levels. 

Figures  7—10  display  the  E-v  curves  for  several  a. 
Each  figure  has  four  curves,  corresponding  to  the  four 
panels  of  PDOS  in  each  of  Figs.  3—6.  The  functional 
behavior  of  these  curves  can  be  understood  qualitatively 
using  Eq.  (2)  and  Figs.  3—6.  If  £  lies  in  the  gap,  the  con¬ 
tribution  from  conduction  bands  is  negative,  but  positive 
from  the  valence  bands.  The  closer  the  PDOS  to  the  £  in 
question,  the  larger  will  be  its  influence.  Applying  this 
argument  to  the  r<|(Cd)  representation,  we  see  that  the 
curves  in  Fig.  7  are  negative  in  the  gap  region  because  the 
PDOS  in  Fig.  3  near  the  bottom  of  the  conduction  bands 
are  much  larger  than  those  near  the  valence-band  top. 
Thus,  on  the  Cd  site,  only  impurities  with  an  s  energy 
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FIG.  9.  The  E-v  curves  for  the  F t,  states  on  a  Te  site. 


below  the  Cd  s  level  ( —  8.99  eV)  will  produce  a  level  in 
the  gap.  However,  we  note  that  in  Fig.  7,  ga(£)=0  for 
models  (b)  and  (c)  just  below  the  valence-band  edge  be¬ 
cause  of  cancellation  between  the  conduction-  and 
valence-band  contributions.  At  this  £  value,  the  E-v 
curve  switches  from  v  =  —  oo  to  u  =  oo  (not  shown);  an 
ideal  vacancy  level  (corresponding  to  Uq  =  00 )  is  located  at 
this  £.  A  similar  consideration,  but  with  the  conduction 
and  valence  bands  interchanged,  leads  to  an  understanding 
of  the  curves  in  Fig.  10.  Using  the  same  principle,  we  can 
easily  understand  why  all  curves  in  Fig.  9  for  the  PilTe) 
representation  are  positive,  but  the  reasons  for  the  large 
displacements  between  these  curves  are  not  easy  to 
deduce.  In  Fig.  8,  the  curve  labeled  a  is  distinctively  dif¬ 
ferent  from  other  curves,  because  the  PDOS  in  panel  (a)  in 
Fig.  4  is  completely  dominated  by  the  conduction  band; 
however,  for  the  other  panels  the  PDOS  just  below  the 
valence-band  edge  are  as  large  as  those  just  above  the 
conduction-band  edge.  This  produces  a  very  sharp  nega¬ 
tive  E-u  curve  for  (a),  but  split  behavior  for  (b),  (c),  and 
(d). 

These  E-v  curves  provide  a  clear  picture  of  how  dif¬ 
ferent  host  band  structures  may  affect  the  deep  levels. 
Numerical  values  for  the  impurity  levels  can  be  obtained 
from  these  figures  by  drawing  vertical  lines  at  the  ap- 
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propriate  impurity  potentials  (i.e.,  differences  between  the 
term  values  listed  in  Table  I),  as  has  been  shown  for 
several  representative  impurities.  To  provide  a  more 
quantitative  comparison.  Table  II  lists  some  calculated 
impurity  levels  £„  and  the  corresponding  changes 
due  to  the  1-eV  change  in  Vg. 

V.  RESULTS  AND  CONCLUSION 

To  summarize  we  recall  that  band  models  (b)  and  (c) 
are  the  same  second-neighbor  ETB  with  two  different  sets 
of  parameters,  and  model  (d)  is  a  first-neighbor  ETB  with 
one  e.\tra  s  orbital  per  atom.  Our  model  [model  (a)]  has 
the  form  of  ETB  but  is  derived  in  a  very  different  manner 
and  includes  all  the  long-range  interactions.  Therefore, 
we  e.spect  that  the  results  from  models  (b)  and  (c)  will  be 
close,  model  (d)  will  have  larger  discrepancies  from  (b) 
and  (c)  than  that  between  (b)  and  (c),  and  model  (a)  will 
differ  even  more.  This  is  evident  from  Figs.  7—10  and 
Table  II.  We  found  the  energies  for  the  r^lCd),  r7(Te), 
and  FijtTe)  states  produced  by  models  (b)  and  (c)  agree 
within  0.1  eV.  For  the  other  states,  i.e.,  FolTe),  r^iCd), 
and  r^lCd),  the  energies  from  (b)  and  (c)  are  qualitatively 
similar,  but  the  difference  can  be  as  large  as  0.4  eV.  The 
largest  discrepancy  between  models  (d)  and  (b)  [or  (c)]  is 
more  than  0.5  eV,  and  that  between  (a)  and  other  models 
is  more  than  1  eV.  The  largest  difference  comes  from  the 


p  levels  on  a  Cd  site.  For  example,  the  filled  p  level  of  C 
on  a  Cd  site  in  model  (a)  is  a  resonance  state  Just  below 
the  valence-band  edge  but  is  a  donor  state  in  the  other 
models.  Similarly,  model  (a)  puts  the  neutral  Te  antisite 
defect  p  levels  at  about  \  and  j  of  the  gap  [£(r7)=0.48 
eV  and  £(rg)=0.95  eV],  while  other  models  assign  them 
as  resonance  states  inside  the  conduction  bands.  We  also 
note  that  the  discrepancies  between  different  models  are 
not  uniform,  but  vary  with  u„.  Consider  FJCd)  for  ex¬ 
ample.  All  four  models  yield  the  same  ordering  and  about 
the  same  energies  for  the  group-III  impurities  Al,  In,  and 
Ga.  However,  as  u  becomes  more  negative,  the  splitting 
between  the  curves  increases,  so  the  discrepancies  become 
larger  [  ~  1  eV  difference  between  models  (a)  and  (d)  for  I 
impurity].  Similarly,  for  the  F 7iTe)  states,  all  four  models 
put  the  Sn  impunty  energies  close  to  the  valence-band 
edge,  but  the  agreement  deteriorates  as  increases. 

Regarding  the  sensitivity  of  energy  levels  to  impurity 
potentials.  Table  II  shows  that  a  1-eV  shift  in  produces 
a  change  in  £„  ranging  from  less  than  0.1  to  0.65  eV. 
Very  little  is  known  about  the  size  or  trends  in  errors  in¬ 
troduced  in  y„  from  the  use  of  atomic  term  values.  How¬ 
ever,  we  know  that  the  discrepancy  of  u„  between  two  dif¬ 
ferent  tables  of  atomic  term  values  can  be  larger  than  2 
eV.  This  discrepancy  translates  into  an  uncertainty  of  less 
than  0.1  to  more  than  1  eV  in  the  impurity  energy  levels. 


T.A.BLE  11.  Defect  energy  levels  £  and  changes  A£  due  to  a  1-eV  change  in  the  impurity-potential  parameter.  All  energies  are  in 
units  of  eV.  F,,  stands  for  ideal  vacancy. 
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which  is  comparable  to  that  due  to  different  host  band 
structures. 

Putting  this  large  uncertainty  in  the  deep  levels  against 
a  band  gap  of  1.6  eV,  we  are  left  with  great  doubts  about 
the  predictability  of  this  oversimplified  theory.  Unfor¬ 
tunately.  the  expenmental  means  available  for  identifying 
microdefects  in  semiconductors  are  still  very  limited,  and 
the  ab  imtio  band  theory  is  still  not  capable  of  accurately 
predicting  the  energy  levels.  Thus,  there  is  a  great  temp¬ 
tation  to  use  simple  theories  like  the  one  carried  out  here 
to  help  with  the  identifications.  To  illustrate  this  point, 
consider  the  following  e.xamples;  Table  II  shows  that  Li 
on  3  Te  site  has  an  s  level  of  0.14  eV  in  model  (a),  so  one 
may  be  tempted  to  relate  it  to  the  acceptor  state  identified 
e.xpenmentally.'^  However,  this  is  not  the  hydrogenic  ac¬ 
ceptor  state  on  a  Cd  site,  as  one  might  anticipate.  One 
might  also  want  to  assign  the  and  t  gap  states  for  the 
Te  antisite  p  levels  on  the  Cd  site  found  from  model  (a)  as 
those  seen  in  e.xperiments.'^- Because  of  the  large  uncer¬ 
tainty  in  the  calculation,  these  results  should  be  regarded 
as  suspicious  surprises  rather  than  theoretical  confirma¬ 
tions. 

The  results  presented  here  should  not  discourage  con¬ 
tinued  research  on  the  ETB  approach,  but  improvement  is 
clearly  needed.  Work  ranging  from  universal''’  *’'^®  to 
specific''*  ' structural  studies  to  our  band  calculations 
and  alloy  studies’ indicates  that  the  ETB  type  of  theory 
is  practical  for  both  bonding  properties  and  electronic 
structures.  The  rea.son  that  ETB  works  well  for  some 
properties,  e.g..  photoemission  spectra  and  bonding  prop¬ 
erties,  but  not  for  impurity  levels,  is  that  the  former  de¬ 
pend  only  on  the  gross  total  density  of  states,  while  the 


latter  have  been  shown  to  be  sensitive  to  the  details  of  the 
panial  densities  of  states. 

To  establish  the  credibility  of  ETB  in  defect  studies, 
one  needs  to  look  at  the  problem  more  seriously.  The 
most  difficult  and  yet  important  task  is  to  develop  a  better 
way  for  determining  the  Hamiltontan  matrix  elements. 
Haas  et  al.*  and  Harrison^^  '*  have  suggested  using  the 
atomic  term  values  as  the  diagonal  matrix  elements.  Our 
work'"*  has  suggested  using  a  universal  long-range  in¬ 
teraction  to  improve  the  accuracy  of  the  conduction 
bands.  Several  studies''^^  '*’^^  have  also  pointed  out  scal¬ 
ing  rules  of  the  matnx  elements.  A  combination  of  these 
ideas  may  lead  to  an  acceptable  model.  Secondly,  both 
the  bonding  and  deep-level  states  of  impurities  should  be 
studied  at  the  same  time  in  order  to  provide  correlated  in¬ 
formation  for  defect  identification.  Finally,  more  realistic 
models  should  be  examined.  Besides  the  substitutional 
site-diagonal  defects,  one  should  consider  the  possibility 
of  interstitial,  paired,  and  even  more  complex  defects. 
One  also  needs  to  deal  with  long-range  impurity  poten¬ 
tials,  possible  charge  shifts,  and  lattice  distortions.  Pro¬ 
gress  in  all  these  areas  can  be  expected  if  the  calculation  is 
constantly  correlated  with  experiments  and  available 
ab  iniiio  theory. 
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The  bond  length  and  energy  changes  of  the  constituents  of  alloys  relative  to  their  pure  crystal 
values  are  calculated  from  an  extension  of  Harrison's  method.  It  is  demonstrated  that  the  already 
weak  HgTe  bonds  are  destabilized  by  adjacent  CdTe.  HgS.  or  HgSe,  but  are  stabilized  by  ZnTe.  It 
is  also  argued  that  dislocation  energies  and  the  hardness  of  semiconductors  vary  as  a  high  inverse 
power  of  the  bond  length  of  the  constituents.  Hence,  the  shorter  ZnTe  bond  as  an  additive  should 
improve  the  structural  properties  of  HgTe  and  CdTe.  Experiments  that  support  these  predictions 
are  noted.  The  electronic  transport  propenies  of  0.1  eV  band  gap  HgZnTe  are  about  the  same  as 
those  of  HgCdTe,  and  the  structural  properties  of  the  Zn  compound  are  superior;  thus,  we 
conclude  that  HgZnTe  is  likely  to  be  the  better  material  for  IR  devices. 


I.  INTRODUCTION 

The  objective  of  this  work  is  to  understand  the  microscopic 
mechanisms  that  govern  the  stability  of  Hg,  .,Cd,Te  al¬ 
loys,  and  then  to  suggest  changes  in  the  material  that  im¬ 
prove  Its  strength  without  adversely  affecting  its  electronic 
behavior.  In  pursuit  of  this  goal  we  have  extended  Harrison's 
bonding  theory,  ‘  which  is  applicable  to  all  tetrahedral  struc¬ 
tured  semiconductors,  to  calculate  bond  length  and  energy 
changes  in  an  alloy — including  charge  shift  and  reconstruc¬ 
tion  effects — relative  to  their  pure  crystal  values,  vacancy 
formation  energies,  dislocation  energies,  and  hardness. 

In  this  paper,  we  review  the  experimental  situation  in  sev¬ 
eral  of  these  areas  and  compare  some  of  the  results  with 
theory  for  all  the  group  IV,  IH-V  compound,  and  II-VI 
compound  semiconductors  and  their  alloys.  For  example, 
the  theory  properly  predicts  the  observed  inverse  ninth  pow¬ 
er  bond  length  id  dependence  of  the  hardness"  of  semi¬ 
conductors  and,  with  no  adjustable  parameters,  their  correct 
magnitude.  We  demonstrate  that  the  weak  HgTe  bond  is 
destabilized  by  ai'oying  it  with  CdTe,  HgS,  or  HgSe;  how¬ 
ever,  the  bond  is  stabilized  by  ZnTe.  Moreover,  because  the 
bond  length  of  ZnTe  12.406  .\\  is  14%  shorter  than  that  of 
HgTe  (2.797  .A)  or  CdTe  1 2. 804  A  I,  the  dislocation  energy  per 
unit  length  and  hardness  of  the  alloys  Hg,  _,Zn,Te  and 
Cd,  _,Zn,Te  are  predicted  in  agreement  with  some  experi¬ 
ments  to  be  significantly  higher  than  those  of  the  compounds 
with  .t  =  0.  Measuremei  :s  indicate  that  the  electron  and 
hole  mooilities  of  Hg,  _  ,Zn,Te  with  .t~0.  16,  correspond¬ 
ing  to  a  0.1  eV  band  gap.  are  comparable  to  those  of 
Hg,  _  ,Cd,Te  for  x  =  0.2.  '  This  is  to  be  expected,  because 
the  electron  effective  mass  in  a  narrow-gap  material  is  small; 
hence,  the  electron  wave  functions  at  ^ IC  are  distnbuted 
over  many  atoms,  with  the  result  that  alloy  scattering  rates 
are  small  at  the  band  edge.  The  hole  mobilities  are  large  for  a 
somewhat  different  reason:  The  valence  band  edge  is  mostly 
composed  of  Tep  states;  however,  there  is  also  some  cation  p 
state  contnbution.  Because  thep  state  energies  of  Hg  and  Zn 


are  nearly  the  same,  the  alloy  scattering  strength  at  the  va¬ 
lence  band  edge  is  small.  We  conclude  that,  because  of  its 
structural  advantage.  HgZnTe  may  be  superior  to  HgCdTe 
for  infrared  devices 

II.  ALLOY  MODIFICATION  OF  THE  BOND  ENERGY 
AND  LENGTH 

A.  Pure  compounds 

The  first  task  is  to  recast  Harrison’s  bond  energy  and  bond 
length  formalism  into  a  structure  suitable  for  generalization 
to  alloys.  Focus  attention  on  the  «th  bond  of  a  pure  zinc 
blende  structured  compound  semiconductor:  for  the  time 
being,  suppress  any  notation  identifying  it.  Then,  in  a  II-VI 
compound,  the  bond  e.nergy  is': 

+  2  +  01 

1  =  0  d 

where  the  first  term  is  the  energy  per  bond  needed  to  transfer 
two  electrons  from  a p  state  on  the  group  VI  atom  (the  anion) 
to  a  p  state  on  the  group  II  atom  (the  cation),  so  that  both 
start  with  four  electrons  (ultimately  in  the  final  bonding  ar¬ 
rangement  there  is  a  net  electron  transfer  from  the  cation  to 
the  anioni,  the  second  term  is  the  promotion  energy  per  bond 
to  form  ip'  hybrids  on  both  atoms,  the  third  term  is  the  bond 
formation  energy  owing  to  the  covalent  and  icnic  terms,  the 
fourth  term  is  Harrison’s  metallization  energy,  and  the  fifth 
term  is  the  repulsive  overlap  energy  that  prevents  the  bonds 
from  collapsing.  The  various  symbols  are  defined  below: 

. ,  is  the  (ii';""}  p  state  correlated  term  values,'* 

’’'m  - 1  ,  1  -  -  I  !/•♦  are  hopping  integrals 

between  two  adjacent  bonds  coupled  through  a 
f„  .  I  are  the  x  state  correlated  term  values, 

—  24.5/tf  "(eV)  is  the  covalent  energy, 
d  is  the  bond  length. 
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e  , ,  =if,|  ,  ,  .'tp,  .  i/4  are  the  sp'  hybnd  energies, 

tt,  1  =  ^  (L’;,  p .  i-'-.,  p ,  ,  i/(f's  -  f.' 

IS  the  contribution  to  the  reduction  of  the  bonding  state  ener¬ 
gy  fp  calculated  by  second-order  pierturbatton  theory  of  the 
bond  in  question  idenoted  by  superscnpt  ol,  owing  to  its  in¬ 
teraction  with  the  three  neighbonng  (/  =  1-3)  antibonding 
slates  shanng  the  same  {j’"rn  } 

=^ii  -ra^/:  =  . 

are  the  probability  amplitudes  of  finding  the  cation  and  an¬ 
ion  rp'-hybrid  contributions  to  the  bonding  and  antibonding 
states,  e  g..  Uf,  _  is  the  probability  amplitude  of  finding  the 
cation’  hybnd  in  the  bonding  state,  etc.. 
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T.\«i  h'  !.  CorTeliiicri  term  values  and  Htmd  energies  for  three  Il-Vl  com¬ 
pounds  Theespenmental  Kmd  iencths.  polaniv.  and  metallic  confnhuti<m 
lo  fhe  h<ind  energy  are  also  plotted 


.Atom 

Correlated  term 

( 

values  leVi 

P 

'P' 

Zn 

-  10 

- 

-  4  120 

-  6  246 

CJ 

-  1611 

- 

-  4  ’S4 

-  5111 

Hg 

-  10.146 

- 

-  4.^72 

-  6.111 

Te 

-  11  620 

- 

-  9.S24 

-  12  273 

Bond  energies 

Compound 

(7^ 

£’.,eVl 

ZnTe 

2.64.1 

0  652 

-  0.625 

-  0.185 

CJTd 

2.S05 

0.710 

-  0  466 

-  0/^44 

HgTe 

2.7^7 

0bS5 

-  0.614 

-  0.480 

'  Reference  12. 


=  is  the  polarity, 

« ,  ==  e.  P';  4-  I/^  is  the  covaiency. 


are  the  bonding  and  antibonding  one  electron  energy  states 
owing  only  to  covalent  and  ionic  interactions, 

are  the  contnbutions  to  the  bond  in  question  arising  from  its 
antibonding  state  interactions  with  a  neighbonng  bonding 
state  I.  and  •/  is  an  adjustable  coefficient  chosen  so  the  mea¬ 
sured  bond  lengths  are  reproduced. 

Hamson  has  wntten  Eq.  i  1)  without  e.xplicit  reference 
to  the  surrounding  bonds.  Equation  j  1 )  in  its  present  form  is 
rather  easily  used  for  alloy  calculations  by  replacing  host 
atom  terms  by  one  appropnate  to  the  alloy. 

The  parameter  y  for  each  compound  is  fi.xed  by  the 
condition  that  there  is  no  net  force  on  atoms  at  their  equilib¬ 
rium  lattice  positions.  This  condition  is  satisfied  if  the  gradi¬ 
ent  of  the  total  energy  with  respect  to  all  the  bond  lengths 
vanishes: 


V,^£^=0.  where  fr  =  V 

The  condition  is  equivalent  to  the  requirement  that  the  net 
bond  tension  for  the  nth  bond  equals  the  repul¬ 

sive  force 

-  T-,,  -  ~  -  0,  i2l 

where 

r,r^— I  -  :.T;  - TYi  =  4<i..TrrTT/tf.  (3i 

dd,  N  .  . \  . 

V  ,  -C  J, 

cd^  r'^1 

V  If'.,.  ,  ^  r,  ,1,  |41 


The  last  identity  in  Eq.  (41  defines  the  indicated  bond  ten¬ 
sions.  Collecting  all  the  results  from  the  earlier  paragraphs, 
one  finds  (dropping  the  bond  designation  nV. 


T1-  =  -  — 
d 


1 1  -  i«j. 


f'i  f'i 


c:  -  e  l 


T'  -  -T-  ■  =  T'  ■ 

*  m  *  rrt  ‘  m  • 


r;  =- 

d 


a-,  V-.  ^  P- 
'  1  -  or,  l«, - ^ - 

el 

ct\  V-  -  K- 
ey  -  e. 


el,- 


(5bi 


c:  .  (5ci 


ri  =Tl  =  Tl 


!  1  -1-  a  ,  Irt,  - 


fi:,  C:  +  p-> 


-  c  . 

Finally  combining  Eqs  the  bond  energy  is 

=  -  1(1  +Q'Ylf':  +  ' 

9  ,  ,  (P'l  .  +  Fj  _  I 

- 1 1  -  or:) 


(F-:  -y  P-^' 


(5dl 


-f-  ,  F,  ,  a-  P'l  _  i  -I — 2 — _ — 2_.  (6) 

Equation  l6i  was  first  denved  by  Harrison  using  a  much 
simpler  argument. '  The  bond  energies  quoted  in  Table  I  and 
on  Fig.  5  are  calculated  from  this  expression,  using  as  input 
the  correlated  term  values  from  Table  I  and  the  expenmental 
bond  lengths. 


B.  Alloy  calculations 

There  are  several  layers  of  sophistication  that  can  be  used 
to  compute  the  bond  energy  and  length  changes  between 
pure  compounds  and  their  allovs  H.mever,  all  the  calcula¬ 
tions  are  motivated  by  the  E.X.-VFS  expenmental  result 
found  first  t'or  Ga,  ,  In,  As  and  later  in  other  compounds.' 
It  vvas  found  that,  while  the  average  bond  length  in  the  alloy 
follows  Vegard's  rule  and  varies  linearly  between  those  of  the 
pure  constituents,  the  individual  bond  lengths  of  GaAs  and 
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InAs  are  nearly  the  same  as  those  of  the  pure  materials — 
changing  by  only  about  4%.  The  shorter  GaAs  bond  in¬ 
creases  Its  length  when  it  is  an  impurity  in  an  InAs  host  by 
about  4^-.  and  the  longer  InAs  bond  shrinks  by  about  4%  in 
the  opposite  extreme  where  it  is  the  impurity  in  GaAs.  The 
lattice  tits  together  by  having  the  metal  sublattice  nearly  re¬ 
tain  Its  fee  structure  with  a  lattice  constant  following  Ve- 
gard’s  rule,  and  each  As  atom  adjusts  its  local  position  and 
accommodates  to  which  metal  atoms  occupy  its  four  neigh- 
bonng  sites.  Thus  there  are  five  local  arrangements;  in  the 
first  an  .\s  atom  is  surrounded  by  four  Ga  atoms  [a 
Gai4iInioi  configuration],  in  the  second  there  are  three  Ga 
atoms  and  one  In  [a  Gai3lln|li  configuration],  and  the  last 
three  configurations  are  Gai2llni21,  Gailllniol.  and 
Gala!lni4l.  In  a  general  cation  A,  ^  alloy  or  anion 

C.-\|  alloy,  the  generalizations  are  obvious,  and  one 

can  examine  .A;4  —  rjjBlrfi:?]  =  1, 2. 3. 4  configurations. 


1.  Method  I 


The  lowest  level  of  approximation  that  takes  account  of 
these  findings  calculates  [from  Eqs.  |1)  and  '21]  a  given 
bond's  energy  and  length  shift  when  its  surrounding  bonds 
have  different  constituents.  The  calculation  is  done  assum¬ 
ing  that  the  surrounding  bonds  retain  their  pure  crystal  bond 
lengths.  This  model  is  somewhat  unrealistic,  in  that  it  ig¬ 
nores  the  strains  that  must  be  present  to  allow  the  lattice  to 
fit  together.  However,  it  is  nonetheless  useful  because  it 
treats  the  largc'^t  terms,  viz.  those  arising  from  charge  trans¬ 
fer  between  the  substituted  species  in  the  alloy,  and  results  in 
analytic  expressions,  w  hile  the  successively  more  complicat¬ 
ed  treatments  must  be  done  numerically.  It  is  also  useful 
later  because,  taken  together  with  the  more  complete  calcu¬ 
lations.  ;i  all'Mvs  one  to  judge  the  relative  contribution  from 
dilTerent  ph;>  steal  effects.  From  Eq.  1 11  in  a  cation-substitut¬ 
ed  allot  one  can  deduce  an  expression  for  the  energy  shift  of 
an  .AC  .bond  w  hen  7  surroundinz  bonds  are  replaced  bv  BC. 
JiE^as 


I 

n 


2[f-:,  u\_  iBCir,,]- 

2[r^  c:.  (ACIH,,]- 
€‘1  -  ejlACl 
J  [C-;.  c:.!BC  ly 
"  4c’J  i  fjl-eiiBC) 

_  lAor,,]- 

-  f- i.ACi 

[Ul  r,  ]- 

[cLiACT'j  r,_]- 
filACi-rj 


|7| 


where  the  p.irtial  dentatite  is  taken  with  respect  to  the  cen¬ 
tral  bond  length  'a  term  designated  with  a  superscript  01  and 
in  this  simple  approximation 

=  7  =  0.1, 2.3.  18) 


•I  AIJIKII  CONflCUHATiOII 


Id  All'Ittt  tMBEODfO  >*  A  VIttTUAI.  CflVfTAt 

Fig.  I.  Schematic  view  of  vanous  bondconlieurations.  fal  A(3lB(  11  configu¬ 
ration.  .bi  .AiJlBii!  coniiguration.  and  ici  .A|3iB;  ll  contiguraiion  embedded 
in  a  virtual  crystal. 


The  bond  length  shifts  are  found  from  Eq.  (21.  The  modifica¬ 
tion  of  these  expressions  for  an  anion-substituted  alloy  fol¬ 
lows  from  symmetry. 

2.  Method  2 

A  more  realistic  calculation  includes  strain  and  permits 
the  lattice  to  fit  together.  The  simplest  non-self-consistent, 
near-neighbor  version  of  this  calculation  proceeds  as  fol¬ 
lows.  In  a  cation-substituted  alloy  A,,,B,  C 
(e.g.  Hg,  _  ,Cd,Te  or  Zn,.,Hg,Te!,  the  cations  are  as¬ 
sumed  to  occupy  their  regular  fee  sublattice  sites,  while  the 
anions  accommodate  to  their  local  configuration.  The  aver¬ 
age  lattice  constant  is  taken  to  follow  Vegard’s  rule; 

d  =  i  1  -  .rid,,  I  AC)  -h  xff„(BC),  (9) 

where  and  d„iBCl  are  the  pure  crystal  lattice  constants. 
Then  in  an  .Ai3iBil)  configuration  [see  Fig.  llai],  the  anion 
will  be  displaced  along  the  BC  bond.  For  definiteness,  sup¬ 
pose;/  i.AC.  S'  J  BCl;  then  the  anion  will  be  shifted  away  from 
the  tetrahedral  site  toward  the  B  atom.  Because 
diuACl'-d  thed  (BC)  bond  will  be  elongated  and  the 

d  (ACl  bonds  will  also  be  locally  elongated  in  a  strained  con¬ 
figuration.  We  can  for  this  case  [see  Fig,  l(ai]  write 

diBCi  =  (l-(5d  (10) 

and  for  the  three  .AC  bonds 

d!AC)=  —  (11) 

A/  3 
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In  an  A(2)B(2)  configuration  [see  Fig.  1(b)],  we  have  for  both 
BC  bonds 

diBC)  =  -jS  =  (12) 

and  for  both  .\C  bonds 

d{Aa=^\  +y<5  +-jS-d.  (13) 

The  values  of  <5  are  determined  by  the  condition  dEp/ 
35  =  0,  or 


Equation  1 14)  is  solved  by 

(II  A.ssigning  the  their  pure  crystal  values. 

1 2)  Supposing  that  the  four  tetrahedrally  configured  an¬ 
ions  are  embedded  in  the  virtual  crystal  [see  Fig.  l(c|]. 
(3 1  For  each  x.  S  is  vaned  and  the  bond  tensions  are 
calculated  and  inserted  into  Eq.  ( 14)  until  a  value  of  5  is 
found  that  satisfies  this  equation. 

The  computer  outputs  are  5  and  the  AC  and  BC  bond  ener¬ 
gies.  This  method  can  be  extended  to  take  explicit  account  of 
the  bond  length  variations  of  larger  clusters,  and  d  can  be 
calculated  explicitly  (instead  of  using  Vegard's  rule)  so  the 
results  can  be  made  to  be  seif-consistent.  Also,  effects  arising 
from  bond  angle  distortions  can  be  included.  The  major  ef¬ 
fects  that  are  currently  included  in  the  calculation  arise  from 
the  bond  length  dependence  of  F,,  and  charge  shifts  driven 


CONFIGURAIION 

CdT* 

_\ _ >4^  MgtJICtflll 


Fit,  2.  Bonu  .cngths  js  a  fundion  of  concentration  x  for  the  Me,  ,Cd.Tc 
alloy  in  difT-.rent  configurations  The  short  dashed  lines  are  tor  the  Cdle 
btinds  and  the  Iong->hori  dashed  lines  represent  the  behavior  rt  the  Hglc 
bonds 


0  3]  as  n  as  oa  oi  as  lo 

Hglr  I 


FtG.  3.  Bond  energy  as  a  function  of  concentration  .t  for  the  Hg(  ,  ,Cd,Te 
alloy  in  different  configurations.  The  short  dashed  lines  are  for  the  CJTe 
bonds  and  the  long-short  dashed  lines  represent  the  behavior  of  the  HeTe 
bonds.  The  solid  lines  are  the  concentration  weighted  averages  of  the  respec¬ 
tive  bona  energies. 


by  hybrid  energy  differences  and  coupled  through  the  metal¬ 
lization  terms.  However.  F,  is  also  sensitive  to  bond  angle 
distortions,  and  these  eifects  are  not  included  in  the  present 
results.  We  expect  their  inclusion  will  modify  the  quantita¬ 
tive  results  by  10%-30%,  but  not  the  trends. 

C.  Results 

The  correlated  atomic  term  values'*  for  the  elements  in  the 
more  common  II-Vl  compounds  are  listed  in  Table  I,  along 
with  the  bond  energies  of  their  nine  compounds.  The  valence 
s  state  energies  for  Zn  are  large,  become  smaller  as  for  Cd, 
and  then  (in  an  unusual  occurrence!  the  trend  reverses  and 
the  Hg  5  levels  are  deeper  again.  This  is  caused  by  the  relativ¬ 
istic  terms  that  become  imponant  in  Hg.  The  p  state  levels 
have  a  similar  trend;  however,  the  differences  among  the 
atoms  are  much  smaller.  The  levels  of  the  anions  are  deeper 
than  those  of  the  cations.  Thus,  there  is  a  net  electron  trans¬ 
fer  from  the  cations  to  the  anions,  which  is  responsible  for 
the  ionic  contribution  to  the  bonding  The  polanty.  also  list¬ 
ed  in  Table  I.  reflects  the  relative  coniribution  of  ionic  and 
covalent  character  to  the  bonds.  Notice  that  is  small  for 
ZnTe  larger  for  CdTe.  and  smaller  again  for  HgTe.  By  con¬ 
trast,  the  metallic  contnbution  is  large  for  ZiiTe  and 
HgTe  and  small  for  CdTe  In  fact.  HgTe  would  not  be  bound 
if  It  were  not  for  the  metallic  contribution  to  the  bond,  not  a 
terribly  surprisi-.g  result,  because  HgTe  is  a  semirnetal.  The 
bond  lengtlis  ot  CdTe  and  HgTe  are  very  nearly  the  same, 
hut  this  IS  an  accident  resultinc  from  the  balance  of  contribu- 
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tions  to  their  bonds;  CdTe  is  dominated  by  the  covalent  and 
ionic  contnbutions  to  its  bond,  while  HgTe  has  bonds  that 
are  mostly  covalent  and  metallic.  We  should  therefore  not  be 
too  surprised  if  their  alloy  displays  some  unusual  bond 
length  variations. 

Figures  2  and  3  illustrate  the  Hg,  _,Cd,Te  alloy  concen¬ 
tration  vanations  of  the  length  and  energy  for  the  five  config¬ 
urations.  These  results  were  calculated  using  the  second 
method  discus.sed  in  Sec.  II  B.  In  Fig.  2  we  see  the  unexpect¬ 
ed  result  that  the  shorter  HgTe  bond  becomes  still  shorter  in 
the  alloy,  while  the  CdTe  bond  length  increases.  Moreover, 
the  changes  are  large  compared  to  the  original  length  differ¬ 
ence.  Thus,  one  can  expect  to  find  local  microstrains,  even  in 
this  nearly  lattice-matched  material.  The  bond  energy  varia¬ 
tions  shown  in  Fig.  3  also  display  an  important  result.  The 
already  weak  HgTe  bond  is  destabilized  by  the  presence  of 
Cd.  The  flat  concentration  variation  of  the  bond  lengths  and 
energies  is  peculiar  to  HgCdTe;  much  more  structured  beha¬ 
viors  are  exhibited  by  other  compiounds. 

Several  considerations  to  keep  in  mind  when  dealing  with 
alloys  are  depicted  schematically  in  Fig.  4  for  a  Hgi3ICd(l) 
configuration.  Due  to  the  charge  shift  alone,  the  HgTe  bond 
shnnks  by  0.030  A,  the  minimum  in  the  configuration  dia¬ 
gram  moves  to  a  smaller  bond  length  and  the  depth  of  the 
w  ell  is  decreased.  Then  the  bond  is  stretched  by  0.01 1  A  to  its 
final  length.  Thus,  the  bond  is  not  centered  at  a  local  mini- 


F«t  4  Svhe.iiiittc  picture;  Ilf  j  pntcniial  ..-onfiguraticn  Cupram  f.ir  the  HgTe 
and  CdTe  Smds  in  .in  Hg!.'  Cd  li  configuration.  For  the  HgTe  tn'nd  the 
cner  ei  and  poMtion  of  the  bond  center  are  calculated  by  method  I  of  Sec. 
If  I)  and  the  final  stretched  position  and  energy  by  method  2.  The  results 
quoted  for  the  CdTe  bond  are  those  determined  hy  method  2. 
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Fig.  5.  a.'  bond  energy  shifts  for  the  II-VI  compounds  in  All )B(3I configu¬ 
rations  e.xpressed  as  a  percentage  of  the  pure  crystal  bond  energy.  The  B 
atom  IS  designated  in  the  line  that  shows  the  shift.  The  pure  crystal  AC  bond 
energies  and  lengths  are  given  at  the  top  of  the  figure.  These  numbers  were 
calculated  by  method  1  in  Sec.  II  B 

mum.  For  e.xample.  this  implies  that  the  force  constant  of 
local  phonon  modes  will  be  related  to  the  second  derivative 
of  this  configuration  diagram  evaluated  at  the  strained  posi¬ 
tion  away  from  the  well  minimum.  This  effect  is  even  more 
pronounced  for  matenals  with  a  larger  lattice  mismatch. 
Notice  that  the  CdTe  bond  is  also  stretched  relative  to  its 
local  minimum  in  this  configuration. 

The  result  for  bond  energy  variations  for  all  the  II-VI 
compound  alloys  is  gathered  in  Fig.  5.  Here  the  energy  shifts 
of  an  .AC  bond  in  a. AillB(3)  configuration  =  3J£i  are 

presented  as  a  percentage  of  the  pure  material  bond  energy 
for  cation  substitute  alloys  along  with  equivalent  results  for 
anion  substantial  alloys,  e.g.,HgS,  _  ,Te, .  The  energy  shifts 
are  calculated  following  the  first  method  presented  in  Sec. 
II  B.  We  have  not  had  the  opportunity  to  run  all  the  com¬ 
pounds  following  the  more  sophisticated  calculational  pro¬ 
cedure.  The  bond  at  the  top  of  a  column  is  the  one  whose 
energy  shift  is  calculated,  and  the  other  constituent  is  identi¬ 
fied  in  the  line  that  indicates  the  energy  shift.  A  negative 
iAE^/Ef,  is  destabilizing  and  a  positive  one  stabilizing.  Fo¬ 
cus  on  HgTe.  It  is  destabilized  by  CdTe.  HgS.  and  HgSe.  but 
is  stabilized  by  ZnTe.  The  CdTe  destabilizes  the  HgTe  bond 
by  18''r  according  to  the  simple  calculation  and  by  15ri> 
according  to  the  better  one.  The  dilference  between  Zn  and 
Cd  occurs  because  in  the  Cd  case  there  is  a  net  electron  trans¬ 
fer  from  the  Cd  to  the  Hg,  w  hile  in  the  case  of  Zn  there  is 
practically  no  such  transfer.  These  extra  transferred  elec¬ 
trons  must  occupy  antibonding  states,  because  the  bonding 
states  on  the  HgTe  are  full,  so  the  net  efl'ect  (including  other 
energy  state  shiftsi  is  to  destabilize  the  HgTe  bond.  This  pre¬ 
dicted  destabilization  of  the  HgTe  bond  in  HgCdTe  alloys 
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Fig.  6.  Theoretical  vs  e.xpcrimentai  hardness.  The  two  values  at  the  top  and 
bottom  of  the  vertical  lines  designated  for  each  compound  were  calculated 
for  two  different  possible  dislocation  arrangements.  The  lower  one  corre- 
spends  to  the  minimum  energy  configuration  and  the  upper  one  to  a  likely 
higher  energy  configuration.  The  experimental  numbers  were  taken  from 
Ref.  2. 

agrees  with  the  d  state  energy  shift  reported  by  Spicer  era/,  at 
last  year's  workshop,'’  and  by  the  o.xygen  uptake  experi¬ 
ments  presented  by  the  same  group  at  this  conference.’ 

III.  DISLOCATION  ENERGIES  AND  HARDNESS 

The  dislocation  energies  and  hardness  of  semiconductor 
are  calculated  by  an  extension  of  Harrison's  method.'  It  is 
demonstrated,  in  agreement  with  expenment,  that  disloca¬ 
tion  energies  per  unit  length,  depending  on  the  covalency, 
are  proportional  tot/  where t/  is  the  bond  length  and  hard¬ 
ness  is  proportional  to  ~  ‘  to  t/  ”  " .  The  low  powers  are  for 
pure  covalent  matenals.  The  hardness  is  related  to  the  inter¬ 
action  energies  among  dislocations.  The  detailed  theory  will 
be  published  elsewhere." 

Figure  6  shows  the  theoretical  hardness  of  several  semi- 
condu  tors  plotted  against  experimental  values."  The  theory 
is  calculated  for  two  different  dislocation  configurations 
caused  by  the  extra  material  pushed  by  the  indentor  into  the 
crystal.  The  lower  end  of  each  vertical  line  is  a  realistic  hard¬ 
ness  for  the  minimum  energy  configuration,  and  the  upper 
ends  of  the  lines  are  the  values  found  for  a  higher-energy 
configurati  n.  In  most  cases,  one  would  expect  real  configu¬ 
rations  to  have  hardnesses  lying  somewhere  between  these 
two  values,  once  the  effect  of  heat  generated  '.ty  t'ne  indention 
process  is  subtracted  from  the  experimental  numbers.  This 
theory  has  no  adjustable  parameters,  yet  it  predicts  the  right 
trends  and  magnitudes. 

IV.  NATIVE  IMPERFECTIONS  GROWN  INTO  ALLOYS 

TTiere  are  two  kinds  of  imperfections  w  e  shall  discuss:  va¬ 
cancies  and  dislocations.  At  present,  the  conclusions  are  ten¬ 
tative,  because  the  detailed  theory  is  still  incomplete.  Va- 
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cancy  formations  energies  and  the  melting  temperature  are 
both  proponional  to  the  bond  energies  of  the  constituents. 
Thus,  one  would  expect  the  vacancy  density  just  below  the 
melting  point  to  be  about  the  same  for  all  materials  with  a 
similar  coordination  number:  in  our  semiconductor  case,  the 
coordination  number  is  four.  Dislocation  energies  are  most¬ 
ly  owing  to  long-range  bond  angle  distortions’  and  accord¬ 
ingly  are  sensitive  to  the  sheer  coefficients,  which  vary 
roughly  as  r/  “  ’  for  many  compounds.  ‘  Based  on  this  general 
argument,  one  would  expect  to  find  fewer  dislocations  in 
materials  with  short  bond  lengths. 

Dislocations  appear  in  bulk  grown  materials  through  a 
sequence  of  steps. We  shall  only  discuss  the  simplest  case, 
where  there  are  no  mechanical  stresses  on  the  growing  mate¬ 
rial.  In  the  temperature  gradient  oehind  the  growth  front, 
there  is  usually  a  supersaturated  vacancy  content.  These  ex¬ 
cess  vacancies  may  diffuse  to  surfaces  where  they  do  no 
harm,  or  they  may  condense  into  vacancy  clusters  that  can 
then  organize  into  dislocation  loops.  If  a  loop  is  not  parallel 
to  the  growth  front,  it  can  climb  toward  the  front  by  absorb¬ 
ing  more  vacancies.  If  it  climbs  faster  than  the  growth  pro¬ 
gresses.  then  the  loop  will  reach  the  front.  It  can  then  spread 
and  grow  along  with  the  crystal.  Observed  from  the  growth 
surface,  such  a  loop  will  appear  to  be  two  edge  dislocations. 
If  this  dislocation  formation  process  can  be  inhibited  in  any 
of  its  steps,  the  resulting  material  will  have  fewer  disloca¬ 
tions. 

Zn  in  Zn,  _  ,Cd,Te  or  Zn,  _,Hg,Te  will  serve  this  func¬ 
tion  in  two  ways.  The  first  is  that  Zn  stabilizes  the  weaker 
bond  in  these  alloys;  accordingly,  vacancy  formation  ener¬ 
gies  in  the  alloy  should  be  larger  than  those  in  the  pure  con¬ 
stituent  CdTe  or  HgTe.  The  melting  temperature  will  also  be 
raised  somewhat;  which  effect  “wins"  is  a  matter  still  to  be 
decided.  The  second  effect  is  clear  cut:  The  shorter  bond 
lengths  of  the  ZnTe  will  introduce  stitf  struts  into  the  system 
that  will  increase  the  dislocation  energies  and  thus  should 
inhibit  dislocation  loop  formation  and  subsequent  climb. 
The  large  dislocation  density  reduction  found  by  Bell  and 
Sen  in  Zn,  _  ,Cd,Te  with  only  4*^  Zn  supports  this  specula¬ 
tion." 

V.  CONCLUSIONS 

We  have  argued  that,  in  an  alloy,  the  bond  energies  and 
bond  lengths  are  modified  from  their  pure  crystal  values. 
Bond  energies  can  be  strengthened  or  weakened  in  an  alloy 
depending  on  the  relative  sp^  state  energies  and  the  conse¬ 
quent  net  electron  shift  between  the  constituents.  For  most 
compound  the  average  bond  length  modifications  for  each 
type  of  bond  are  generally  small  compared  to  those  suggest¬ 
ed  by  Vegard's  rule,  in  agreement  with  experiment.  How¬ 
ever,  HgCdTe  is  an  exceptional  ca^e.  because  the  bond 
lengths  for  HgTe  and  CdTe  are  nearly  the  same  (by  acci- 
dentl.  The  bond  length  and  energy  shifts  also  depend  sensi¬ 
tively  on  the  local  configuration  and  the  bonds  are  locally- 
strained.  These  effects  must  be  taken  into  consideration  in  a 
proper  theory  of  alloy  vacancy  formation  energies,  phonon 
frequencies,  dislocation  energies,  etc. 

The  conclusions  most  relevant  to  the  initial  question 
posed  in  the  introduction  are  that  for  the  narrow  gap  allows: 
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Cd  destabilizes  the  already  weak  HgTe  bond  and  has 
little  etfect  on  the  dislocation  energy. 

Zn  both  stabilizes  the  HgTe  bond  and  increases  the  dis¬ 
location  energy. 

Zn  also  decreases  the  dislocation  density  of  CdTe  to 
improve  substrates. 

HgZnTe  electron  and  hole  mobilities  are  about  the 
same  as  those  of  comparable  band  gap  HgCdTe. 
HgZnTe  is  likely  to  be  a  superior  material  to  HgCdTe 
for  detectors  operating  in  the  8-12  fim  regime. 
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Abstract 

The  objective  of  this  paper  is  to  identify  the  principal  microscopic 
phenomena  controlling  dislocation  densities  in  bulk  grown  semiconductors. 

Then,  based  on  this  understanding,  a  strategy  for  selecting  materials  to  reduce 
disLocation  densities  is  offered.  The  relevant  quantities  are  calculated  from 
an  extension  of  Harrison's  bonding  theory,  which,  with  our  improved  accuracy 
relates  properties  of  the  solids  to  the  constituent  atoms*  valence  electron 
energy  states  and  wave  functions.  We  report  on  the  alloy  composition  varia¬ 
tion  of  bond  energies,  bond  lengths,  charge  redistribution  among  constituents, 
vacancy  formation  energies,  dislocation  energies,  and  hardness.  Several 
III-V  and  II-VI  compound  semiconductors  are  treated  Including,  GaAs,  GalnAs, 
HgCdTe,  and  ZnHgTe. 


(a)  This  work  was  supported  In  part  by  AFOSR  Contract  F49620-81-K0012  and 
DARPA  Contract  MDA-903-33-C-03 08 . 
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Introduction 


According  to  a  currently  accepted  model  (1)  of  the  mechanism  leading  Co 
dislocations  in  bulk-grown  material  at  a  given  temperature  and  temperature 
gradient  behind  the  growth  front;  their  density  Is  governed  by: 

•  The  vacancy  density  that  depends  on  the  formation  energies, 

•  The  competition  between  vacancy  annealing  rates  and  vacancy 
interaction  caused  clusters, 

•  The  condensation  rate  of  these  vacancy  clusters  into  dislocation 
loops,  and 

•  The  subsequent  growth  race  of  these  loops. 

The  objective  of  this  paper  is  Co  identify  some  of  the  principal  microscopic 
phenomena  controlling  these  features  in  semiconductors  compounds  and  their 
pseudobinary  alloys.  If  any  of  the  foregoing  steps  can  be  inhibited,  then 
there  will  be  fewer  dislocations  in  bulk-grown  crystals. 

A  model  of  the  bonding  of  cetrahedrally  coordinated  semiconductor  com¬ 
pounds  due  to  Harrison  underlies  this  wotk  (2).  He  has  derived  expressions 
for  the  bond  energies  and  strain  coefficients  of  pure  semiconductor  compounds 
in  terms  of  the  conscituenc  atom's  valence  state  wave  functions  and  energies. 

There  are  four  contributions  to  the  bond  energies: 

•  A  covalent  energy  arising  from  Che  interaction  between  sp^  hybrids 
on  adjacent  sites,  which,  according  to  a  universal  rule  deduced 

by  Harrison,  is  related  to  the  inverse  square  of  the  bond  length,  d, 

•  An  ionic  energy  which  is  proportional  to  the  energy  difference 
between  the  sp-  hybrid  energies  of  the  anion  and  cation, 

•  A  metallic  energy  arising  from  the  interaction  between  filled 
bonding  and  unfilled  antibonding  states  on  adjacent  bonds,  and 

•  An  overlap  repulsion  energy  which  is  taken  to  vary  as  d 

The  shear  coefficients,  which  play  an  important  role  in  dislocation  energies, 
are  shown  to  vary  for  covalent  solids  as  Che  covalent  energy  per  cell  unit 
volume  or  as  The  Ionic  energy  arising  as  it  does  from  coulomb  inter¬ 

actions,  is  insensitive  to  bond  angle  distortions  and  depends  only  on  bond 
lengths.  Thus,  in  the  ionic  semiconductors  the  shear  coefficients  are  reduced 
from  Chose  of  equivalent  bond  length  covalent  compounds,  but  also  t.he  power 
law  dependence  on  bond  length  increases  to  for  Che  extreme  ionic  limit  (3). 

The  metallic  interaction  is  bond  angle  sensitive  and  contributes  to  the  d"^^ 
power  law. 

We  have  modified  Harrison's  theory  in  several  ways  C3)  (^).  First, 

Instead  of  using  atomic  term  values  as  the  input  energies  in  the  bonding 
calculation,  we  use  correlated  atomic  energies  chat  take  account  of  all  the 
atomic  electronic  energy  changes  associated  with  a  state  change.  In  our  prior 
publications  these  energies  were  calculated  from  pseudopotencials  (3).  In 
this  paoer  we  have  modified  the  procedure  and  taken  the  ground  state  energies 
froro  free  atom  loniration  experiments,  and  only  calculate  the  excitation 
energies  from  the  pseudopotent ials .  Second,  the  Harrison's  theory  has  been 
extended  to  allovs  in  which  each  tvpe  of  tr.e  Individual  bond's  energy  and 
length  changes  and  the  net  allcv  substitution  energies  are  calculated.  A 
theory  of  dislocation  energies  and  hardness  of  semiconductors  has  been  advanced. (i) 
In  this  paper  the  second  item  will  be  emniiasiced  since  item  one  is  less 
relevant  to  this  conference  and  these  will  be  published  elsewncre. 
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Bond  Length  and  Energy  Changes  In  an  Alloy 


Single  Impurity 

Start  by  considering  Che  simplest  defect,  one  cation  of  a  host  seaicon- 
ductor  compound  AC  is  replaced  by  an  isoelectronic  impurity  I,  as  sketched 
in  Fig.  1.  This  will  be  generalized  to  an  alloy  later  ans  the  cation 

substituted  case  follows  by  syiametry.  To  study  the  principal  effect 

of  core  lattice  distortions,  we  allow 
the  positions  of  the  first  shell  of 
atoms  labeled  by  bond  length  d^  to  move, 
but  fix  the  atoms  of  the  second  neigh¬ 
bors  and  beyond  at  their  pure-crystal 
positions.  Eventually,  the  effects  of 
long  range  strain  fields,  generated  by 
a  point  distortion,  must  be  added  to 
this  core  distortion.  Here  we  allow 
only  three  different  bond  length®. 
dl  -  (l-6)d,  dj  -  (1  +  2/36  +6^)^'^d. 
and  d^  *  d,  where  d  is  the  bond  length 
of  the  pure  host  AC  compound,  and  6  is 
a  scale  parameter. 

There  are  four  bonds  with  energy 
Ej,  twelve  with  energy  E2,  and  3b  with 
energy  E3.  These  energies  can  be 
written  in  closed  form  in  the  Harrison 
theory.  The  energies  and  differ 
from  the  host  bond  energy  becaiise 
of  bond  length  and  angle  di.stortions 
as  well  as  difference  in  the  ionic 
energy  of  I  and  C  atom  Induced  charge 
transfer.  The  energy  £3  differs  from 
Eq  only  by  the  charge  transfer  coupled 
through  the  metallization  terns.  There 
are  two  Interesting  energies  we  can  calculate.  The  first  Is  the  energy  ^ 
required  to  substitute  the  I  atom  for  a  C  atom,  l.e.  bring  a  free  I  atom  from 
infinity  and  replace  a  C  atom  that  is  taken  from  the  crystal  and  removed  to 
infinity. 

ds  +  12dE2  +  36flE^  -  (c^  -  Cj)  (X) 

where  i£.  «  E.  *  E  ,  J  •  1,2,3,  and  and  are  the  free  atom  valence  elec¬ 
tron  binding  Energies.  The  second  is  the  bond  energy  change  of  the  impurity 
AI  bond  relative  to  the  bond  energy  of  a  pure  AI  compound,  denoted 

i  ~  *c^ 

This  energy  tells  us  If  the  AI  bond  is  stabilized  (fij  <  0)  or  destabilized 
(At  >  0)  when  It  Is  in  an  AC  host.  If  we  define  ^  as  the  bond  energy  dif¬ 
ference  between  the  lA  and  CA  bonds  each  in  their  respective  pure  crystals 

^  z  BE(LA)  -  BE(CA)  (3) 


Schematic  representation  of  the 
bonds  around  an  atooi  of  Interest 
(designated  by  Che  square).  The 
first,  second,  and  third  neigh¬ 
bor  bonds  are  designated  42* 
and  d^ 


then  one  can  vnrite 


(4) 


where  t  f*  extra  energy  difference  caused  by  strain  and  charge  transfer. 


The  energyes  E,  and  bond  lengths  d,,  j  =  1,  2,3,  are  calculated  by  mini¬ 
mizing  ^5.  If  we  let  dj  =  dg  (1  -  65)  where  dj  is  the  bond  length  in  the 
pure  AI  lattice,  Chen  Che  approximate  expression 


(5) 


^^3  B- 


8(Cii  -  C^,) 


27  B, 


+  (small  terms  due  to  charge  transfer) 


can  be  deduced*  The  bulk  moduli!  B  of  the  host  and  Bj  of  the  impurity  lattices 
and  a  shear  coefficient  -  C12  the  host  appear  in  the  expression.  While 
this  expression  is  approximate*  the  effects  of  the  various  types  of  strains 
can  be  visualized  using  It.  Our  detailed  calculations  are  done  using  the  full 
theory.  If  one  fixes  the  bond  lengths  »  d^  •  dj  •  dg,  then  6  *  and  the 
strain  energy  Asc  large  and  positive.  If  one  allows  only  bond  length  dis¬ 
tortions  (C]^i  -  Cj^2  “  then  6  is  reduced  and  the  bonds  tend  to  adjust  so 


dj_  ftsdj.  This  reduces  considerably  from  the  undistorted  lattice  case. 


If  dx  >  dg  Chen  both  bonds  are  stretched  somewhat,  -  d^  >  0  and  d2  -  dQ  >  0. 
However,  this  configuration  produces  large  bond  angle  distortions.  W’hen  the 
shear  coefficients  are  turned  on  the  lattice  relaxation  is  modified  and  the 
bond  lengths  cannot  adjust  as  much,  so  di  differs  from  by  a  larger  amount. 


energy  either  as  a  bond  length  or  bond  angle  distortion. 


If  we  calculate  A  from  the  full  theory,  then  we  predict  values  that  are 
too  large.  This  occurs  because  the  theory  predicts  -  Cj^2  '^hich  agree  well 
with  experiment  but  it  predicts  bulk  moduli!  with  the  proper  trends  from  one 
compound  to  mother  but  which  are  about  a  factor  of  2.4  coo  small.  If  we  use 
the  experimental  values  of  the  strain  coefficients  in  Eq,  (5)  then  good  agree¬ 
ment  is  found  with  the  experiments  (5)  on  and  2nSe».^Tex*  The 

results  are  quoted  in  Table  I. 


Table  I.  Bond 

Leneths  in 

A  for  Impurities 

In  the  Indicated 

Host 

Ga 

in  InAs 

In  in  GaAs 

Se  in  2nTe 

Te  in  Znfe 

Experiment 

2.487 

2.587 

2.1.90 

2.595 

Eq.(5)  and 
experimental  BiC 

2.499 

2.547 

2.510 

2.570 

Full  Theory 

2.538 

2.518 

2.570 

2.512 

Alloy 


In  an  3^C  alloy  the  four  cations  around  a  given  C  anion  can  be 


1— X*  -- 

arranged  in  rive  different  configurations,  denoted  by  A(4-T).  BC*"},  7]  *  0,1, 
2,3.4,  An  A(2)  B(2)  configuration,  for  example,  is  one  in  which  the  C  atom 
has  two  A  atom  and  two  B  atom  neighbors.  It  is  possible  to  solve  the  full 
alloy  problem  for  large  clusters,  but  for  now  we  have  restricted  the  cluster 
to  near  nelghoors  only,  and  in  Fig,  1  the  central  atom  is  now  an  anion  C  and 
d^  *  “  d^ff  are  taken  to  be  an  effective  medium  bond  length  that  is  deter¬ 

mined  seif- *cns istencly .  The  different  b.’nds  of  cvpe  1  no  longer  have  the 
same  length.  The  procedure  is  as  follows:  First  assume  a  value  for  d^jf, 
say  Che  virtual  crvstal  bond  lencth  (l-x>  d  f.xD  xJ  vbc).  Next  calculate 
the  di  values  for  the  ve*"ioi.s  conf  igur  it  ions  ov  m'-nimirihg  cne  energy  of  the 
con  f  icur  It  ;  ■*n.  Ih^  n,  ♦.r:.:'  ^ration  iveract  the  various  a  values  to  find  a 
new  d.,.-f.  Ihen  ,  •.•rocedure  until  it  converges. 
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Figure  2:  Bond  length  (a)  and  energy  (b)  of  Gaj^^^tn^As 
as  functions  of  x.  The  dashed  curves  are  the  average  values 
for  the  designated  bonds,  and  the  solid  curved  are  the  alloy 
averages. 

Because  Che  bond  angle  distortion  terms  are  unphysically  large  we  have  left 
them  out  of  the  present  caiculacion.  Hence  bond  length  shifts  will  be  some¬ 
what  too  soali. 

Results  for  the  bond  energy  and  bond  length  shifts  with  concentration  are 
given  for  Gaj^.^IoyVs  in  Fig.  2  and  in  rig.  3  examples  of  two 

different  behavior  patterns.  The  predicted  trends  for  CalnAs  agree  with 
experiment  but  Che  bond  Length  changes  are  coo  small  as  expected  with  Che  bond 
angle  distortion  terms  absent.  Notice  :haC  the  longer  bond  in  this  case  de¬ 
creases  In  length  and  the  shorter  one  increases  as  expected.  However,  in  the 
BIoP  case  the  charge  shift  terms  are  so  large  that  the  longer  In?  bond  has  a 
minimum  as  a  fucntion  of  composition  rather  than  a  oonoconic  decrease.  Be¬ 
cause  of  Che  large  bond  length  difference  between  BP  and  In?  there  Is  probably 
a  oLissability  gap  In  this  alloy  that  prevents  these  compositions  from  being 
prepared.  However,  one  may  be  able  to  examine  Che  variation  of  the  anomalous 
behavior  of  the  InP  bond  in  a  BP  host  (z  «  1). 

The  Hg^.^Cd^Te  system  is  completely  anamolous.  The  bond  lengths  of  CdTe 
(2.805*  A)  and  Hgte  (2.797' A)  are  nearly  the  same  by  an  accident.  CdTe  bonds 
are  dominated  by  covalent  and  ionic  interactions,  while  HgTe  is  sore  covalent 
and  the  metallic  terms  are  important.  The  mix  of  interaccions  In  the  ttro 
cases  leads  Co  Che  same  bond  lengths.  When  an  alloy  is  made  the  charge  shift 
terms  dominate,  and  they  cause  the  longer  CdTe  bond  to  become  still  longer  and 
Che  HgTe  to  contract  by  amounts  that  are  large  compared  to  the  pure  crystal 
difference.  Horeover,  the  already  weak  Hgte  has  its  bond  energy  reduced  still 
more  by  the  presence  of  Cd.  Since  the  melting  point  of  HgTe  increases  as  Cd 
is  added  and  the  local  strength  of  HgTe  bonds  adjacent  to  a  Cd  decrease, 
vacancy  densities  will  Increase.  All  these  and  other  observed  phenomena  are 
predicted  by  Che  theory.  A  complete  catalog  of  results  requires  more  pages 
than  we  have  been  allocated  la  this  article,  but  we  have  tried  to  offer  a 
repcesentati ve  group. 


Figure  3:  Bond  length  (a)  and  energy  (b)  of  as 

functions  of  x.  The  dashed  curves  are  the  ax^erage  values 
for  Che  designated  bonds  and  the  solid  curves  are  the  alloy 
averages . 

Conclusions 

We  have  denonstrated  in  agreezeot  with  experlaienc  that  dislocation  energies 
and  hardness  of  sesiconduccors  are  proportional  to  the  shear  coefficients  which 
vary  roughly  as  d”^.  Thus  to  decrease  dislocations  In  a  given  seiDiconductor 
one  should  inhibit  their  forsation  by  introducing  some  means  of  shortening 
bonds.  This  should  be  accomplished  vithout  also  decreasing  the  ratio  of  the 
vacancy  formation  energy  to  Che  melt  temperature.  While  detailed  calculations 
CO  support  the  following  contentions  are  still  incomplete*  a  sec  of  criteria 
on  an  impurity  (denoted  1)  In  a  host  semiconductor  (denoted  H)  that  are  likely 
CO  satisfy  these  conditions  are:  the  bond  length  of  the  Impurity  is  smaller 
Chan  Chat  of  Che  host  d^  <  d^»  the  smaller  bond  energy  Is  stabilised  in  the 
alloy  min  |AEbH*  1  ^  these  circumstances,  the  average  bond  length 

will  shrlnx  and  in  tne  vicinity  of  each  impurity  the  four  surrounding  bonds 
nearly  have  the  length  of  the  Impurity  but  are  stretched  slightly.  The  next 
neighboring  host  bonds  are  also  stretched.  The  net  effect  is  an  arrangement 
t/hich  Is  more  rigid  than  the  unperturbed  lattice  and  consequently  the  local 
shear  strain  energy  increases*  causing  dislocation  energies  to  Increase. 

Examples  of  this  case  are  Bi^^Oa^As,  and  Zn^.^Hg^Te.  Both  B  in  GaAs  (6)  and 
Zn  in  HgTe  (7)  have  proven  to  be  effective  in  reducing  dislocation  densities. 
Unfortunately*  B  has  a  low  solubility  in  GoAs  and  it  is  not  clear  chat  enough 
can  be  gotten  in  to  make  it  easy  to  prepare  dislocation-free  material. 

A  second  case  where  an  improvement  occurs  Is  if  dj  >  dj^,  and  again 
min  ^  this  case  each  of  the  four  longer  impurity  atom 

bonds  are  cocsprecsed  by  the  surrounding  host  bonds  and*  tnore  importantly,  the 
twelve  next  neighbor  host  bonds  arc  also  compressed.  Once  again  for  sowill 
Impurity  concentrations  the  net  dislocation  energy  should  be  increased.  How¬ 
ever,  in  this  case  Che  cfiect  Is  ct-nnpeting  erclnst  a  net  bend  lengthening 
trend  of  the  alloy  which  Is  tending  to  make  it  less  rigid  against  a  shear,  so 
at  higher  concentrations  the  occhaaism  should  cease  to  function.  An  example 


of  this  case  Is  Gai«^InyA.9,  where  In  has  proven  to  be  effective  in  dislocation 
reduction  of  GaAs  (8). 

We  have  demonstrated  that  our  modification  of  the  Harrison  bonding  theory 
accruately  predicts  the  observed  change  of  bond  lengths  in  semiconductor  alloys 
and  offers  guidance  to  coeans  for  dislocation  reduction. 
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ABSTIL\CT 

Several  recent  theoretical  studies  of  the  local  structure  of  semiconductor  alloys  are 
summarized.  First,  dilute  limit  calculations  of  local  bond  lengths  and  mixing  enthalpies 
are  discussed.  These  calculations  include  effects  due  to  both  bond  length  and  bond- 
angle  distortions,  as  well  as  local  chemical  rearrangements.  Then,  a  new  statistical 
theory  of  concentrated  alloys  is  described.  Deviations  from  random  alloy  distributions 
(microclusters)  are  predicted. 


INTRODUCTION 

This  paper  summarizes  our  recent  theoretical  studies  directed  toward  understand¬ 
ing  the  microscopic  structures  of  pseudo-binary  semiconductor  alloys  A^jOi.^C.  \V"e  first 
present  a  detailed  calculation  of  the  local  bond  length  relaxation  in  the  dilute  limit 
X  — •  0,  i.e.  the  case  where  an  A  atom  is  substituted  for  a  B  atom  in  a  BC  compound. 
The  mixing  enthalpy  parameter  Q  is  found  to  be  related  to  small  excess  substitution 
energies.  These  excess  energies  are  calculated  directly  through  a  minimization  pro¬ 
cedure.  Thus,  the  accuracy  of  the  predicted  0  is  not  limited  by  trying  to  find  small 
differences  between  large  numbers.  The  theory  is  then  generalized  to  concentrated 
alloys  using  statistics  based  on  combinations  of  tetrahedral  clusters  of  five  atoms.  Our 
results  predict  that  microclustering  occurs  in  a  majority  of  alloys.  We  conclude  by  iden¬ 
tifying  systematic  correlations  between  the  theory  and  several  experiments. 

Before  discussing  the  calculation,  it  is  useful  to  provide  some  background  about 
the  structure  of  these  alloys.  It  was  customary  to  assume  that  these  alloys  have  two 
sublattices  in  which  the  C  atoms  occupy  one  sublattice,  and  A  and  B  atoms  are  ran¬ 
domly  distributed  on  the  other.  This  picture,  referred  to  as  the  virtual  crystal  approxi¬ 
mation  (\'CA),  implies  that  the  nearest-neighbor  (nn)  bond  lengths  in  the  alloy  are  the 
concentration  weighed  average  values,  i.e.  d^c  =  dge  =  d  =  x  +  (1-x)  d^^  where 
the  values  with  a  superscript  (0)  denote  the  pure-crystal  values.  On  the  other  hand, 
according  to  Pauling’s  covalent  radii  approximation  (CRA),  the  local  bond  lengths 
retain  their  respective  pure-crystal  values,  i.e.  d^(-  =  d|^  and 

*^BC  — 

If  we  define  =  (d  -  dy(^)  /  d  and  5  =  (d  -  d^^)  /  d,  then  the  ratio  6/6q  in  VCA 
is  zero,  but  in  CR.\  it  is  1.  However,  Mikkelsen  and  Boyce^'^  found  from  their  EXAFS 
experiment  on  Gaxlni.x-'^^  nn  bond  lengths  do  not  fit  either  VC.-\  or  CRA. 

Instead,  they  found  the  value  of  S/Sq  to  be  close  to  3/-4.  Since  then,  similar  exporimeiiis 
have  been  done  for  a  number  of  zinc-blende  pseudo-binary  alloys,*''  and  the  3/1  rule 
ap[)cars  to  be  cpiite  general. 
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DILUTE  LIMIT 


The  dilute  limit  is  the  easiest  case  but  is  still  not  trivial.  Its  solution  provides 
both  end-point  results  (x  =  0  and  1),  as  well  as  insight  into  the  extension  to  the  concen¬ 
trated  alloy  case.  \  complete  description  of  this  case  is  being  reported  elsewhere;  here 
we  summarize  the  essential  results.  The  substitution  energy  Ag  for  an  A  atom  replacing 
a  B  atom  in  a  BC  compound  is  calculated  and  minimized  to  find  the  relaxed 
configuration.  Ag  can  be  written  as  Ag  =  4  (AE^  +  AEg  +  AEgh),  where  AE(,  is  the 
binding  energy  difference  between  the  AC  and  BC  compounds,  AEg  is  the  strain  energy, 
and  AE(,h  is  a  chemical  energy  shift.  All  AE'  s  are  energies  per  bond.  Then, 
AE  =  AEg  -I-  AE(,h  is  the  excess  energy  per  bond  for  the  impurity  substitution.  AEg  is 
calculated  by  dividing  the  crystal  into  two  regions.  Outside  R  (which  is  the  distance  of 
the  second-shell  atoms  to  the  impurity),  the  distorted  crystal  is  treated  as  an  elastic  con¬ 
tinuum  with  a  radial  displacement  field  which  is  inversely  proportional  to  the  square  of 
the  radial  distance,  so  AEJ°“*'^  =  l/4RCu“,  where  C  is  an  effective  shear  coefficient, 

C  =  TT  (1.6  (Cji  -  C,2)  -f  4.8  C44),  (1) 

and  u  is  the  magnitude  of  the  displacement  at  R.  Inside  R,  the  strain  energy  aJ'")  is 
treated  with  a  valence  force  field  (VFF).^^)  Finally,  the  chemical  energy  shift  AE^h  is  cal¬ 
culated  from  Harrison’s  model  and  arises  from  changes  in  the  metalization  energies 
caused  by  different  bond  lengths  Ad  s  dgc  -  d;^c  and  covalent  energies 
AV3  =  V3(AC)  -  V3mc).  Note  that  do  =  (d^^  -  d^?))  /  d^^  and 
6  ~  (d|^^  -  dy^c)  /  d^^  in  this  dilute  limit,  so  the  excess  energy  AE  can  be  expanded  up 
to  second  order  in  6,  u,  and  AV3  For  a  given  pair  A  and  B,  AE  is  an  explicit  function  of 
8  and  u.  Minimization  of  AE  with  respect  to  5  and  U  leads  to  the  equilibrium  local 
bond  length  d^^c  and  energy  AE.  Then,  AE  is  used  to  estimate  the  mixing  enthalpy 
parameter  0  in  the  mixing  enthalpy  AH  =  x(l-x)n  by 

n  =  2  (AE  (A  in  BC)  -b  AE  (B  in  .\C)).  (2) 


A  systematic  comparison  with  other  models  based  on  strain  energy  alone  shows 
that  an  increase  of  the  range  of  the  fixed  boundary  R  increases  the  relaxation  of  d^p,  i.e. 
it  causes  S/Sq  to  increase.  The  inclusion  of  the  bond  angle  restoring  force,  on  the  other 
hand,  reduces  the  relaxation.  It  turns  out  that  a  delicate  cancellation  of  these  two 
effects  causes  a  simple  spring  model  pointed  out  by  Shih  et  al.  (SSIIS)^^*  to  yield  accu¬ 
rate  results.  In  this  model  6  =  6g  /  (1  +  1/3  q/o|),  where  a  and  «[  are  the  bond¬ 
stretching  force  constants  for  the  host  (BC)  and  the  impurity  (AC)  crystals.  W'ith 
a  Oj,  this  model  predicts  6  =  3/46g  for  a  zinc-blende  alloy.  Although  our  full  pertur¬ 
bation  theory  (FPT)  and  the  VFF  model  of  Martin  and  Zunger  (MZ)^''*  predict  d^^Q  with 
an  average  absolute  deviation  comparable  to  the  experimental  uncertainty  of  0.01,\.  the 
simple  spring  model  is  even  better. 

We  note  that  while  in  our  theory.  MZ  and  SSI  IS.  the  H  values  are  directly  calcu¬ 
lated  without  any  adjustable  parameters,  our  theory  and  SSIIS  agree  with  the  exiicri- 
ment  as  well  or  even  slightly  better  than  the  one-parameter  theories.*’®^  Although  our 
theory  predicts  a  negative  H  value  for  all  three  (Ga.  .M)  alloys,  the  magnitude 
(R  =  -0.17  kcal/mole)  is  t(K>  small  to  account  for  the  ordering  of  G:1t,.\I,_j.\s  grown  at 
600  to  700  "C  found  recently The  calculated  R  values  also  provide  guidance  in 
separating  miscible  from  immiscible  alloys.  In  a  random  alloy,  the  criterion  for  alloy 
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mixing  for  all  x  is  T  >  T^.,  where  the  critical  temperature  T  is  given  by  =  n/2Rg 
with  Rg  being  the  universal  gas  constant.  Figure  1  is  a  plot  of  T^/To  against 
I  ^0  I  /  I  i  ’  "^here  T2  is  the  lower  of  the  two  constituent’s  melting  temperatures,  and 
=  1.63x^  with  Xm  being  the  ratio  of  rms  bond  length  amplitude  fluctuation  to  the 
bond  length  at  To.  The  simple  spring  model  gives  T^/T2  =  ('^o/'^m)^  ^  indicated  by  the 
solid  curves. 

There  is  an  empirical  rule  stating  that  a  miscibility  gap  will  occur  if  |  Sq  | 
between  two  alloy  components  exceeds  7.5%.  However,  if  is  plotted  against 

I  60  1  1  •'he  simple  spring  model  would  not  exhibit  a  smooth  simple  quadratic  curve,  and 
our  theoretical  points  would  be  much  more  scattered.  This  suggests  that 
I  *^0  I  /  I  I  >  1  is  a  better  criterion  than  |  I  >  0.075.  Figure  2  also  clearly  shows 
the  chemical  effects,  namely  negative  and  positive  chemical  energies  AE,.),  for  cation  and 
anion  substitutions  respectively.  The  full  theory  and  the  experiments  correlate  within 
the  experimental  uncertainties.  The  simple  SSHS  model  clearly  is  an  excellent  universal 
representation.  However,  T,,/T2  varies  faster  than  quadratically  for  larger  [  5o/'^m  I 
values,  as  born  out  from  both  the  experimental  data  and  the  full  theory. 


CONCE.NTR.\TED  ALLOYS 

Turn  now  to  the  concentrated  alloy  case.  First,  an  improved  statistical  model  is 
required.  We  have  extended  regular  solution  theory  based  on  pair  energies  to  one  for 
five-atom  clusters.  For  an  alloy,  the  building  blocks  are  clusters  of 

A(m)B(  l-m)C,  where  m  ranges  from  0  to  4.  For  a  given  alloy  concentration  x  and  for  a 
given  set  of  energies  associated  with  these  clusters,  we  have  derived  expressions  for 
the  cluster  population  distribution  x^,  =  ITn,  /  N,  where  N  is  the  total  number  of  unit 
cells  and  n^  is  the  averaged  number  of  cells  with  A(m)B(4-m)C  clusters.  The  partition 
function  Z  is  obtained  using  a  steepest  descents  argument  which  then  yields  the  mixing 
Flelmhotz  free  energy  AF.  The  result  reduces  to  Guggenheim’s  tetrahedron  case^^)  if 
pair  potentials  (for  the  second-neighbors)  are  assumed.  Anotlier  major  difference  is  that 
we  only  need  to  solve  a  single  quartic  equation,  while  Guggenheim  needed  to  solve  four 
simultaneous  quartic  equations. 

The  key  to  the  problem,  however,  lies  in  the  calculation  of  the  energies  If  one 
assumes  that  the  size  of  the  tetrahedra  for  all  m-clusters  at  a  given  alloy  concentration 
takes  on  the  corresponding  VGA  values  but  allows  the  central  C  atom  to  relax,  then  the 
energies  as  functions  of  x  behave  like  those  shown  in  Figure  2(a).  There  are  at  least  two 
major  flaws  in  this  result.  First,  the  energies  are  too  large  and  would  correspond  to  R 
values  many  times  the  experimental  values.  Second,  at  x  =  0.75,  0.5.  and  0.2.  these 
energies  imply  compound  formation  for  A3BiC4,  A2B.2C4,  and  AiB3C^,  respectively, 
which  is  opposite  to  the  known  tendency  for  spinodal  decomposition  of  Ga^In,.,(.\s  at 
low  T.  However,  if  the  local  cell  volume  of  each  cluster  is  allowed  to  be  in  mechanical 
pquilibrium  with  a  continuous  medium  with  an  effective  shear  coefTicient 
C  =  xCl?(  +  (1-x)  C..cy,  where  the  C  value  for  the  pure  material  is  given  by  Eq.  (1), 
then,  the  corresponding  energies  a.s  a  function  of  x  are  given  in  Figure  2(b),  which 
now  yielfis  a  reasonable  value  of  mixing  enthalpy  and  correctly  predicts  the  tciidency 
toward  spinodal  decomposition  at  low  temperature.  With  this  set  of  energies,  one  can 
then  calculate  the  cluster  distribution  Xj,,,  ami  compare  them  with  the  corresponding 
values  for  a  random  alloy,  i.e.  x,jf' =  ( ln,)x'"(l -x) ' Figure  .'J  shows  the  deviation 
from  randomness  Ax^,  =  x^  -  x,jf)  as  a  function  of  x  for  four  arbitrarily  chosen  growth 
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temperatures.  It  shows  that  Ax„  can  be  several  percent  at  ordinary  growth  tempera¬ 
tures  and  there  are  deviations  from  the  symmetric  behavior  about  x  =  1/2  predicted  by 
the  pair  energy  model. 

It  is  also  interesting  to  examine  the  sets  of  bond  lengths  d^c^  correspond¬ 

ing  to  the  values  of  Figure  2(b).  The  spread  in  the  d  values  among  different  clusters 
is  found  to  lie  within  the  width  of  EXAFS  lines. 

The  existence  of  microclusters,  whose  populations  deviate  from  those  of  a  random 
alloy,  will  impact  on  many  experimental  results.  These  include  phase  diagrams,  EXAFS 
line  positions  and  widths,  magnetic  quadruple  splittings,  modulation  spectroscopy 
widths,  Raman  frequencies  and  intensities  for  different  modes,  etc.  We  have  demon¬ 
strated  that  the  theory  agrees  with  EX/VFS  and  measured  mixing  enthalpy  parameters. 
It  should  be  tested  against  more  experiments  to  further  our  understanding  of  the  struc¬ 
tural  properties  of  semiconductor  alloys. 
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SEMICONDUCTOR  PSEUDO-BINARY  ALLOYS; 

BOND-LENGTH  RELAXATION  AND  MIXING  ENTHALPIES 
A.-B.  Chen,  Physics  Department 
Auburn  University,  Alabama  368149 
A.  Sher,  SRI  International 
Menlo  Park,  California  914025 

Harrison's  bonding  theory,  the  valence  force  field  (VFF),  and  an 
elastic  continuum  are  combined  to  study  the  substitution  energies  and 
local  bond  lengths  d^  of  isoelectronic  impurities  in  semiconductors. 

Explicit  expressions  for  A  and  d  are  derived,  which  enable  us  to  absorb 

3  1 

measured  elastic  constants  into  the  calculation  and  to  study  the  chemical 
effects  arising  from  differences  in  the  covalent  radii  and  polarities. 
Several  models  based  on  VFF  alone  are  also  derived  for  comparison.  The  full 
theory  and  at  least  five  VFF  models  are  found  to  produce  impurity  bond 
lengths  in  excellent  agreement  with  experiment.  The  substitution  energies 
are  shown  to  provide  good  estimates  of  the  mixing  enthalpies  {i  of  pseudo¬ 
binary  alloys  and  to  predict  miscibility  gaps  properly.  The  chemical  shifts 
in  n  are  found  to  be  negative  for  most  cation  alloys  but  positive  for  anion 


substitutions . 


I  INTRODUCTION 


The  discovery  of  a  bimodal  distribution  of  the  nearest-neighbor  bond 
lengths^  in  Ga^In^_^As  has  sparked  considerable  interest  in  the  bonding 
nature  of  semiconductor  alloys.  This  finding  has  changed  the 
conventional  picture  of  the  alloy  crystal  bond  configuration,  which  has  far- 
reaching  implications  about  the  electronic  structure,  structural  stability, 
and  thermodynamics  of  these  materials.  Because  of  the  complexity  of  both 
the  structural  and  the  potential  disorder  in  these  alloys,  ab-initio  band- 
structure  techniques  have  not  yet  evolved  to  a  stage  suitable  for  direct 

7 

calculations.  Therefore,  we  have  extended  Harrison’s  bonding  theory  to 

5  6 

study  the  alloy  structural  properties.  ’  In  this  paper,  we  apply  an 
intermediate  version  of  the  theory  to  the  dilute-limit  case  of  an 
isoelectronic  impurity. 

A  particularly  useful  application  of  the  theory  is  its  perturbation- 

expansion  form,  in  which  measured  elastic  constants  are  incorporated  to 

obtain  accurate  results.  This  form  is  also  useful  for  comparison  with  other 

3  a  o 

previously  published  models'^’  ’  that  are  based  on  the  valence  force  field 
(VFF)^^  model  alone.  Thus,  all  the  factors  influencing  bond-length 
relaxation,  e.g.  strains,  boundary  conditions,  and  chemical  effects,  can  be 
studied.  The  ability  to  incorporate  the  chemical  effects  is  one  major 
difference  between  this  theory  and  other  VFF  models. 


The  rest  of  the  paper  contains  the  following  sections:  Section  II 
describes  a  theory  for  calculating  impurity  substitution  energies.  Section 


III  casts  the  theory  into  perturbation  form  and  combines  it  with  a  valence 
force  field  and  an  elastic  continuum.  Several  VFF  models  are  derived  in 
Section  IV.  The  modifications  of  numerical  results  due  to  chemical  effects 
on  local  bond  lengths  and  alloy  mixing  enthalpies  are  summarized  and 
discussed  in  Section  V,  Conclusions  are  drawn  in  the  last  section,  Sec¬ 
tion  VI. 


II.  IMPURITY  SUBSTITUTION  ENERGY 

Consider  the  problem  of  substituting  an  isoelectronic  atom  A  for  a  B 
atom  in  a  zinc-blende  compound  BC  (e.g.,  In  substitutes  for  Ga  in  GaAs,  as 
shown  in  Figure  1).  In  general,  the  bond  lengths  d^ ,  d^,  •*' 

first-,  second-,  and  third-shell  bonds  surrounding  the  impurity  are 
different  from  the  equilibrium  values  of  either  the  pure  BC  compound, 
denoted  as  d,  or  the  "impurity"  compound  AC,  denoted  as  d^  =  d  (1-6^).  If  A 
starts  being  a  free  atom,  and  B  also  ends  being  a  free  atom,  then  the  energy 
difference  between  the  final  and  initial  states  is  defined  as  the 
substitution  energy  and  is  given  by 
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A  sketch  of  the  flatted  picture  of  a  single  impurity  A 
in  a  BC  compound.  The  positions  and  displacements  for 
those  atoms  labeled  are  used  in  Appendix  B. 


is  the 


where  and  Cg  are  free-atom  energies  for  A  and  B  respectively, 
total  energy  of  the  semiconductor  with  a  defect  as  sketched  in  Figure  1,  and 
Epure  that  of  the  pure  host  BC  crystal.  Equation  (1)  can  be  written  as 


^s  "  ^^def"^dl3^~^^A''^B^*^^dis”^pure^ 


^R'^^A"^B^*^dis’ 


where  we  have  added  and  substracted  a  term  which  is  the  total  energy 

of  a  distorted  BC  compound  with  all  the  atoms  held  at  the  positions 
specified  in  Figure  1,  except  the  central  atom  is  a  B  atom.  Clearly,  » 
is  a  replacement  energy,  and  the  distortion  energy  A^^^  » 
Edig“Epure  energy  required  to  deform  a  pure  BC  crystal  from  its 

equilibrium  lattice  to  that  specified  in  Figure  1.  Ag  contains  all  the 
chemical  terms  that  arise  from  different  bond  lengths  and  polarities  between 
AC  and  BC  bonds. 

7 

Aq  can  be  treated  most  easily  by  Harrison’s  bonding  theory.  In  this 

K 

theory,  the  energy  per  bond  relative  to  the  vacuum  state  is 


E.  »  2g. +V  +6g  +6e  , 
b  b  o  m  m' 


where  is  the  energy  of  the  bonding  state  constructed  from  the  two  hybrid 
orbitals  facing  each  other  along  the  bond  direction 


(4) 


A  C 

with  and  e^  being  the  energies  of  the  anion  and  cation  hybrids  orbitals 

respectively.  The  antibonding  energy  e  has  the  same  form  as  in  Eq.(4) 

except  with  a  plus  sign.  is  called  the  covalent  energy,  which  is  the 

total  electronic  Hamiltonian  matrix  element  between  the  two  hybrids  in 

C  A 

question,  and  the  polar  energy  is  the  difference  =  1/2  (e^^  -e^  ).  The 
♦  “ 

E  and  E  ,  the  metallization  energies,  are  the  shifts  of  the  bonding 
m  m 

level  caused  by  interactions  with  the  neighboring  anti-bonding  states,  where 
and  -  indicate  whether  the  common  adjacent  atom  is  a  cation  or  anion.  For 
example,  this  term  for  an  AC  bond  labeled  by  d^  in  Fig.  1  due  to  an 
antibonding  state  labeled  by  d^  is  given  by 

e®'J(2,1)  -  luJ(1)l^lU®(2)l^  |V^|^/[eJ(1)-E®(2)], 


c  c  c 

where  A  and  B  denote  AC  and  BC  bonds  respectively,  V  -  l/i*  (e„  -e  )  with 

1  s  p 

c  c 

Eg  and  Ep  being  the  s-  and  p-  term  values  of  the  common  adjacent  atom  C. 

Up  (1)  is  the  probability  amplitude  of  finding  an  electron  in  the  hybrid 

orbital  of  the  C  atom  in  the  bonding  state  of  an  AC  bond  with  a  bond  length 

d  ,  whereas  U  (2)  is  the  corresponding  probability  amplitude  for  the 
1  a 


anti  bonding  state  of  a  BC  bond  of  bond  length  d^.  Finally  is  a  repulsive 
pair  potential  required  to  prevent  the  crystal  from  collapsing  and  to 


(5) 


guarantee  a  correct  equilibrium  bond  length. 


The  local  perturbation,  Eq.(5),  is  applicable,  because  the  square  root 
of  the  numerator  is  much  smaller  than  the  separation  between  the  anti¬ 
bonding  and  bonding  levels,  and  the  valence  band  is  completely  filled,  so 
the  interaction  between  the  bonding  states  only  spreads  the  levels  into 
bands  without  affecting  the  center  of  gravity  of  the  occupied  states. 

The  replacement  energy  Aj^  of  Eq.(2)  can  now  be  written  explicitly; 

The  distortion  energy  of  Eq.(2)  now  involves  only  BC  bonds  of  different 

bond  lengths.  It  can  be  treated  with  exactly  the  same  procedure  for  any 

given  set  of  bond-length  distributions.  Thus,  a  straightforward  energy 

minimization  procedure  can  be  carried  out.  The  accuracy  of  this  procedure, 

however,  depends  in  turn  on  the  accuracy  of  scaling  rules  for  and  and 
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the  input  parameters.  At  present,  Harrison’s  model  with  1/d  and  <* 

1) 

1 /d  and  his  universal  parameters  are  only  semi-quantitative.  We  are 
improving  the  quantitative  nature  of  the  theory  so  that  the  full  theory  will 
yield  accurate  predictions  of  the  structural  and  thermal  properties  of 
semiconductor  defects  and  alloys. 


III.  PERTURBATION  EXPANSION,  VALENCE  FORCE  FIELD,  AND  ELASTIC  MEDIUM 


As  pointed  out  earlier,  a  perturbation  expansion  of  the  theory  is 
instructive.  This  is  feasible  because  the  differences  AV2  -  V^(.AC)-\l^(BC) 

and  ■  V^CAO-V^CBC)  are  small  compared  to  each  individual  value  for  many 

of  the  Isoelectronic  impurities  in  III-IV  and  II-VI  compounds.  To  this  end, 
Eq.(6)  is  rewritten  as 

-6[e*‘®{2,1)-e®*®(1,1)]+6[e;J'®(1,2)-E®’®(2,2)] 
m,C  m,C  m,C  m,C 

-6Ce®'®(1.2)-e®'®(2,2)], 
in,c  ni,c 


A  B 

where  E^^  ( 1 )  and  E^^  ( 1 )  are  energies  per  bond  in  Eq.(3)  for  AC  and  BC 
compounds  respectively,  with  the  relaxed  bond  length  d^  »  d  (1-6).  The 
difference  between  these  energies  E^^  and  the  corresponding  values  at  their 
respective  equilibrium  bond  lengths  dj  and  d  are  just  the  strain  energies 
per  bond  in  uniform  deformation: 


eJ(1)  -  Ej(dj)+2/3  Bjdj(d^-dj)^, 


E®(1)  -  E®(d)+2  /3  Bd(d.-d)^, 


(8) 


where  and  B  are  the  bulk  moduli  for  the  "impurity"  AC  and  host  BC 
crystals.  The  rest  of  the  terms  in  Eq.(7)  are  all  due  to  changes  in 
caused  by  the  differences  AV^  and  AV^.  We  shall  use  Harrison's  scaling 

»T 

rules  to  deduce  them.  Expanding  Eq.(7)  to  second  order  in  AV^  and  Ad  = 
d^'d^ ,  we  write 


6Ce^J(2.1)-ef*J(1,1)]  =  f.Ad-g.  AV  +h  (Ad)^ 
m,C  m,L  1  1  d  1 


-WjAdAV^+UjCAV^)'^. 


where  f^,  g^,  and  so  on,  are  appropriate  derivatives  evaluated  for  the 
"impurity"  crystal  AC.  When  similar  expansions  are  made  for  the  rest  of  the 
terms  in  Eq.(7),  it  becomes  [with  d^  =  d(1-6),  d^  =  d(1-6^)] 

(1/4)Ap  -  AEj^+(fj-f)Ad-(g^-g)AV^+(hj+h)(Ad)^ 

-(Wj+W)  AdAV2+(Uj+U)(AV^)^ 

Bd^6^, 


where 


AE,  - 


E;(d,)-Ej(d)-(1/4)(e,-eJ 


Is  just  the  difference  in  the  binding  energy  per  bond  between  the  "BC"  and 

"AC"  crystals.  In  Eq.(IO),  the  coefficients  f,  g,  without  a  subscript  are 

those  for  the  host  BC  system.  It  is  convenient  to  define  an  excess  energy 

AE  -  A  /M-aE  ,  which  is  the  extra  energy  per  bond  required  for  the  impurity 
S  D 

substitution  over  and  above  the  binding  energy  difference  between  the  BC  and 
AC  crystals.  The  binding  energy  difference  accounts  for  much  of  the 
substitution  energy;  however,  the  correction  measured  by  the  excess  energy 
can  be  significant.  The  excess  energy  results  from  strain  energies  and 
chemically  driven  charge  redistributions  around  the  defect.  Using  Eqs.(2) 
and  (10)  and  defining  F  =  fj-'f  and  G  =  gj-g,  we  can  write  AE  up  to  second 
order  in  AV^and  Ad  as 


AE  -  2/3  Bjd^{6-6^)^-2/I  Bd^A^+FAd-GAV^ 
+H(Ad)^-WAdAV2+U(AV2)^+(1/H)A^^g. 


(12) 


where  H-h^+h,  W=Wj+w,  and  U  =  u^  +  u. 

To  treat  the  distortion  energy  we  divide  the  crystal  into  two 

regions.  Inside  a  sphere  of  some  radius  R  measured  from  the  impurity,  the 
strain  energy  is  taken  to  be  the  valence  force  f ield^*^value: 


A 


(in) 

dis 


2  / 

S  a,[A(d  *3  )]  2  2]  B.  v[A(3^*'3.)]^, 

Sd"^  i  ^  ^  ^  Sd'^  i  >  j 


(13) 
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where  i  sums  over  all  the  bonds  Inside  R,  and  the  pairs  in  the  B"terms 

include  those  that  have  adjacent  atoms  inside  R  and  on  the  boundary.  The 

parameters  a  and  B  are  force  constants  to  be  considered  later.  A(d*dj)  = 

d,  d  .-3,  *3 .  measures  the  change  of  the  dot  product  between  bond  vectors 

i  J  i  j 

due  to  distortions.  Outside  R,  we  assume  an  elastic  continuum  with  radial 
displacements  Inversely  proportional  to  the  square  of  the  distance  from  the 
center.  The  elastic  energy  in  this  medium  can  be  shown  to  be  (see 
Appendix  A) 


A 


(out) 

dis 


2 

RCu 


(14) 


where  the  effective  shear  coefficient  is  C  »  4ir[0.4(C^  ^ 

is  the  displacement  at  R.  In  view  of  the  fact  that  the  bonds  d^  and  d^  are 
coupled  through  the  chemical  terms  in  Eq.(IO),  the  smallest  logical  radius  R 
is  the  second-shell  atomic  distance,  namely  R  =  2/2d//J.  Atoms  on  this 

boundary  have  displacements  of  the  forms  u  =  d(Y,Y,o)//3,  ...  Thus,  u  =  /2 

y  d//3*  and  the  elastic  energy  in  the  continuum  is 


(out)  4  (15) 

dis  3  \  3 


The  distortion  energy  represented  by  Eq.(14)  contains  six  different 

contributions  (see  Appendix  B):  the  bond-stretching  energy  of  the  four 

2  2  2  2 
first-shell  bonds,  6a6  d  ,  the  B-terms  from  the  first-shell  bonds,  B6  d  , 


11 


the  a-terms  from  the  second-shell  bonds,  2a.(.5*2y)  d  ,  the  8“terms  between 

2  2 

the  first-and  second-shell  bonds,  2B(6+Y)  d  ,  the  6~terms  among  the  second- 
2  2 

shell  bonds  2B6  d  ,  and  finally  the  B”terms  between  the  second-shell  bonds 

2  2 

and  those  in  the  continuum,  (l/2)6d  C(36  +  X^ Y)  +  (6  +  X2^^  ^1  “ 

110  /2/(19  /il)  and  -  2-iiO  /2/(11  /TT ) . 

To  assemble  all  the  contributions  to  Eq.(12),  we  need  to  consider  the 
assignments  of  the  elastic  constants  and  the  force  constants  a  and  6  in 
VFF.  While  the  experimental  values^  ^  of  *^®  used  for 

the  elastic  constants,  a  and  6  have  to  be  deduced.  If  Martin's  original 

1 0  3 

procedure  (also  followed  by  Martins  and  Zunger  )  is  used,  then  Eq.(13) 

alone  will  not  produce  the  correct  (experimental)  bulk  moduli.  There  are 

small  corrections  due  to  Madelung  terms,  which  are  hard  to  treat  in  the  case 

of  non-uniform  distortions.  A  simpler  procedui^e  is  adopted  here.  We  use 

the  experimental  bulk  moduli  for  Bj  and  B  in  Eq.(12)  and  experimental 

elastic  constants  to  calculate  C  of  Eq.(l^),  and  then  force  a  and  B  in  the 

VFF  to  produce  the  correct  bulk  moduli  B  and  shear  coefficients 

7 

Such  an  approach  is  also  consistent  with  Harrison's  bonding  theory  and 
other  approaches  in  which  the  Coulomb  forces  are  automatically  incorporated 
in  the  band  and  bond  energies,  and  do  not  need  to  be  redundantly  treated. 
With  our  procedure,  the  bulk  modulus  is  simply  related  to  the  force 
constants  by  B  »  (3a  +  B)/(^/Jd).  Table  1  lists  our  a  and  B  values.  We 
want  to  point  out  in  advance,  however,  that  the  numerical  results  deduced 
from  our  and  Martin's  sets  of  a  and  B  do  not  introduce  differences  more 
than  the  present  experimental  uncertainties  in  the  local  bond  length 
(-  0.01A)  and  the  mixing  enthalpies  (>  0.5  kcal/mole). 


Table  1.  The  bond  lengths  d  (in  A,  valence  force  constants  a  and  6  (in  N/m) 
shear  coefficients  C  of  the  continuum  (in  10^^  erg/cm^), 
melting  temperatures  Tj^CK),  and  Liedermann  ratios  X^,  for  the 
compounds  used  in  this  paper. 


Compound 

d 

ct 

6 

C 

O'  ♦ 

Xm 

AlP 

2.367 

44.323 

8.068 

122.396 

1773 

0.070 

AlAs 

2.451 

40.849 

8.717 

112.695 

1873 

0.073 

AlSb 

2.656 

34.073 

6.900 

85.351 

1323 

0.062 

GaP 

2.360 

44.764 

10.737 

145.921 

1510 

0.064 

GaAs 

2.448 

39.235 

9.159 

121.844 

1738 

0.071 

GaSb 

2.640 

31.876 

7.347 

89.372 

985 

0.055 

InP 

2.541 

40.363 

6.543 

91.785 

1343 

0.059 

InAs 

2.622 

33.203 

5.752 

78.816 

1215 

0.061 

InSb 

2.805 

28.557 

4.891 

60.721 

798 

0.049 

ZnS 

2.342 

40.429 

5.273 

89.272 

2123 

0.081 

ZnSe 

2.454 

32.200 

4.562 

82.687 

1788 

0.080 

ZnTe 

2.637 

29.445 

4.659 

62.430 

1511 

0.071 

CdTe 

2.806 

26.569 

2.722 

38.453 

1371 

0.067 

HgTe 

2.798 

26.396 

2.746 

40.363 

943 

0.056 

Ref.  30 


Using  the  above  procedure  and  adding  all  contributions,  the  excess 
energy  per  bond  from  Eqs.(lO)  and  (12)  is  the  full  perturbation  theory  (FPT) 
result 


AE  =  C3aT(6-5  )^/2+6t(5-6„)^/2+c»(6  +  2y)^/2 
I  o  10 


+  6(6+Y)^/2+66^/'»+6(36  +  X^Y)^/8  +  B(6  +  X2Y)^/'l]d^ 


+  72  CY^d^/(3/3)+AE^^, 


where  the  chemical  contribution  is  written  as 


-  F  .  Ad+H(Ad)‘'+AE^ 
ch  ch  p 


I  F^^(26+Y)d+|  H(26+Y)^d^+(UAV2  -GAV^) , 


(17) 


where  AF  .  =  F-WAV_  and  AE  =  UAV^-GAV..  F  .  is  a  chemical  force,  which, 
ch  3  P  3  3  ch 

when  it  is  positive,  tends  to  push  the  C  atom  away  from  the  impurity  atom  A. 

This  force  arises  from  the  difference  in  the  bond  tensions  induced  between 

the  AC  and  BC  bonds  adjacent  to  C  because  the  neighboring  anti-binding 

states  are  different  from  those  of  their  respective  host  states.  AE  is  due 

P 

to  the  difference  in  the  polarities  AV^  alone  and  is  independent  of  the 
displacement.  Finally  H  can  be  regarded  as  a  chemically  induced  force 


-  - ' '  -  - 


V-V- 'V- 
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constant  which,  when  It  is  positive,  tends  to  restrain  the  lattice  from 
distortion  and  increases  the  elastic  energy. 

The  equilibrium  requirements,  3(AE)/36  =  0  and  3(AE)/3Y  =  0,  then  lead 
to  the  solution  Y  -  Q6,  and  6  is  given  by 

6  -  (6  +  6')  /{1->-[a(1-2Q)+B(17/^-\Q)->-l6H(1-2Q)/9]/(3aj*Bj)}. 

O  i 

where  the  constant  X  is  1 +3X^ /*<+X2/2,  and 

with  Q  -  2J/K,  J  =  a+XB/2+8H/9,  and  K 

-  Ha+2/2Cd/(3/5)+(1+X^/4+X2/2)e+8H/9. 

IV.  VALENCE  FORCE  FIELD  MODELS 

In  this  section  we  consider  several  models  based  on  the  valence  force 
field  (VFF).  These  models  have  been  used  frequently  to  explain  the  impurity 


3  14  8  9 

bond  relaxation  .  We  shall  first  derive  the  explicit  expressions  for 

these  models  and  then  connect  them  with  the  existing  results. 


Model  A :  Third-Shell  Atoms  and  Beyond  Are  Fixed  at  Their  Pure  Crystal 
Positions 


Let  the  bond  lengths  surrounding  the  impurity  again  be  d^  =  d(1-6) 

and  let  the  second-shell  atoms  have  radial  displacements  of  the  forms  (d//3) 

(Y.T.o)  etc.  Beyond  and  including  the  third-shell,  all  the  other  atoms  are 

held  at  their  pure-crystal  positions.  There  are  nine  different  contribu- 

butions  to  the  strain  energy  in  VFF  (see  Appendix  B):  the  a-terms  from  the 

four  bonds  surrounding  the  impurity:  6aj(6-6^)  d  ,  the  B'terms  among  the  six 

2  2 

pairs  of  these  bonds,  6j(5"6q)  d  ,  the  g-terms  between  the  first-shell  and 

2  2 

second-shell  bonds,  2B(6+Y)  d  ,  the  a-terms  from  the  second-shell  bonds, 

2  2  2  2 
2a(5+2Y)  d  ,  the  6"terms  among  the  second-shell  bonds,  2  g  6  d  ,  the  g-terms 

2  2  2  2 

between  the  second-shell  and  third-shell  bonds:  (9/2)g6  d  +  g(6+2Y)  d  ,  the 

2  2 

a-terms  from  the  third-shell  bonds,  8aY  d  ,  the  g-terms  among  the  third- 
2  2 

shell  bonds,  4gY  d  ,  and  the  g-terms  between  the  third  and  fourth-shell 
2  2 

bonds  66Y  d  .  Thus,  the  excess  energy  (in  this  case  1/4  times  the  strain 
energy)  becomes 


2  1  O  M  1  0  2 

1  2  1  2  1  ? 
+2  a(6+2Y)  +  g6  +  ^B(6+2Y) 

2  S  2  2 

+2aY  4sy  ]d  . 


(20) 
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The  minimization  of  AE  w.r.t  6  and  Y  leads  Eq.(20)  to  Y  *  -6/4,  and 


“  6q/[ 1 +(a+176/2)/ (6aj+6j) ] . 


We  note  that  there  is  some  ambiguity  in  the  third  contribution  listed 
above  for  the  6“terras  between  the  impurity  and  host  bonds.  The  value  of  B 
could  be  chosen  as  one  of  these  combinations  B.Bj.  1/2  (B+Bj)f  /BBj  or  other 
proper  combinations.  Because  the  values  of  B  and  Bj  are  comparable  and  B 
values  are  much  smaller  than  a  (see  Table  1),  the  results  for  6  and  AE  are 
not  too  sensitive  to  the  choice.  There  is  also  some  ambiguity  in  the  values 
for  for  the  "undistorted"  crystal.  The  -d^/3  used  is  the 

simplest  choice.  A  different  choice  will  not  affect  the  results  for  6  at 
all,  but  will  make  AE  slightly  different.  In  fact.  Model  A  was  firstly  used 

3 

by  Martins  and  Zunger.  However,  their  expression  for  6  is  different  from 

Eq.(21),  because  they  made  different  choices  of  the  two  quantities  Just 

mentioned.  Nevertheless,  Section  IV  will  show  that  these  two  expressions 

yield  very  similar  results.  These  ambiguities  do  not  occur  in  the  full 

theory  in  Section  III,  where  the  impurity-host  interactions  are  taken  into 

account  naturally  by  the  replacement  energy  A„  [see  Eq.(IO)]. 

K 

Model  B;  Second-Shell  Atoms  Connect  to  a  Fixed  Boundary 

This  model  corresponds  to  Y  =  0  in  Model  A.  So  we  have 

AE  =  [|  aj(6-6Q)^^i6j{6-6^)^+p6^-i|  66^]d^ 

17 


g 

The  spring  model  recently  discussed  by  Shih  et  al.  corresponds  to 
Eq.(25)  with  o  -  Oj ,  so  6/6^  »  3/^. 


Model  D.  VFF  With  the  Continuum  Connected  to  tf.  >  Second-Shell  Atoms 


Model  D1 


In  this  case,  AE  only  contains  the  first  five  contributions  listed  for 
Case  A  plus  the  elastic  energy  in  the  continuum.  However,  the  B-terms 
between  the  second-  and  third-shell  bonds  are  modified  because  atoms  outside 
R  in  the  continuum  now  have  radial  displacements  proportional  to  the  inverse 
of  the  square  of  the  radius.  The  result  is 


AE  -  [•|a,(6-6  )^+-r6x(5"5  )^+-5a(6  +  2Y)^+^S(6+Y)^ 
2  1  o  M  1  0  2  2 


■g6(36  +  A  + 


(26) 


where  A,  and  A 2  ^re  the  same  as  the  constants  that  appear  in  Eq.(l6).  The 
corresponding  equilibrium  condition  can  be  shown  to  be 


6  -  6Q/{1+[a(1-2Q)  +  196/'»-B(l+3A^/4  +  A2/'*)Q]/(3aj  +  8j/2) , 


(27) 
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Q  =  (2a+B  +  3X,8/‘*+ApB/2)/('<a+2/2Cd/3v^+8  +  i,X^6  +  2^2^^  • 


Model  D2 


A  comparison  between  Eq.(26)  and  the  full  perturbation  (FPT)  theory 
Eq.(l6)  shows  two  major  differences.  First,  all  the  chemical  terms  are 
absent  in  Eq.(26).  Secondly,  the  terms  ( 1 /2)6j (6-6Q)^d^+( 1 /4) 6  6^d^  in 
Eq.(l6)  become  (1/4)  6j  (6  -6Q)^d^+(1/2)  B  6^d^  inEq.(26).  This  difference 
in  the  strain  energy  will  mask  the  true  effects  of  chemical  forces  if  S  from 
Eq.(26)  is  compared  with  FPT.  A  better  way  to  study  the  chemical  effects  is 
to  use  the  following  equation: 


2  1  2  X  0  2  2 


./2  .23 
" 373  CT  d  , 


which  is  Eq.(l6)  with  all  the  chemical  terms  neglected.  The  corresponding  6 
becomes 


6  -  6Q/{1+Ca(1-2q)+-J^B-S(l4^^4^2^®^''^^“l'’®I^^’ 


with  Q  still  given  by  Eq.(28). 


Model  E:  Continuum  Connected  to  the  First-Shell  Atoms 

In  this  case  Y  =  “3  /3  6/ (8/2)  and  AE  only  Includes  the  first  three 
contributions  listed  In  Model  A  plus  the  strain  energy  of  the  continuum: 


1  2  2 
+  -jjCdA  Dd'". 


The  relaxation  parameter  is  given  by 

6  »  Cd+(1-^^)^6]/(3aj+^6j)}  • 


We  note  that  the  continuum  model  used  to  estimate  the  bond-length  relaxation 

9 

by  Baldereschi  and  Hoffield  corresponds  to  Eq.(32)  without  the  6  terms, 


which  yields  6/5  -  to  0.5,  rather  than  the  proper  values  around  0.7  to 


V.  ALLOY  MIXING  ENTHALPY 


The  impurity  substitution  excess  energies  AE  provide  a  first  estimate 

of  the  mixing  enthalpies  of  pseudo-binary  alloys.  Most  current 

thermodynamics  theories  of  semiconductor  alloys  are  based  on  an  extension  of 

1 2 

the  binary  solution  model.  In  this  model,  the  mixing  Helmholtz  energy  of 
an  alloy  is  defined  as 


^^m  “  ^alloy"^^^AC*^^BC^’ 


where  y  =  1-x,  and  and  are  the  respective  free  energies  of  the  pure 
AC  and  BC  compounds  at  the  same  temperature.  Because  the  C  atoms  occupy  a 
sublattice,  the  nearest  neighbors  of  A  and  B  atoms  in  the  alloy  are  the  C 
atoms.  Thus  the  pair  potentials  that  enter  the  binary  solution  theory  are 
now  the  second-neighbor  interactions.  Let  N,,,  N,„,  and  N__  be  respectively 
the  numbers  of  the  second-neighbor  AA,  AB,  and  BB  pairs,  with  corresponding 


pair  interaction  energies,  e,,,  e,„,  and  Cnn •  For  tetrahedral 

AA  AB  BB 

semiconductors,  there  are  a  total  of  6N  second-neighbor  pairs  for  a  crystal 
containing  N  unit  cells.  Denote  the  ratios  N^g,  and  Ngg  to  6N  as  r^^, 

r,_  -  r,  and  r^-  respectively.  Then  those  ratios  are  related  to  the  alloy 
composition  by  r^^  «  x-r/2  and  rgg  =  y-r/2.  The  mixing  free  energy  has  two 


terms 


AF  = 


AE  -TAS, 


where  the  mixing  energy  is  given  by 


'aa  V'ab''ab*'bb'’bb  «AA*»=BB> 


6NrA€ , 


where 


Ac 


The  mixing  entropy  AS  can  also  be  written  from  a  simple  generalization  of 

1 2 

the  random  distribution.  For  modest  pressure,  AE  is  the  same  as  the 
mixing  enthalpy  AH^. 


Now  the  pair  interaction  energies  can  be  approximately  related  to  the 
impurity  substitution  energies  by 


A^(A  1„  BC)  . 


(37) 


and 


Ag(A  in  BC)  »  12 


(38) 


Thus  Ae  of  Eq.(36)  becomes 


Ae  -  4,, [A  (A  in  BC)+A„(B  in  AC)] 
2^  s  s 

«  ^  [AE{A  in  BC)+AE(B  in  AC)]. 


Usually,  the  experimental  AH  is  written  as 

m 


AH^  »  x(1-x)n, 


which  is  equivalent  to  assuming  a  random  distribution,  i.e.  r  =  2x(1-x)  . 
Using  this  expression  for  r  and  comparing  Eq.(i<0)  and  (35),  we  see  that  the 
mixing  enthalpy  parameter  fi  is  given  by 


n  -  2[AE(A  in  BC)  +  AE(B  in  AC)]. 


This  connection  provides  a  further  check  of  the  theory. 


VI.  NUMERICAL  RESULTS  AND  DISCUSSION 


A.  Chemical  Terms 


Table  2  lists  6^  =  l-d^/d,  6  =  1-d^/d,  the  excess  energy  (per 
bond)  AE  for  the  full  theory  and  its  corresponding  VFF  Model  D2,  anc  the 
terms  derived  from  the  metallization  energies:  6'  (Eq.(19),  F  ,  H,  AE  , 

O  Un  p 

and  AEj,^  (Eq.(17)].  The  appropriate  derivatives  f,  g,  h,  [see 

Eq.(9)]  are  computed  using  the  atomic  term  values  that  we  have  generated 

1  ^  5 

form  impurity-level  and  structural  studies. 

For  substitutions  involving  the  cation  pair  (Ga,  Al),  F^^  has  the  same 

sign  as  6  ,  which  means  that  F-.  prevents  relaxation  and  thus  tends  to 
o  w  n 

increase  the  strain  energy.  The  chemical  forces  H  are  also  significant.  As 

a  result,  all  six  cases  involving  this  pair  have  nearly  equal  d^  and  d^, 

i.e.  the  small  bond  length  differences  are  made  even  smaller  in  the  alloy. 

The  excess  energies  all  become  negative,  mainly  because  AE^  is  negative. 

For  the  systems  Involving  the  (Ga,  In)  and  (In,  Al)  pairs,  F^^^  has  the 

opposite  sign  from  6^,  so  6^  and  6^  have  the  same  sign.  The  chemical  force 

favors  bond  distortion.  However,  because  H  is  positive  and  it  introduces  an 

increase  in  the  denominator  of  Eq.(l8),  most  of  the  effect  of  6^  is 

cancelled.  For  cases  Involving  (Ga,  In),  the  polarity  contributions  AE^  are 

2 

all  negative.  The  fj^^Ad  term  is  negative,  but  HAd  is  positive,  so  they 

cancel  to  a  certain  degree  and  leave  AE  lowered  primarily  because  of  AEp. 

While,  AEp  is  still  negative  for  the  (In,Al)  substitutions,  its  magnitude  is 

2 

reduced  considerably.  The  other  chemical  energies  f^^Ad  +  HAd  can  be  as 


Table  2.  Comparison  between  the  full  theory  and  the  corresponding  VFF  Model  D2 
to  study  the  effects  of  chemical  terms.  All  AE's  are  in  units  of 
kcal/mol-band . 


Full  Theor\ 


Host 

Model  D2 

6 

AE 

1 

AlP 

0.003 

0.002 

0.001 

1 

GaP 

-0.003 

-0.002 

0.001 

.a 

AlAs 

0.001 

0.001 

0.000 

A1 

GaAs 

-0.001 

-0.001 

0.000 

Ga 

AlSb 

0.006 

0.004 

0.005 

A1 

GaSb 

-0.006 

-0.004 

0.005 

In 

GAP 

-0.077 

-0.052 

0.959 

Ga 

InP 

0.071 

0,056 

0.734 

In 

GaAs 

-0.071 

-0.048 

0.752 

Ga 

InAs 

0.066 

0.052 

0.592 

In 

GaSb 

-0.062 

-0.043 

0.554 

Ga 

InSb 

0.059 

0.046 

0,445 

In 

AlP 

-0.074 

-0.053 

0,761 

A1 

InP 

0.068 

0.053 

0.674 

In 

AlAs 

-0.070 

-0.048 

0.705 

A1 

InAs 

0.065 

0.052 

0.576 

In 

AlSb 

-0.056 

-0.039 

0.440 

A1 

InSb 

0.053 

0.042 

0.368 

Cd 

ZnTe 

-0.064 

-0.048 

0.432 

Zn 

CdTe 

0.060 

0.050 

0.314 

Hg 

CdTe 

0.003 

0.002 

0.001 

Cd 

HgTe 

-0.003 

-0.002 

0.001 

Hg 

ZnTe 

-0.061 

-0.045 

0.392 

Zn 

HgTe 

0.058 

0.048 

0.286 

As 

AlP 

-0.03 

-0.026 

0.179 

P 

AlAs 

0.034 

0.025 

0.185 

As 

GAP 

-0.037 

-0.025 

0.226 

P 

GaAs 

0.036 

0.027 

0.211 

As 

InP 

-0.032 

-0.023 

0.136 

P 

InAs 

0.031 

0.024 

0,128 

Sb 

AlAs 

-0.084 

-0.058 

1.024 

As 

AlSb 

0.077 

0.059 

0.919 

Sb 

GaAs 

-0.078 

-0.052 

0.908 

As 

GaSb 

0.073 

0.055 

0,823 

Sb 

InAs 

-0.070 

-0.051 

0.603 

As 

InSs 

0.065 

0.051 

0.551 

Sb 

AlP 

-0.122 

-0.085 

2.007 

P 

AlSb 

0.109 

0.085 

1.855 

Sb 

GAP 

-0.119 

-0.075 

2.132 

P 

GaSb 

0.106 

0.083 

1,806 

Sb 

InP 

-0.104 

-0.072 

1.383 

P 

InSb 

0.094 

0.077 

1.191 

Se 

ZnS 

-0.048 

-0.036 

0.231 

S 

ZnSe 

0.046 

0.037 

0.221 

Te 

ZnSe 

-0,075 

-0.056 

0.550 

Se 

ZnTe 

0.069 

0.054 

0.532 

Te 

ZnS 

-0,126 

-0.092 

1.565 

S 

ZnTe 

0.112 

0.091 

1.446 

ohdO'-^^N) 

H(N/m) 

AEp 

AE.^ 

-  -  - — -  - - 

0.150 

4.581 

-0.016 

-0.016 

-0.150 

4.581 

-0.016 

-0.016 

0.243 

5.733 

-0.020 

-0.021 

-0.243 

5.733 

-0.020 

-0.021 

0.389 

5.632 

-0.054 

-0.053 

-0.389 

5.632 

-0,054 

-0.053 

0.699 

3.778 

-0.188 

-0.219 

-0,699 

3.778 

-0.188 

-0.206 

0.804 

4.778 

-0.257 

-0.283 

-0.804 

4.778 

-0.257 

-0.265 

0.352 

5.201 

-0.363 

-0.308 

-0.352 

5.201 

-0.363 

-0.287 

0.769 

3.506 

-0.035 

-0.087 

-0.769 

3.506 

-0.035 

-0.083 

0.942 

4.437 

-0,048 

-0.111 

-0.942 

4.437 

-0.048 

-0.099 

0.689 

4.979 

-0.061 

-0.073 

-0.689 

4.979 

-0.061 

-0.061 

0.202 

-0.484 

-0.005 

-0.064 

-0.202 

-0.484 

-0.005 

-0.072 

-0.278 

-0.753 

-0.018 

-0,026 

0.278 

-0.753 

-0.018 

-0.026 

0.075 

0.002 

0.052 

0.037 

-0.075 

0.002 

0.052 

0.035 

-0.085 

0.717 

-0.005 

0.008 

0.085 

0.717 

-0.005 

0.008 

-0.181 

1.078 

-0.011 

0.012 

0.181 

1.078 

-0.011 

0.014 

-0.057 

0.919 

-0.003 

0.008 

0.057 

0.919 

-0.003 

0.009 

-0.815 

0.644 

-0.180 

0.002 

0.815 

0,644 

-0.180 

0.027 

-1.599 

0.927 

-0,363 

-0.106 

1.599 

0.927 

-0.363 

-0.061 

-0,824 

0.855 

0.171 

-0.009 

0.824 

0.855 

-0.171 

0.008 

-0.944 

0.645 

0.241 

0.074 

0.944 

0.645 

-0.241 

0.123 

-1.868 

0.930 

-0.505 

-0.046 

1.868 

0.930 

-0.505 

0.093 

-0.922 

0.854 

-0.214 

0.059 

0.922 

0,854 

-0.214 

0.123 

0.077 

0.645 

0,003 

0  .000 

-0.077 

0.645 

0.003 

0.001 

0.028 

0.635 

0.008 

0.024 

-0.028 

0.635 

0.008 

0.025 

0.101 

0.644 

T  022 

0.041 

-0.101 

0.644 

0.022 

0.051 

-0,013 

-0.013 

-0.018 

-0.018 

-0.039 

-0.039 

0.742 

0.530 

0.472 

0.330 

0.247 

0,160 

0.679 

0.596 

0.602 

0.485 

0.369 

0.310 

0.373 

0.247 

-0.018 

-0.018 

0.429 

0.322 

0.187 

0.194 

0.240 

0.228 

0.144 

0.138 

1.060 

0.984 

0.929 

0.904 

0.645 

0.613 

2.127 

2.030 

2.244 

2.084 

1.501 

1.379 

0.231 

0.222 

0.574 

0.557 

1.606 

1.496 
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large  as  AE^.  But  the  overall  reductions  of  AE  are  only  a  fraction  of  those 


for  the  (Ga,  In)  cases.  For  the  several  II-VI  systems  studied,  both  and 


H  are  small,  and  the  net  changes  in  6  have  the  same  sign  as  6^.  However, 


because  6  is  small  in  the  (Cd,  Hg)  substitutions,  F„  actually  causes  a 
o  u  n 


reversal  in  which  the  short  bond  length  gets  shorter  and  the  longer  one  gets 
longer.  This  is  the  only  exceptional  case  of  this  type  found  for  all  the 
systems  studied.  The  change  of  AE  due  to  chemical  terms  in  the  (Hg,  Zn) 


substitution  is  also  peculiar — it  increases  mainly  because  AE^  is  positive. 


Next  examine  the  anion  substitutions.  For  the  groups  involving  the  (P, 
As)  pair,  the  chemical  shifts  are  all  small,  but  the  trend  is  less  toward 
relaxation  and  larger  AE.  The  groups  involving  (As,  Sb)  and  (P,  Sb)  pairs 


behave  very  similarly:  F^^  are  significant  and  are  opposing  relaxations. 


i.e.  <Sq  and  <S^  have  opposite  signs.  At  the  same  time,  the  H  values  are 


several  times  smaller  than  those  for  the  corresponding  III-V  cation 


substitution  case.  Thus,  most  of  6'  translates  into  a  real  reduction  of 

o 


the  ratio  6/6^,  and  consequently  introduces  extra  strain  energy.  Although 


the  AE  energies  are  significant  and  negative,  F„  Ad  are  positive,  and  the 
p  V  n 


net  AE^^  can  be  either  positive  or  negative.  However,  the  induced  strain 


energy  due  to  reduction  of  the  6/6^  makes  all  AE  positive  for  these  two 


groups  of  systems.  For  II-VI  systems,  all  the  chemical  effects  again  are 
small,  but  the  net  chemical  changes  on  AE  are  slightly  positive. 


The  above  discussion  can  be  summarized  in  Figure  2,  where  the  excess 
energies  AE  calculated  from  the  full  perturbation  theory  (FPT)  and  Model  D2 

Q 

are  plotted  normalized  to  the  results  of  the  simple  spring  model  (SSHS)°,of 
Shih  et  al .  i.e.  Eq.(2ii)  with  aj=  a  ,  so  AE/(3/8  ad^)  =  The  calculated 


.2 


lOx  6q  »  (dim^ntionlMv) 


The  excess  energies  AE  over  3/8  ad^  calculated  from  the 
full  perturbation  theory  and  its  corresponding  VFF 
Model  D2.  The  solid  curve  corresponds  to  Eq.  (24) 


with  a  “ 


1 .  In  fact 


AE  rises  faster  for  6^  ^  o  than  for  6^2  o,  mainly  because  a/a^  - 

if  the  relation^’ -  (dj/d)^  with  S  of  order  of  3  to  5  is  used  in 

Eq.(24),  we  obtain  a  percentage  correction  of  56^/^  to  the  SSHS  results, 

which  explains  the  skewed  behavior  of  the  curve.  It  is  also  clear  that  AE 

2 

rises  faster  than  6  for  larger  6  .  However,  the  zeroth  order  theory  of 
0  0 

SSHS  is  clearly  an  excellent  representation  of  the  global  features  of  AE. 

The  results  from  Model  D2  are  closer  to  the  parobolic  form  than  those  from 
FPT.  The  figure  clearly  shows  the  general  trends;  the  chemical  terms  cause 
negative  shifts  in  AE  for  cation  substitutions  and  positive  shifts  for  anion 
impurities.  It  is  also  clear  that  the  chemical  shifts  can  be  very  large. 
These  effects  will  have  important  consequences  on  the  alloy  mixing 
enthalpies  to  be  discussed  later. 


B.  Impurity  Bond  Length 

Table  3  lists  the  impurity  bond  lengths  d^  calculated  from 
different  models,  while  a  comparison  of  theory  and  the  available 

1  1  li  1  <5 

experimental  data  ’  is  presented  in  Table  1|.  The  actual  size  of 

changes  in  d^  Induced  by  the  chemical  terms  can  be  seen  by  comparing  Model 
D2  with  the  full  theory.  Except  for  the  systems  involving  the  substitution 
of  (Ga,As)  and  (P,Sb)  pairs  (where  changes  range  from  0.01  to  0.03A),  all 
the  chemically  induced  changes  are  less  than  0.01  A.  Comparison  among 
Models  A,  B,  and  C  shows  that,  while  extending  the  boundary  helps  the 
relaxation  (compare  Model  B  to  A),  i.e.  6/5^  is  closer  to  1,  the  inclusion 
of  the  bond-bending  forces  (the  6  terms)  (compare  Model  B  with  C)  prevents 
it.  The  simple  spring  model  (Model  C)  which  contains  neither  of  these 
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Table  3.  Calculated  Impurity  local  bond  lengths  (in  A)  from  the  full  theory 
and  several  valence  force  models  discussed  in  Section  III  and 
their  comparison  with  the  values  calculated  by  Martins  and 
Zunger  (Ref .  3) 


A 

B 

C 

D1 

D2 

E 

FPT 

MZ 

Ga 

AlP 

2.362 

2.363 

2.362 

2.362 

2.362 

2.364 

2.364 

A1 

GaP 

2.365 

2.364 

2,365 

2.365 

2.365 

2.365 

2.363 

Ga 

AlAs 

2.449 

2.449 

2.449 

2.449 

2.449 

2.450 

2.452 

Al 

GaAs 

2.450 

2.450 

2.450 

2.450 

2.450 

2.450 

2.447 

Ga 

AlSb 

2.645 

2.646 

2.644 

2.645 

2.644 

2.649 

2,653 

Al 

GaSb 

2.651 

2.650 

2.652 

2.652 

2.652 

2.650 

2.644 

In 

GAP 

2.477 

2.462 

2.492 

2.479 

2.483 

2,435 

2.487 

2.474 

Ga 

InP 

2.406 

2.421 

2.402 

2.395 

2.399 

2.409 

2.396 

2.409 

In 

GaAs 

2.559 

2.544 

2.573 

2.561 

2.565 

2.518 

2.570 

2.556 

Ga 

InAs 

2.492 

2.506 

2.486 

2.482 

2.485 

2.496 

2.481 

2.495 

In 

GaSb 

2.747 

2.734 

2.760 

2.749 

2.754 

2.710 

2.750 

2.739 

Ga 

InSb 

2.683 

2.697 

2.678 

2.673 

2.676 

2.686 

2.680 

2.683 

In 

AlP 

2.487 

2.472 

2.494 

2.490 

2.493 

2.447 

2.498 

2.480 

Al 

InP 

2.412 

2.427 

2.408 

2.401 

2.405 

2.415 

2.400 

2.414 

In 

AlAs 

2.561 

2.546 

2.572 

2.563 

2.568 

2.523 

2.575 

2.553 

Al 

InAs 

2.493 

2.506 

2.487 

2.483 

2.487 

2.497 

2.480 

2.495 

In 

AlSb 

2.754 

2.741 

2.763 

2.756 

2.760 

2.721 

2.765 

2.746 

Al 

InSb 

2.693 

2.705 

2.689 

2.685 

2.687 

2.696 

2.683 

2.693 

Cd 

ZnTe 

2.756 

2.740 

2.760 

2.760 

2.763 

2.720 

2.770 

2.755 

Zn 

CdTe 

2.673 

2.688 

2.676 

2.660 

2.665 

2.671 

2.658 

2.674 

Hg 

CdTe 

2.800 

2.801 

2.800 

2.800 

2.799 

2.801 

2.790 

Cd 

HgTe 

2.804 

2.803 

2.804 

2.805 

2.805 

2.804 

2.813 

Hg 

ZnTe 

2.750 

2.735 

2.754 

2.753 

2.757 

2.715 

2.758 

2.748 

Zn 

HgTe 

2.671 

2.685 

2.674 

2.659 

2.664 

2.671 

2.662 

2.673 

As 

AlP 

2.425 

2.418 

2.429 

2.427 

2.428 

2.406 

2.427 

2.422 

P 

AlAs 

2.392 

2.399 

2.387 

2.387 

2.389 

2.394 

2,390 

2.395 

As 

GAP 

2.417 

2,409 

2,424 

2.417 

2,420 

2.396 

2.416 

2,414 

P 

GaAs 

2.386 

2.393 

2.380 

2.382 

2.383 

2.389 

2.386 

2.387 

As 

InP 

2.596 

2.589 

2.599 

2.598 

2.600 

2.579 

2.598 

2.595 

P 

InAs 

2.561 

2.568 

2.558 

2.557 

2.558 

2.563 

2,560 

2.562 

Sb 

AlAs 

2.584 

2.566 

2.597 

2.587 

2.592 

2.539 

2.577 

2.574 

As 

AlSb 

2.506 

2.522 

2.496 

2.495 

2.498 

2.511 

2.514 

2.510 

Sb 

GaAs 

2.569 

2.553 

2.584 

2.571 

2.576 

2.524 

2.546 

2.564 

As 

GaSb 

2.501 

2.516 

2.489 

2,492 

2.495 

2.508 

2.525 

2.505 

Sb 

InAs 

2.747 

2.730 

2.754 

2.750 

2.754 

2,705 

2.733 

2.739 

As 

InSb 

2.669 

2.683 

2,663 

2.658 

2.662 

2.672 

2.683 

2.667 

Sb 

AlP 

2.555 

2.529 

2.569 

2.561 

2.568 

2.488 

2.550 

2.542 

P 

AlSb 

2.440 

2.462 

2,426 

2.425 

2.430 

2.447 

2,448 

2,444 

Sb 

GAP 

2.526 

2.503 

2,551 

2.529 

2.537 

2.461 

2.504 

2.519 

P 

GaSb 

2.431 

2.451 

2.414 

2.418 

2.422 

2.440 

2.454 

2.436 

Sb 

InP 

2.712 

2.687 

2.720 

2.719 

2.725 

2.654 

2.702 

2.700 

P 

InSb 

2.599 

2.619 

2,591 

2,585 

2.590 

2.604 

2.611 

2.597 

Se 

ZnS 

2.420 

2.409 

2.421 

2.424 

2.426 

2.396 

2.426 

2.420 

S 

ZnSe 

2.367 

2.376 

2.365 

2.360 

2.364 

2.370 

2.363 

2.367 

Te 

ZnSe 

2.586 

2.569 

2.588 

2.589 

2.592 

2.540 

2.591 

2.584 

Se 

ZnTe 

2.501 

2.517 

2.497 

2.490 

2.494 

2.504 

2.495 

2.502 

Te 

ZnS 

2.543 

2.513 

2.544 

2.552 

2.558 

2.478 

2.557 

2.539 

S 

ZnTe 

2.406 

2.429 

2.400 

2.390 

2.396 

2.410 

2.396 

2.407 
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terms,  evidently  represents  a  delicate  cancellation  of  these  effects  and 
predicts  results  close  to  those  of  the  full  perturbation  theory  and 
experiment.  Although  the  d^  values  of  Model  C  are  often  very  close  to  those 
of  Model  A,  there  are  cases  [e.g.,  Ga(P,  Sb)]  in  which  Model  C  can  differ 
from  Model  A  by  0.025  A.  Model  A  produces  about  the  same  d^  values  as  Model 
D1 ,  where  the  maximum  difference  in  d^  la  only  0.015  A.  Martins  and  Zunger^ 
used  the  same  model  as  Model  A,  however,  their  analytic  expression  for  6  is 
different  from  that  given  by  Eq.(21).  Nevertheless  numerical  results 
indicate  these  two  calculations  agree  to  0.01  A.  The  slightly  different 
forma  of  strain  energies  used  in  Model  D2  and  D1  only  introduce  a  small 
change  in  d^  with  the  largest  difference  being  less  than  0.01  A.  The  first- 
shell  continuum  model  (Model  E)  allows  too  little  relaxation,  so  while  the 
other  models  produce  a  ratio  5/6^  ranging  from  0.6  to  0.8,  Model  E  only 
ranges  from  0.^  to  0.6.  The  reason  that  the  fixed  boundary  in  Model  A  works 
is  that  the  effective  shear  coefficient  C  (see  Table  1)  characterizing  the 
strain  energy  in  the  elastic  continuum  is  large.  However,  Model  3  is  too 
rigid,  and  does  not  provide  enough  buffer  between  the  impurity  bond  and  the 
fixed  boundary. 

The  comparison  of  the  theoretical  results  with  the  available 

experimental  data  in  Table  indicates  that  Models  B  and  E  are  the  least 

accurate.  Models  A,  D1  ,  MZ  and  the  full  theory  are  comparable  in  that  all 

have  an  average  absolute  deviation  of  0.012  '  which  is  about  the 

experimental  uncertainty  in  EXAFS.  The  agreement  between  theory  and 

experiment,  however,  is  not  uniform.  The  most  surprising  result  in  Table 

is  that  the  simple  spring  model.  Model  C,  and  its  cruder  version  used  by 
8 

Shih  et  al .  (a  -  a,  so  6/6  «  0.75t  labeled  as  SSHS)  has  the  smallest 


variance  in  ,  about  0.006  A.  We  know  this  does  not  imply  that  the  simple 
spring  model  represents  the  real  picture  of  bond-length  relaxation.  Nature 
is  just  playing  tricks  on  us  again.  For  example,  if  we  let  all  the  shear 
coefficients  be  zero,  i.e.,  g-C-Oin  our  model,  then  as  the  range  of  the 
boundary  is  gradually  extended,  the  local  bond  length  will  eventually  relax 
to  the  impurity  bond  length  d^  -  dj ,  or  6  «  6^.  This  can  be  seen  in  Model  A 
from  Eq.  (21),  where  S  reduces  to  6Q/(1+a/6otj) ,  rather  than  6Q/(l+a/3aj)  as 
predicted  by  Model  C,  and  in  Model  D  from  Eqs.(27)  and  (30),  6  becomes  6^, 
if  the  continuum  is  taken  to  be  shearless.  Considering  that  various  effects 
are  included  that  may  mask  the  absolute  accuracy  of  d^  predictions  (e.g. 
while  low- temperature  bond  lengths  are  used  in  the  calculation,  the  room- 
temperature  values  of  elastic  constants  are  adopted),  the  agreement  of 
various  models  with  experiments  in  Table  should  be  regarded  as  excellent. 
There  are,  however,  many  other  impurity  systems  in  which  the  simple  spring 
model  predictions  differ  considerably  from  other  models,  as  is  shown  in 
Figure  3,  where  6  is  plotted  against  6^  for  the  full  theory.  Those  points 
that  deviate  significantly  from  the  0.75-slope  line  are  the  systems  with 
(As,  Sb)  or  (P,  Sb)  substitutions.  Additional  EXFAS  measurements  on  these 
systems  are  needed  to  test  these  predictions. 


C.  Mixing  Enthalpies 

Table  5  lists  the  mixing  enthalpy  parameters  n  (in  kcal/mole)  for 
a  number  of  alloys  estimated  from  Eq.(ill)  for  all  the  models  considered 
along  with  other  theoretical"^'  and  experimental  values.  As  already 


discussed,  the  chemical  terms  reduce  the  excess  energies  in  the  cation 


Table  5.  Mixing  enthalpy  parameter  ?.  (in  kcal/mole)  estimated  from  the  full 

perturbation  theory  (FPT)  and  several  valence  force  models  discussed 
in  Section  III  and  comparison  with  experiments  and  other  theories. 


A 

B 

C 

D1 

D2 

(Ga.ADP 

0.00 

0.01 

0.00 

0.01 

0.01 

(Ga  y  A1 ) As 

0.00 

0.00 

0.00 

0.00 

0.00 

(Ga,.Al)Sb 

0.02 

0.03 

0.02 

0.02 

0.02 

(Ga,In)P 

3.76 

4.79 

3.0. 

3.29 

3.39 

(Ga, In)As 

2.97 

3.76 

2.36 

2.61 

2.69 

(Ga,In)Sb 

2.22 

2.83 

1.77 

1.95 

2.00 

(In,Al)P 

3.24 

4.22 

2.77 

2.78 

2.87 

(In,.Al)As 

2.86 

3.65 

2.32 

2.49 

2.56 

(In.ADSb 

1.81 

2.33 

1.49 

1.57 

1.61 

(Cd,2n)Te 

1.80 

2.43 

1.73 

1.43 

1.49 

(Hg,Cd)Te 

0.00 

0.00 

0.00 

0.00 

0.00 

(Hg,2n)Te 

1.63 

2.20 

1.56 

1.30 

1.36 

Al(P,As) 

0.81 

1.03 

0.65 

0.71 

0.73 

Ga (P , As ) 

0.95 

1.18 

0.70 

0.86 

0.87 

In(P,.\s) 

0.60 

0.78 

0.52 

0.51 

0.53 

Al(As,Sb) 

4.31 

5.45 

3.38 

3.80 

3.88 

Ga(As,Sb) 

3.77 

4.69 

2.81 

3.40 

3.46 

In(As,Sb) 

2.61 

3.39 

2.23 

2.24 

2.31 

Al(P,Sb) 

8.60 

10.99 

6.99 

7.54 

7.73 

Ga(P,Sb) 

8.54 

10.61 

6.36 

7.72 

7.88 

In(P,Sb) 

5.87 

7.64 

5.08 

4.99 

5.15 

Zn(S,Se) 

1.04 

1.39 

0.98 

0.85 

0.90 

Zn ( S , Se ) 

1.04 

1.39 

0.98 

0.85 

0.90 

Zn(Se,Te) 

2.47 

3.27 

2.23 

2.09 

2.16 

Zn(S,Te) 

7.02 

9.34 

6.45 

5.80 

6.02 

E 

FPT 

MZ^ 

dl’’ 

Exp/ 

0.01 

-0.05 

0.00 

-0.07 

0.02 

0.02 

0.03 

0.11 

0.0 

0.03 

-0.15 

0.02 

0.02 

0.03 

0.0 

5.24 

2.54 

4.56 

3.63 

2.94 

3.25 

3.5 

4.14 

1.60 

2.49 

2.81 

2.42 

1.25 

1.65 

2.0 

3.0 

3.09 

0.81 

2.53 

1.85 

1.83 

1.47 

1.9 

4.60 

2.55 

3.93 

2.17 

3.60 

2.81 

2.37 

2.5 

2.50 

1.36 

2.06 

1.46 

1.45 

0.6 

2.45 

1.24 

2.12 

1.97 

1.63 

1.34 

0.00 

-0.07 

0.7 

1.4 

2.23 

1.50 

1.91 

1.81 

1.48 

3.0 

1.14 

0.76 

1.32 

0.94 

1.15 

0.98 

0.66 

0.12 

0.4 

1.0 

0.84 

0.57 

0.72 

0.58 

0.52 

0.4 

5.92 

4.09 

5.22 

3.67 

4.58 

3.35 

2.76 

4.0 

4.5 

3.67 

2.52 

2.89 

2.29 

2.17 

6.65 

2.25 

2.9 

12.00 

8.32 

11.66 

8.66 

8.04 

5.76 

1.49 

0.90 

1.49 

0.90 

3.63 

2.26 

2.91 

3.11 

2.12 

3.12 

1.55 

9.72 

6.20 

^Ref.  3.  Col.  A  of  Table  III. 
^Ref.  16 


impurities  and  increase  them  for  anion  impurities.  The  corresponding 
changes  in  P  are  the  differences  between  the  FPT  and  Model  D2.  We  note  that 
the  reductions  of  Q  for  the  (Ga,  In)  alloys  are  very  large  (>  1  kcal/mole) 
and  also  significant  for  (In,  Al)  alloys.  However,  the  increases  in  n  for 
the  anion  substitutional  alloys  are  not  as  large.  Also  the  p  differences 
between  Models  D1  and  D2  are  less  than  10  percent.  Model  A  produces  p 
values  about  20  percent  larger  than  Model  D1 ;  Model  B  in  turn  is  20  percent 
higher  than  Model  A;  and  Model  E  is  10  percent  higher  than  Model  B.  The  P 
values  in  the  simple  spring  model.  Model  C,  are  seen  to  be  about  the  same  as 
Model  01,  although  the  differences  among  systems  can  be  positive  or 
negative.  Although  MZ  used  the  same  strain  model  as  Model  A,  their  p  values 
do  not  agree  with  our  Model  A  values  because  their  way  of  estimating  p  is 
different.  In  fact,  MZ's  values  are  closer  to  Model  B  than  A. 

To  distinguish  the  quantitative  nature  of  different  theoretical  models, 
we  note  that  there  are  also  important  factors  that  may  mask  the  comparison 
between  theory  and  experiment  for  R.  One  important  factor  is  that  the 
mixing  enthalpies  extracted  from  phase  diagram  analysis  are  sensitive  to 
sample  and  experimental  conditions.  These  contain  contributions  from 
various  nonideal  structures  such  as  vacancies,  impurities,  dislocations, 
grain  boundaries,  and  surface  conditions,  in  addition  to  the  ideal  AH 

m 

considered  here  for  solid  solutions.  Thus,  our  theoretical  AH  should 

m 

represent  a  lower  bound.  Another  complication  comes  from  the  version  of  the 

theory  of  solid  solution  adopted.  The  theory  used  for  analysis  so  far 

assumes  a  regular  solid  solution  with  second-neighbor  pair  interactions  as 

20-22  21 

was  outlined  in  Section  V.  Recent  experimental  and  theoretical 


studies  have  suggested  the  possibility  of  compositional  clustering  or  long- 

range  correlations  in  alloys.  Extending  the  theory  to  include  such  effects 

will  alter  the  simple  results  for  AH  from  Eq.C^ll).  Moreover,  there  is 

oil 

evidence  from  the  composition  variation  of  the  alloy  hardness  and  from  the 

25 

optical  phonon  frequencies  that  the  shear  coefficients  of  alloys  increase 
near  the  center  of  the  composition  range.  This  will  cause  the  effective 
continuum  shear  coefficient  C  in  Eq.(l4)  to  be  composition  dependent,  which 
will  cause  £j  to  increase.  Despite  these  uncertainties,  useful  comparisons 
across  the  board  in  Table  5  can  still  be  made. 

Based  on  the  above  considerations,  we  can  conclude  that  Models  B,  E, 

and  MZ  predict  n  values  that  are  too  high.  We  should  emphasize  that  all  the 

n  numbers  for  models  from  A  through  MZ  are  directly  calculated  without  any 

adjustable  parameters.  The  fact  that  Models  A,  C,  and  D1 ,  D2,  and  the  FPT 

agree  with  the  experiment  as  well  as  or  even  slightly  better  than  the  one- 
1618 

parameter  theories  ’  (DL  and  FM)  is  already  quite  an  accomplishment.  The 

1 7 

few  numbers  taken  from  Van  Vechtan’s  calculations  (VV)  indicate  that  the 
dielectric  model  predicts  results  at  larger  varience  with  experiments. 

There  are  two  important  implications  about  the  FPT  that  can  be  drawn  from 
Table  3-  First,  the  theory  predicts  a  small  but  negative  J5  value  for 
several  alloys.  This  not  only  means  that  there  is  no  miscibility  gap  in 
these  alloys  but  also  Implies  a  tendency  toward  ordering,  in  which  the 
substitutional  atoms  tend  to  be  surrounded  by  different  second-neighbor 
species.  For  stoichiometric  compositions,  this  implies  a  tendency  toward 
compound  formation.  Secondly,  the  FPT  tends  to  predict  smaller  £5  values 
than  observed  experimentally,  which  should  be  expected  according  to  our 
discussion.  To  the  extent  that  the  FPT  predicts  the  correct  AH  values  for 


6XD 

an  ideal  solution,  the  difference  aH  -AH  may  be  attributed  to  imperfect 

mm 

conditions  and  deviations  from  the  ideal  solution  theory. 

Finally,  the  calculated  £1  values  in  Table  5  provide  some  guidance  in 

2  6  27 

separating  the  completely  miscible  alloys  from  immiscible  ones.  ’  '  In  a 

28 

true  random  alloy,  the  criteria^  for  alloy  mixing  at  a  temperature  T  is 

that  T  ST  ,  where  the  critical  temperature  T  is  given  by  n/(2R  )  with  R 
c  egg 
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being  the  universal  gas  constant.  For  an  alloy  to  be  miscible 

throughout  the  whole  concentration  range,  the  requirement  is  that  both  the 

melting  temperatures  and  of  the  pure  AC  and  BC  compounds  are  greater 

than  T  .  Table  6  lists  the  values  of  T  associated  with  the  £J  values  in  the 
c  c 

FPT,  the  ratios  T  /T.  and  T  /T^,  and  the  average  absolute  values  of  6  for 
c  1  c  2  o 

the  alloys  considered  in  Table  5.  In  Table  6,  T^  is  set  equal  to  zero  if  q 

is  negative  and  is  chosen  to  be  the  lower  value  of  the  two  melting 

temperatures,  so  the  criterion  for  not  having  a  miscibility  gap  is  T^/T2<  1. 

There  is  an  empirical  rule‘s  stating  that  this  will  be  satisfied  if  the 

lattice  mismatch  |  between  the  two  alloy  components  is  less  than  7.5 

percent.  However,  we  find  that  (see  Appendix  C)  a  more  precise  rule  is  that 

Is  I  <  6  ,  where  5  =  1.63  v  and  y  is  the  ratio  of  the  rms  bond  length 
o  m  m  *m  *m 

amplitude  fluctuation  to  the  bond  length  at  the  smelting  temperature  T^.  The 

30 

values  of  T  for  the  compounds  involved  and  the  associate  v  values 
m  m 

estimated  from  Eqs.(C2)  and  (C3)  are  tabulated  in  Table  1.  The  model  used 

2 

in  Appendix  C  yields  T^/T^  =  suggests  that  it  is  instructive 

to  plot  T  /T  as  a  function  of  |6  I  /6  ,  as  is  shown  in  Figure  4  for  the 
c  2  '  o  '  m 

T  calculated  from  FPT.  This  plot  is  similar  to  the  AE  vs  6  curve  in  Fig. 
c  o 

2,  because,  in  fact,  n  is  proportional  to  the  sum  of  the  AE  values  of  the 
two  constituents  [see  Eq.  (Hi)].  However,  if  T^/T^  is  plotted  against 


(t  V 
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Table  6.  Comparison  of  the  Critical  Temperature  T^  of 

Mixing  and  Melting  Temperatures  of  the  Constituents  T^ 
and  T2»  in  the  Order  of  Their  Appearance  in  the 
Parentheses.  Also  Tabulated  are  the  Averaged 
Absolute  Values  of  |6o|  and  the  Ratio  |(So|/6ni 


System 

I6nl  (%) 

0 

U 

H 

T^/T, 

Tc/T? 

1  6o  l/6m 

(A1 ,Ga) As 

0.1 

0 

0 

0 

0.009 

(Al,Ga)P 

0.3 

0 

0 

0 

0.029 

(Al, Ga)Sb 

0.6 

0 

0 

0 

0.067 

(Al,In)Sb 

5.5 

342 

0.25 

0.42 

0.679 

(Ga,In)Sb 

5.8 

204 

0.21 

0.25 

0.716 

(Al , In)As 

6.8 

547 

0.29 

0.45 

0.687 

(Ga, In)As 

6.9 

403 

0.23 

0.34 

0.697 

(Al,In)P 

7.1 

642 

0.36 

0.48 

0.732 

(Ga,In)P 

7.4 

639 

0.43 

0.48 

0.763 

In(P,As) 

3.2 

144 

0.11 

0.12 

0.330 

Al (As.P) 

3.5 

191 

0.10 

0.11 

0.307 

Ga(As,P) 

3.7 

236 

0.14 

0.15 

0.352 

In(As,Sb) 

6.8 

635 

0.52 

0.79 

0.840 

Ga (As » Sb) 

7.6 

924 

0.53 

0.94 

0  844 

Al (As.Sb) 

8.1 

1030 

0.56 

0.78 

0.810 

In(P,Sb) 

9.9 

1450 

1.08 

1.82 

1.222 

Ga(P,Sb) 

11.3 

2180 

1.25 

2.21 

1.256 

AKP.Sb) 

11.6 

2095 

1.19 

1.58 

1.116 

(Cd,Hg)Te 

0.3 

0 

0 

0 

0.033 

(Zn,Hg)Te 

6.0 

377 

0.25 

0.40 

0.659 

(Zn,Cd)Te 

6.2 

312 

0.21 

0.23 

0.564 

Zn(S,Se) 

4.7 

226 

0.11 

0.13 

0.362 

Zn(Se,Te) 

7.2 

569 

0.32 

0.38 

0.615 

Zn(S,Te) 

11.9 

1560 

0.74 

1.03 

1  .01  7 

THEOnV  <b)  EXPERIMENT 


alone,  the  FPT  points  are  much  more  scattered,  and  those  of  SSHS  would  not 
even  exhibit  a  smooth  simple  quadratic  form  because  the  lower  melting 
temperature  is  not  a  smooth  function  of  I  6^  I  •  This  result  suggests 
that  I  6  1/6  <  1  is  a  better  criterion  than  |  6  |<  0.075. 

Figure  3  also  clearly  shows  the  chemical  effects:  all  the  cation- 
substitution  alloy  points  lie  below  the  solid  curve  and  all  the  anion- 
substitution  alloys  have  (T^/T^)  values  on  or  above  the  curve,  corresponding 
to  negative  and  positive  shifts  in  AE  due  to  the  chemical  terms.  Again,  the 
curve  based  the  SSHS  model  is  an  excellent  universal  representation.  From 
the  figure,  we  see  that  all  (P,Sb)  alloys  should  have  miscibility  gaps  and 
all  (As,  Sb)  alloys  are  predicted  to  be  miscible,  although  on  the  border¬ 
line,  because  the  actual  mixing  enthalpies  are  larger  than  these  ideal 
calculated  values.  The  figure  also  shows  that  Zn(S,  Te)  has  a  miscibility 

gap  but  a  smaller  value  of  T  /T-  than  the  (P,  Sb)  alloys,  despite  the  fact 

c  & 

that  its  lA^lvalue  is  larger.  All  these  predictions  are  consistent  with  the 
available  experimental  evidence. 


VII.  SUMMARY  AND  CONCLUSION 

In  this  paper,  a  simple  theory  of  defect  substitution  energies  is 
formulated.  The  substitution  energy  is  compactly  separated  into  a 
replacement  energy  Ap^  and  a  distortion  energy  of  the  pure  host  crystal  [see 
Eq.(2)].  However,  a  rigorous  application  of  this  theory  requires  an 
improvement  in  certain  quantitative  aspects  of  Harrison's  bonding 
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theory,  particularly  the  elastic  constants.  The  most  interesting 
application  of  this  theory  presented  in  this  paper  is  its  perturbation  form 
which  enables  us  to  absorb  the  measured  elastic  coefficients  into  the 
calculation,  and  more  importantly,  to  study  the  chemical  effects.  The 
origin  of  chemical  influences  on  impurity  bond  relaxation  can  be  attributed 
to  three  mechanisms  [see  Eq.(17)]:  a  chemical  force  that  either  helps 
or  hinders  lattice  relaxation,  depending  on  whether  it  has  the  same  or 
opposite  sign  from  the  bond  length  difference  d-d^  between  the  host  and 
impurity,  a  chemical  energy  that  depends  on  the  difference  of  the  polarities 
between  the  impurity  and  host  bonds,  AV^,  and  an  effective  elastic  force 
constant  H  that,  when  positive,  also  tends  to  restrain  the  lattice  from 
distortion.  To  study  the  effect  of  boundaries  between  the  core  atoms  around 
an  Impurity  and  the  rest  of  the  elastic  medium,  various  models  based  on  the 
valence  force  field^*^  are  derived  and  their  results  are  compared  with  the 
full  perturbation  theory  and  available  exerimental  data.  We  found  at  least 
five  models,  including  the  FPT,  that  produce  the  correct  impurity  bond 
lengths  with  variances  for  the  compounds  studied  about  equal  to  the 
experimental  uncertainties  in  EXAFS  ’  ’  (~  0.01  A).  However,  some  modeis 

are  oversimplified  and  will  certainly  not  predict  other  properties  equally 
well.  However,  more  experimental  lattice  constant  measurements,  to  further 
test  the  theory,  particularly  on  (As,  Sb)  and  (P,Sb)  substitution  systems 
for  which  there  are  larger  differences  between  different  models,  are  needed. 
It  would  also  be  instructive  to  see  if  the  predicted  reversal  for  HgCdTe  is 
found. 

The  excess  energies  of  impurity  substitution  are  also  shown  to  provide 
good  estimates  of  the  mixing  enthalpies  q  of  pseudo-binary  alloys.  The 


chemical  shifts  are  found  to  have  a  negative  net  contribution  to  n  for  most 
cation  substitutions,  but  positive  contributions  for  anion  substitutions. 

The  chemical  reduction  of  in  (Ga,  In)  alloys  is  larger  than  1  kcal/mole 
(30  percent  to  100  percent).  Several  VFF  models  and  the  full  perturbation 
theory  produce  results  for  fl  that  are  as  good  as  the  best  theories  with  one 
adjustable  parameter.  However,  the  full  theory  tends  to  yield  answers  on 
the  low  side  of  the  experimental  values,  which  we  argue  is  as  it  should  be, 
because  there  are  nonideal  structures  that  also  contribute  to  n.  The 
calculated  Q  values  and  the  melting  temperatures  are  used  to  predict  the 
existence  of  alloy  miscibility  gaps,  and  the  results  correlate  well  with 
experiments. 

Finally,  we  wish  to  comment  on  the  accuracy  of  the  theories  that  are 
connected  to  the  present  model.  The  perturbation  theory  has  already  been 
stretched  beyond  its  expected  region  of  validity  and  predicts  d^  to  within 
experimental  uncertainties  C  0.01  A)  even  for  cases  with  large  bond  length 
differences  (6^  “  0.1).  The  accuracy  can  only  be  improved  if  the  full  non¬ 
perturbation  theory  outlined  in  Section  II  is  used.  This  calculation  is 
needed  for  the  strong  substitution  cases  that  were  not  considered  in  this 
paper:  examples  are  (B,  Ga),  (B,  In),  (N,P),  (N,As),  (N,Sb),  (0,S),  (0,Se) 
and  (0,Te)  substitutions.  Although  we  believe  that,  for  the  properties 
treated,  the  model  with  a  continuum  attached  to  the  second-shell  is  as 
accurate  as  the  perturbation  theory  used,  it  remains  to  be  seen  if  this  is 

true  for  other  properties,  especially  strain  coefficients.  Finally,  the 

B  6 

present  theory  has  been  extended  to  study  alloys  by  embedding  clusters  in 
an  effective  medium.  This  enables  us  to  study  the  bond  length  and  energy 
variations  throughout  the  whole  concentration  range.  However,  a 
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quantitative  calculation  still  awaits  an  improvement  of  the  accuracy  of 
Harrison's  theory.  A  similar  procedure  is  also  being  extended  to  a  study  of 


the  alloy  electronic  structure,  for  which  a  cluster  CPA  (coherent-potential 

*'  31 

approximation)  involving  both  potential  and  structural  disorder  will  be 

used. 
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Appendix  A 


ELASTIC  ENERGY  IN  CONTINUUM 


In  Section  3.  the  elastic  energy  outside  a  sphere  of  radius  R  centered 

at  the  impurity  is  assumed  to  be  a  continuum  with  a  radial  displacement 

■*''  2  22'' 
uar/r.  If  the  displacement  at  R  is  u^,  then  u  (r)  »  u^  (R  /r  )r.  The 

energy  density  in  the  continuum  is  given  by 


^1  2  2  2 

6c(r)  =  -  C,,(e  +e  +e  )  +  C,_(e  e  +e  e  +e  e  ) 
2  11  XX  yy  zz  12  yy  zz  xx  yy  xx  zz 


1  ?  ?  ? 

+  :^C,,,,(e  +e  +e  ), 
2  *14  xy  yz  zx  * 


where 


3u 


^x  -  ir  =  R\(r2-3x2)/r5. 


'yv  S  y 


2  5 

+  =  -6R  u  xy/r  , 


Appendix  B 


DISTORTION  ENERGY 


In  this  appendix,  we  count  the  detailed  contributions  of  the  bond 
stretching  terms  A(rj^*rj^),  and  "bond  bending"  terms  A  (r^^'r^)  for  i*j  in  VFF 
CEq.(13)]  that  enter  Eq.(l6)  in  FPT  and  in  the  VFF  Models  in  Section  III. 


a  and  6  terms  from  the  first-shell  bonds 


The  four  bond  vectors  pointing  away  from  the  central  impurity  according 
to  Fig.  1  are  r^  -  (1-6 , 1-6 , 1-6 )d//T.  -  (1-6,-1-6,-1-6)d//J  . . .  Thus 

A(r^*r^)  -  -26d^  and  A(r^*r2)  "  ~  •|6d^.  The  a-terms  contribute  4*30^ 


j2 \2  /o 


("26d‘^)‘^/8d*^  »  60^6*^ d^,  and  the  g-terms  contribute  6»36j(-  6d)^/8d^  = 

2  2 

Bj6  d  .  If  an  A  atom  is  replaced  by  a  B  atom,  as  was  done  in  FPT,  the  Qj. 
and  Bj  are  replaced  by  a  and  6  respectively. 


a-terms  from  the  2nd-shell  bonds,  8-terms  between  the  first-  and 


secoqd-shell  bonds  and  among  the  2nd-shell  bonds 


For  these  terms  we  need  to  consider  the  four  bond  vectors  pointing  away 
from  C.  They  are  r^  -  (-1+6,-1+6,-1+6)d//T,  r^  =  (-l+iS.  1+6  +  Y,  1+6  +  Y)  d/ZJ, 
r  -  (1>6+Y,  -1+6,  1+Y+6)d//T,  ...  Then  A  =  |  (6+2Y)d^, 


(6+Y)d^,  and  A(r2*r2)  ”  ■|6d^.  Thus  the  o-terms  from  the 

second-shell  bonds  become  i(*3*3oi  [^(6+2Y)  d^]^/8d^  =  2a  (6+26)^d^,  the  0- 

2  2  2  2 

terra  between  the  first-  and  second-shell  bonds  are  ^*3*36  [j(6+Y)d  ]  /8d  = 
2  2 

2(y+6)  0d  ,  and  the  g-terms  among  the  second-shell  bonds  are  ^•3*36 

2  2  2  2  2  2 
6d  )  /8d  =  26^ed  . 


g-terras  for  the  third-shell  bonds,  8-terms  between  the  second-  and 
third-shell  bonds  and  among  the  third-shell  bonds  adjacent  to  the 
second-shell  atoms 


For  these  terms  we  need  to  consider  the  bond  vectors  pointing  away  from 

B  in  Figure  1.  They  are  r^  •=  (1-6,  -1-Y-5 ,-1-Y-6 )d//J,  r^  ■ 

d+Y',  1+3Y'-Y,  1+3Y'-Y)d//3,  r^  =  (-1-Y",  1+3Y"-Y,  -1 +Y"-Y)d//3' and  r^^  = 

2 

(-1-Y",  -1+Y"-Y,  1+3Y"  -Y)  d//I.  Thus  we  have  A  (r^ •  ^2^  ~  ~  Y  ^36+5Y'), 

2  2 
A(r2'^3)  “  I"  (6  +  2Y-5Y").  ACr^'r^^)  »  -  |  d^Y" ,  A(r^*r2)  “  1“  (■2T-Y'+  3Y"), 

2  2 

^^^3*^3^  "3  andA(r^*r^)  =  •!-  (I'lY'-'IY).  For  Model  A,  Y'  =  r"  = 

2 

0,  so  the  a  terras  of  the  third  shell  bonds  become  4-3»3a  [2(Ar2*r2)  + 

2  2  2  2 

(Ar^«rp  ]/8d  -  8aY  d,  the  B  terms  between  the  2nd-  and  3r’d-shell  bonds  are 

'*-3-3BC(A  ^2’^  3)^]/8d^  -  I  Bd^[6^+|  (6+2Y)^]  =| 

2  2  2  2 

66  d  +6(6+2Y)  d  ,  and  the  6  terras  among  the  third-shell  bonds  adjacent  to 
the  second-shell  atoms  are  4«3«3  6[2(Ar^ •r3)^/(8d^)+(Ar^*rj^)^]=4BY^d^.  For 


continuumm  and  the  only  contribution  from  this  group  are  the  8  terms  between 
the  second-  and  third-shell  bonds.  Since  the  displacements  in  the  continuum 


are  proportional  to  1/R‘^,Y'  -  8/2' Y/(19/T9)  and  y"  -  8/2Y/(1  l/TT) .  Th 

these  6  terms  become 


— -  BC|  (36+5Y')^4^  (6--2y-5Y")^] 

2d  ^  ^ 


I  BC(36+ 


j40/2 

19/T9 


Y) 


2 


_  _ _  ^0/?  ,.v2-..2. 

"^2(6  +  2Y-Ti^Y)  Jd 


B  terms  for  bonds  adjacent  to  the  third-shell  atoms 


These  terms  only  enter  Model  A,  so  r’  •  r"  -  0.  There  are  two 
different  groups,  one  like  those  adjacent  to  C  and  another  like  those 
meeting  at  C".  The  four  bond  vectors  pointing  away  from  C  are  r^= 

(-1,  -1+Y,  -1+Y)  d/ZJ,  r^  ”  (-1,1,1)  d//3,  r^  -  (1,-1,1)  d//T,  and  r^  = 
(1,1,-1)  d//5.  Thus  the  only  contribution  from  this  group  is  ^l-3*3B 

[A(r«r2)]^/8d^  «  2BY^d^.  The  four  bond  vectors  around  C"  are:  r^ 

-  (-1+Y,  -1-Y,  -1)d//T,  r^  =  (-1,1,1)d//3,  r^  =  (1,-1+Y,  1+Y)d//3,  and  r 
(1 , 1 ,-1 )d//3,  which  only  results  in  the  first  order  term 
»  j  Yd^.  Thus  the  group  contributes  to  '1*3*3B[2  (^r^‘r^)^]/8d^  «=  ^BY^d^ 


2  2 

and  the  combined  contribution  from  these  two  groups  is  68Y  d  . 


Appendix  C 


Starting  with  Eq.(2i»)  and  using  the  SSHS  model  a  -  Oj ,  one  finds  the 
mixing  enthalpy  parameter  £j  to  be 


a 


where  N  is  Avagadro's  number  and  a  =  (1/2)  (a.^+Oo^).  Then  relate  the  mean 
O  Aw  Ow 

2 

square  bond  length  fluctuation  <z  >  at  the  melting  temperature  to  T^for  a 
compound  by  equating  the  average  potential  energy  per  unit  cell  to  half  of 
the  thermal  energy: 


<P.E.>  s  =  •?(2-3kT„) 

2  2  m 


where  k  is  the  Boltzmann  constant.  Defining  a  Liedermann  ratio  of  melting 
X  by 


2,  J/2 


and  choosing  the  mixing  criterion  to  be  T  /T  <1,  where  T  now  is  the 

c  m  m 

smaller  value  of  the  two  melting  temperatures  of  the  constituent  compounds 
we  require  that 


nk 


2  2 

IJR  ax  d 
g  *m 


3  (‘^AC-^O^ 


2  h2 


<1. 


or 


where  6 


1.63  Xm  ^hd  U  I  is  the  percentage  bond  length  difference. 


REFERENCES 


•  Visiting  professor  at  Stanford  University. 


1.  J.  C.  Mikkelsen  and  J.B.  Boyce,  Phys.  Rev.  Lett.,  49,  1412  (1982); 
Phys.  Rev.  B28,  7130  (1983). 


2.  A.  Zunger  and  J.  E.  Jaffe,  Phys.  Rev.  Lett.  51,  662  (1983). 


3.  J.  L.  Martins  and  A.  Zunger,  Phys.  Rev.  B30,  6217  (1984), 


4.  A.  Balzarotti,  A.  Kisiel,  N.  Motta,  M.  Zimmal-Starnawska,  M.  T.  Czyzyk, 
and  M.  Podgorny,  Phys.  Rev.  B30,  2295  (1984). 


5.  A.-B.  Chen  and  A.  Sher,  "Bonding  Theory  of  Semiconductor  Alloys," 

Microscience  3  (A  multiclient  publication  of  SRI  International,  Menlo 
Park,  California,  1984). 


6.  A.  Sher,  A.-B.  Chen,  and  W.  E.  Spicer,  13th  International  Conference  on 
Defects  in  Semiconductors,  August  12-17,  1984,  Coronado,  California, 

p.43. 


7.  W.  A.  Harrison,  Electronic  Structure  and  Properties  of  Solids  (W.  H. 
Freeman  and  Co.,  San  Francisco,  1980);  "The  Bonding  Properties  of 
Semiconductors,"  Microscience,  3,  p.  35  (A  multiclient  publication  of 


52 


SRI  International,  Menlo  Park,  California,  1983);  and  Phys.  Rev.  B27, 
3592  (1983). 

K.  Shih,  W.  E.  Spicer,  W.  A.  Harrison,  and  A.  Sher,  Phys.  Rev.  B 
(to  appear  in  Jan  15,  1985). 

A.  Baldereschi  and  J.  J.  Hopfield,  Phys.  Rev.  Lett.  28,  171  (1972). 

R.  M.  Martin,  Phys.  Rev.  i»005  (1970). 

The  experimental  elastic  constants  were  taken  from  the  tabulated 
values  in  Refs.  10  and  25  except  AlP  and  AlAs,  for  which  we  used  the 
estimated  values  by  J.  0.  Wiley  in  Semiconductors  and  Semimetals,  Vol. 
10,  edited  by  R.  K.  Wlllardson  and  A.  C.Beer  (Academic  Press,  New  York 
1975),  p.  13^1. 

R.  Kibuchi,  Physica  103B,  4l  (1981). 

A.-B.  Chen  and  A.  Sher,  unpublished. 

J.  B.  Boyce  and  J.  C.  Mikkelsen,  unpublished  data  for  Zn(Se,  Te).  We 
thank  them  for  their  communication  prior  to  publication. 

Extrapolated  results  quoted  in  Ref.  3. 


G.  B.  Stringfellow,  J.  Cryst.  Growth  27,  21  (1974);  58,  194  (1982). 


I 


HD-RICC  999  SEMICONDUCTOR  ALLOV  THEORVIU)  AUBURN  UNIV  ALA  DEPT  OF 
PHVSICS  A  CHEN  27  SEP  89  AF0SR-TR-8S-B199 
AFOSR-84-0282 

UNCLASSIFIED  F/G  28/12 


2/2 


J.  A.  Van  Vechten,  in  Semiconductor  Handbook,  Vol.  3.  edited  by  S.  P. 
Keller  (North-Holland,  Amsterdam,  1980),  p.  1 


P.  A.  Fedders  and  M.  W.  Muller,  J.  Phys.  Chem.  Solids  ^5,  685  (19811). 


M.  B.  Panish  and  Ilegems,  Prog.  Solid  State  Chem.  7,  39  (1972);  A. 
Langier,  Rev.  Phys.  Appl.  8,  2959  (1973);  and  Ref.  16. 


T.  P.  Pearshall,  R.  Carles,  and  J.  C.  Portal,  Appl.  Phys.  Lett.  £2,  1136 
(1983). 


K.  Kakimoto  and  T.  Katoda,  Appl.  Phys.  Lett.  1|0,  826  (1982). 


P.  Parayanthal  and  F.  H.  Poliak,  Phys.  Rev.  Lett.  1922  (19811), 


G.  P.  Strivastova,  J.  L.  Martins,  and  A.  Zunger,  Phys.  Rev.  B  (in 
press) . 


N.  A.  Groyunova,  A.  S.  Barshchvevskii ,  and  D.  N.  Fretiakov,  in 
Semiconductors  and  Semimetals,  Vol.  il.  Edited  by  R.  K.  Wlllardson  and 


A.  C.  Reer  (Academic  Press,  New  York,  1968),  Ch.  1. 


S.  S.  Mitra  and  N.  E.  Massa,  in  Handbook  on  Semiconductors,  Vol.  1 
edited  by  W.  Paul  (North-Holland,  1982),  p.  29. 


L.  M.  Foster,  J.  Electrochem.  Soc.  Solid  State  Sci.  and  Techn.  _1^ 

1662  (19711). 


L.  M.  Foster  and  J.  F.  Woods,  J.  Electrochem.  Soc.  Solid-State  Sci.  and 
Techn.  1919,  50^  (1972). 


R.  Kubo,  Statlstloal  Mechanics  (Univ.  of  Tokyo,  1971),  p.  336. 

The  critical  temperature  will  be  reduced  by  a  few  percent  If  more 
elaborate  models  are  used.  See  .  W.  Christian,  Transformations  In 
Metals  and  Alloys,  2nd.  edition,  edited  by  D.  W.  Hopkins  (Pergamon 
Press,  New  York,  1981),  p.  19^. 

American  Institute  of  Physics  Handbook.  3rd  ed.  (McGraw-Hill,  New 
York,  1972),  p.  9-16. 

K.  C.  Hass,  R.  J.  Lempert,  and  H.  Ehrenrelch,  Phys.  Rev.  Lett.  52,  77 


